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ABSTRACT 

FORAST  is  a  programming  language  designed  for  use  on  ORDVAC  and  BRLESC, 
the  high-speed  digital  computers  of  the  Ballistic  Research  Lahoratories . 
Programs  written  in  this  language,  with  minor  limitations,  may  be  executed 
on  either  computer.  BRL  Report  No.  1172,  tl]  #  describes  FORAST  in  its 
generality  and  was  written  primarily  for  professional  programmers.  This 
report  is  intended  for  the  novice.  Fundamental  concepts  and  details  of  the 
language  are  illustrated  in  many  examples  so  that  the  novice  is  taught  how  to 
program  and  obtain  practical  solutions  for  a  variety  of  mathematical  problems. 
Intended  as  a  supplement  to  [l.]  ,  this  report  does  not  illustrate  the  full 
generality  of  the  language.  Some  of  the  material  is  repetitious  but  amplified 
and  several  references  are  made  to  d]  . 
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I.  INTRODUCTION 


The  Ballistic  Research  Laboratories*  high-speed  computers^  ORDVAC  and 
BRLESC,  are  used  to  obtain  solutions  of  mathematical  problems.  In  general, 
the  numerical  solutions  are  obtained  by  instructing  the  high-speed  computers 
to  EVALUATE  FUNCTIONS  that  approximate  the  solutions  desired.  Hence,  to  obtain 
solutions  from  high-speed  computers,  the  computers  are  instructed  to  carry  out 
the  detailed  operations  required  to  evaluate  particular  functions  of  interest. 


Conventional  mathematical  notation,  which  permits  considerable  freedom 
of  expression  and  in  which  problems,  solutions,  and  functions  are  generally 
expressed  in  a  form  governed  by  human  convenience,  is  not  fully  acceptable  by 
high-speed  computers.  Indeed,  the  instructions  to  computers  that  specify  the 
essential  operations  and  the  numerical  quantities  involved  in  the  operations 
are  ultimately  expressed  in  a  PRIMITIVE  LANGUA.GE  that  is  understood  by  the 
computer.  For  most  high-speed  digital  computers,  the  primitive  language 
consists  of  combinations  of  only  two  characters,  0  and  1,  the  basic  characters 
of  the  binary  number  system.  Specific  combinations  of  these  characters  must 
be  constructed  to  represent  the  sequence  of  operations  and  the  quantities 
required  to  obtain  a  numerical  solution.  The  human  task  of  constructing  the 
required  combinations  of  binary  characters  is  tedious,  difficult  and  highly 
susceptible  to  human  error.  To  simplify  this  task,  HIGHER-LEVEL  LANGUAGES  are 
designed  so  that  computers  can  be  instructed  by  professionals  with  programs  in 
their  primitive  languages  to  automatically  translate  them  into  their  respective 
primitive  languages.  In  contrast  to  the  PRIMITIVE  LANGUAGE,  the  HIGHER-LEVEL 


LANGUAGES  admit  combinations  of  the  letters  of  the  alphabet,  the  decimal  digits, 
and  some  mathematical  symbols  such  as  +  -  =  .  etc.  The  programs  which 

automatically  accept  the  higher-level  language  and  carry  out  the  translation 
are  caxxeu  uui''irxxiCX\D,  aSaEMB i ,  xxvHMoijrtiurvo,  eoc.,  and  geiiexcixxy  aumxo  bume 
of  the  freedom  of  mathematical  notation  together  with  some  English  word  state¬ 
ments  for  controlling  sequences  of  operations.  Although  the  higher-level 


languages  do  indeed  relieve  the  human  of  uiaiiy  tedious  details ,  a  novice  may 
still  feel  burdened  with  the  necessity  for  specifying  particular  details  re¬ 
quired  in  a  high-level  language. 
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F0RA3T  is  a  high-level  language  designed  lOr  use  on  QRDVAC  and  BRLE3C. 
BRL  Report  No.  1172,  [1]  ,  describing  FORAST  in  its  generality  was  intended 
as  a  reference  manual  for  professional  programmers.  This  report  is  intended 

■vriY^*^  e  a  n  ■?  q  vi+  c 

L>0  dXiVU.  _ !_»_> 


for  the  novice  and  is  confined 
illustrated  by  numerous  examples,  teaching  the  novice  how  to  program  a  variety 
of  mathematical  problems.  Each  example  includes: 


a  statement  of  a  problem; 
a  method  for  obtaining  a  solution; 


CL 

a 

a 


chart  indicating  the  seQ_U0 


nee  of  the  necessary  operations; 


program  written  in  FORAST; 

tabulated  copy  of  the  results  obtained  from  the  computer. 


Many  of  the  methods  used  in  the  examples  were  chosen  to  illustrate  particular 
concepts  of  the  language  and  do  not  necessarily  represent  the  most  efficient 
solutions  in  the  given  situations. 


r\ 


o 


II.  EVALUATION  OF  FUNCTIONS 


Solutions  to  many  mathematical 
recording  of  mathematical  functions, 
operations : 


problems  can  be  reduced  to  evaluation  and 
Below  are  listed' the  four  associated 


1. )  SUBSTITUTION  of  numerical  values  for  each  argument  of  the  functions; 

2. )  FERFOPmMING  the  j.n(ij.CcL"u0  d  arithmeti  c  opox’o.'tions j 

3. )  RECORDING  numerical  values  of  relevant  functions; 

4. )  CONTROLLING  the  sequence  of  operations  involved  in  1,2,  and  3- 


More  specifically,  to  obtain  a  solution  to  a  problem  from  a  high-speed  computer, 
we  instruct  the  computer  to  substitute,  to  evaluate,  to  record  and  to  control 
all  operations  required  to  obtain  a  solution. 


In  a  given  problem,  all  these  operations  must  be  planned  in  advance  and 
the  computer  must  receive  specific  instructions  as  indicated  on  a  Flow  Chart. 
A  general  Flow  Chart  is  illustrated  schematically  in  FIGURE  1. 


FIGURE  1. 


GENERAL  FLOW  GH.ART  FOR  EV.4LUATING  FUNCTIONS 
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Before  describing  the  fundamental  concepts  of  the  language^  we  consider 
next  the  major  components  of  the  computers  and  the  function  of  each  as  regards 
their  use  in  obtaining  a  solution  of  a  problem. 


FIGURE  2. 

IVIAJOR  COMPOh^riTS  OF  HIGH-SPEED  DIGITAL  COMPUTERS 


The  components  shown  in  FIGURE  2.  are  typical  for  most  high-speed  digital 


Ild-lIJ.  V. 


of 


computers.  The  in^ut  devices  generally  consist  of  electro— mec 
reading  paper  cards^  magnetic  or  paper  tape.*  The  output  devices  use  similar 


*  To  record  information  on  paper  cards  or  paper  tape^  the  cards  or  tapes  are 
perforated^  (i.e.^  holes  are  punched  in  them^,  wherein  specific  combinations 
of  punches  define  specific  characters.  On  magnetic  tape,  specific  combi¬ 
nations  of  small  magnetized  areas  correspond  to  specific  characters . 
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mrt  <~i  v-i  ^  V 

iiiCCLlXO  UU  _^XXiiL»  U1 


r»/^T»TrM  i+qv^c?  llc?^  r*  /^TnV\  •^■»nQ+'-{/^v^r'  /•n'P 

viCLiiy  j.  o  i^x  wiio  v>^J. 


these  devices  together  with  high-speed  printers  and  other  optical  or  photo¬ 
graphic  devices  for  output.  The  input-output  devices  may  also  serve  as 
external  storage  units.  Magnetic  tape  units  may  and  are  used  as  input  or  out¬ 
put  devices  or  even  as  external  storage.  Hence  input  and  output  devices  are 
used  to  get  information  into  and  out  of  the  internal  storage  as  indicated  by 
the  arrows  of  FIGURE  2.  In  general,  the  information  that  passes  through  the 
inpu't  U.0V1C0S  is  r0coTu.0u.  in  in‘t0rn3,l  soOj/3,g0  s-nd  't.u.0  inf* o xijia't  1  on  ihEi  is  x*0- 
corded  on  output  devices  is  a  copy  of  selected  information  that  exists  in 
internal . storage.  Information  recorded  in  internal  storage  is  retained  there 

+  +V»^  +V»q  r»  /^wrM  i+^v* 

CL  O  CLO  CLC.  iD  J.  J.  C;  Cl.  ^  C:CI.  J.O  -LO  UCLO-Xi^  VA  J.  ^  LAli  C/  J. -C  U  lilC  C.  U ICX  J.O  J.ilOl.'A  LAC.  ^ 

ed  to  replace  it  with  other  information.  The  arithmetic  unit  is  used  to  perform 
ordinary  arithmetic,  i.e.,  addition,  subtraction,  multiplication  and  division. 


The 


I  T  +  * 


Q^+T'xro+C 


.q1 1 


in 


instructions  that  have  been  recorded  in  internal  storage. 


Knowing  the  capabilities  of  the  major  units  and  the  fundamental  concepts 
of  the  high-level  language,  a  programmer  designs  a  set  of  instructions  (expressed 
in  the  high-level  language)  that  will  when  executed  by  the  computer  produce  a 
solution  to  a  given  problem.  The  set  of  instructions  and  accompanying  data 
(called  a  program)  are  recorded  on  paper  cards.  These  cards  are  inserted  in 
a  card  input  device.  Without  loss  of  generality,  one  may  assume  that,  after 
activating  a  few  switches  on  the  control  panel  of  the  computer,  an  equivalent 
copy  of  the  Instructions  that  were  recorded  on  cards,  is  recorded  in  internal 
storage.  The  instructions  are  then  executed  and  res\ilts  produced  accordingly. 
Actually,  the  FORAST  COMPILER  which  is  temporarily  recorded-  in  internal  storage 
activates  the  card  input  device,  and  translates  the  high-level  language  to  the 
primitive  machine  language.  It  is  the  primitive  machine  language  corresponding 
to  the  high-level  language  that  is  recorded  in  internal  storage.  In  interpre¬ 
ting  and  attempting  to  transform  the  high-level  language,  the  compiler  can  and 
does  detect  some  violations  of  grammatical  errors;  in  such  cases, it  records 
pertinent  information  on  the  output  device  that  identifies  the  particular 
violation  and  does  not  permit  the  computer  to  attempt  to  execute  the  program. 
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In  suBjuiary,  corresponding  "to  a  given  problem,  we: 

1. )  plan  a  solution  and  exhibit  the  plan  in  the  form  of  a  flow-chart; 

2.  )  write  a  set  of  instructions  in  the  high-level  language  to  execute 

the  plan; 

3-)  record  the  set  of  instructions  and  accompanying  data  on  paper  cards; 

4  1  eTiV»Tn-?+  +V»a  >“>0 1 • 

5.)  obtain  the  computer  output  and  tabxilate  the  results. 

Although  not  explicitly  stated  previously,  it  is  important  to  recognize  that 
once  a  program  has  been  tried,  tested  and  proved,  it  can  be  re-used  for  other 
sets  of  numerical  data. 
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EXAMPLE  1. 


To  illustrate  some  of  the  general  concepts  stated  thus  far,  and  to 
Introduce  some  fundamentals  of  the  language  we  consider  next  Example  1. 


X.  ,  y.  ,  where  i  =  1,2,3, 
1  '  1 


i.e.,  a  table  of  numerical  values  for 


X 

1 

30 

n 

-o 


Y 

1 

40 

6 


REQUIRED; 


For  each  pair  of  values,  X^,  Y^,  compute  and  record 


i.e.,  we  want  to  produce  a  table  with  entries 

X  Y  Z(X,Y) 

1  1  1.4l4 

30  4o  50 

-8  6  10 


for  an  indefinite  number  of  pairs,  (X,Y). 
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To  plan  a  solution  for  this  problem,  we  tacitly  assume  that  the  given 
numerical  values  of  X  and  Y  will  be  recorded  in  a  "prescribed  form"  on  paper 
cards,  one  pair  of  values  on  each  card,  and  utilizing  as  many  cards  as  are 
necessary  for  the  complete  table  of  values •  Although  this  is  not  the  most 
efficient  method  for  recording  the  data,  it  is  a  practical  assumption  that  is 
readily  fulfilled.  One  could  assume  two,  three,  or  more  pairs  recorded  on  each 


card;  howe 


T  T*  r»Q 

V.4.  I 


Next,  we  make  use  of  a  fundamental  characteristic  of  card-input  devices. 
Cards  containing  information  to  be  used  in  the  solution  of  a  problem  are 
stacked  in  the  input -hopper  of  the  input-device  with  the  card  at  the  bottom  of 
the  stack  "engaged"  under  a  "read"  station.  (See  FIGURE  7-)  When  the  input- 
device  is  activated  by  the  computer,  the  information  punched  in  the  card  under 
the  "read'*  station  is  interpreted  and  the  information  transmitted  to  internal 
storage.  During  the  process  of  interpretation  and  transmission,  the  card  auto¬ 
matically  moves  to  a  storage  bin  and  the  "next"  card  in  the  stack  moves  under 
the  read  station  ready  for  the  next  activation  of  the  input  device.  Hence, 
cards  that  have  been  read  accumulate  in  the  storage  bin  and  are  not  accessible 
to  the  computer  for  re-reading  unless  manually  re-inserted  in  the  input-hopper. 
Cards  are  generally  read  only  once  for  a  given  problem  since  the  information 
recorded  on  them  is  recorded  In  internal  storage  indefinitely  for  as  many 
future  references  as  desired.  A  similar  sequence  of  card-processing  applies 
to  card-output  devices,  i.e.,  the  automatic  passage  of  cards  from  input-hopper 
to  "punch"  station  and  then  to  a-  storage  bin.  To  utilize  this  fundamental 
sequence,  we  plan  to  have  a  pair  of  values  read  and  substituted  for  X  and  Y, 
instruct  the  computer  to  evaluate  the  corresponding  function,  record  the  result 
on  the  card  output  device  and  then  repeat  the  process  using  the  next  pair,  etc. 
Hence  our  plan  can  be  exhibited  in  the  form  of  the  simple  Flow-Chart  as  shown 
in  FIGURE  5. 
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"WORDY"  CHART 


CONCISE  SYMBOLIC  CHART 


l.iT  /^TTT^TTl  'T 

rxuuruli 

A  SIMPLE  PLAN  TO  OBTAIN  A  SOLUTION  FCR  MAMPLE  1. 
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The  "vordy"  chart  is  intentionally  "wordy"  for  the  beginner.  The  concise 
symbolic  chart  is  generally  used  and  implies  all  that  is  stated  explicitly  in 
the  "wordy"  chart.  Note  the  capitalized  words:  SUBSTITUTE  and  its  equivalent^ 


READ:  EVALUATE  and  the  eauivalent  exollcit  definition. 


=  7x2 


+  Y 


RECORD 


/ I _ _ • 1 L.  T  •RTm  f  T  TMTO  ^  J  T  niT/^TT _  _  \ 

^^ex'uerna.j.j.y ;  a.na  ii/y  eq,u±va,j_eriu^  rnxi^x  .  v-^nxi^x  ciiia  ctre  bynunumuus .  ; 

Note  also  that  both  charts  indicate  that,  after  recording  results  correspond¬ 
ing  to  a  given  pair  of  values,  the  next  processing  desired  is  indicated  by  a 

+.*.  +Vi/^  A  t-tV.  cs  v+  -I'v.  +..\  "Wrt  — iv..*. 

Lr>.^  LrlH^  X  i  ^  J.  1 ICL  _l.  XJXXXllX  V^CXU-OCO  OliC  llC.A.Lr  _pa,  J.  i  Lr<^  “ 


cessed  in  the  same  manner  as  its  predecessor.  The  process  is  to  be  continued 
until  the  cards  in  the  input-device  are  exhausted. 


Next,  we  write  instructions  in  the  FORAST  language  to  instruct  the  computer 
to  carry  out  the  desired  processes  shown  on  the  flow-chart.  These  instructions 
are  written  on  the  standard  CODING  FOP.M  shown  in  FIGURE  4. 
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FCRAST  CODING  FORM 


PROBLEM  C'906 _ 

CODER  M.J.RomanellDATE  1  Jul  65  PAGE  1 


on  tyh©  Cod.m§  Pojrni  m*©  punoliGu.  on  cni'u.s  whicii 
can  then  be  submitted  to  the  computer.  The  cards  corresponding  to  the  in¬ 
structions  written  on  the  Coding  Form  of  FIGURE  4.  are  shown  in  FIGURE  5. 
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The  instruction  cards  of  FIGURE  5>  together  with  a  few  data  cards,  (the 
cards  on  which  we  assumed  the  X*s  and  Y’s  would  be  recorded),  are  shown  in 
FIGURE  6.  These  cards  represent  the  complete  package  to  he  submitted  to  the 
computer  to  obtain  the  desired  solutions 
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FIGURE  6. 


Recall  that  we  assumed  the  given  input  data,  (the  numerical  values  of  the 
X’s  and  Y’s),  woiild  be  recorded  on  paper  cards  in  a  "prescribed  form".  The 
FORAST  programming  language  provides  for  instructing  the  computers  to  accept 
data  represented  on  cards  in  many  forms.  A  few  specific  forms  are  considered 
standard,  (i.e.,  those  forms  that  are  most  convenient  for  scientific  calcu¬ 
lations  and  used  frequently),  and  only  departirres  from  standard  forms  must  be 
c;T.pr>-i-rn pH  in  r  TTiRnnp-n  provided  by  the  language.  The  prescribed  form  chosen 


<>..V 


for  this  example  is  referred  to  as  standard  "floating-point"  form.  This  form 
is  similar  to  that  generally  called  scientific  notation,  where  numerical  values 
are  expressed  in  a  two-part  form;  one  part  called  the  "coefficient",  the  other 
is  called  the  "exponent".  For  example,  the  quantity  J2.h^  can  be  expressed  as: 


7.2U5  X  10  or  .72U5  X  10^  or  72U5.  x  lO"^  ,  etc. 

The  base  ten  is  generally  omitted  and  the  two-part  form  used  consists  of  the 
signed  coefficient  and  the  signed  exponent.  The  exponent  is  restricted  to 
integers.  The  computers  represent  mmibers  internally  in  a  similar  manner,  using 
binary  notation  instead  of  decimal.  In  ORDVAC,  the  base  is  2,  binary;  in  BRLESC, 
the  base  is  I6,  sexadecimal.  The  use  of  other  bases  and  forms  is  permissible  and 
means  for  transforming  from  one  base  or  form  to  another  base  or  form  is  provided. 


For  input  data,  some  standard  floating-point  forms  used  are: 


Exponent 

a. ) 

Sign, 8  decimal  digits 

Sign, 2  decimal  digits 

h. 

Sign, 7  decimal  digits 
and  a  decimal  point 

Sign, 2  decimal  digits 

c.  ) 

Sign, 11  decimal  digits 

(no  exponent) 

a 

"  •  y 

St  crn  "1  0  HpR-impT  Hicri'hc 

and  a  decimal  point 

(no  exponent) 
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In  the  forms  where  the  decimal  point  of  the  coefficient  is  not  specified^ 


(forms  "a"  and  '*c"),  it  is  assumed  to  be  before  the  first  decimal  digit 
punched.  In  each  of  the  above  forms,  12  columns  of  a  card  are  used.  Other 
variations  of  these  forms  are  permitted;  however,  in  EXAMPLE  1.  we  used  form 
"a”  for  both  input  and  output  representations.  As  is  conventional  in  mathe¬ 
matical  notation,  the  positive  sign  may  be  omitted. 

To  introduce  some  of  the  fundamental  rules  of  the  FORAST  language,  we 
-rAfp-r  fVip  T’pafl<=T'  +.n  st.anflaT’fl  end i  no-  fn-rm  shown  "in  FTOTIRE  L. 


Observe  the  division  of  the  form  into  four  general  headings,  namely; 
LOCATION,  (1-6);  ORDER  TYPE,  (7-10);  FORMJLA  STATEMENTS,  (II-76);  and  GAJ5D 
IDEN,  (77-80).  The  nxmibers  written  here  in  parentheses,  and  written  under 
the  general  headings  on  the  coding  form  correspond  to  the  colimms  of  the 


paper  cards  on  'vdiich  the  information  will  be  recorded.  The  space  allocated 
for  eolumns  (I-6),  labelled  LOCATION,  is  reserved  for  names  to  refer  to  state¬ 
ments  and  pertinent  data.  Columns  (7-IO),  labelled  ORDER  TYPE,  are  reserved 
for  classifying  the  information  recorded  in  columns  (1I-76).  Only  a  small 
number  of  standard  names  are  permitted  in  this  space.  A  list  of  such  names 
will  be  furnished  later.  Columns  (II-76)  are  reserved  for  the  formulas,  English 
word  statements,  and  other  data  required  to  describe  the  processes  necessary  to 
obtain  a  desired  solution.  Columns  (76-80),  labelled  CARD  IDEN  are  reserved 
for  arbitrary  identification  of  individioal  cards.  Although  not  mandatory,  the 
cards  are  generally  ordered  in  n\mierical  order  for  easy  reference  so  that  if 
and  when  the  compiler  detects  a  violation  when  it  reads  a  card,  the  identifi¬ 
cation  recorded  in  columns  (76-80)  of  the  card  is  printed  to  identify  the  card 
containing  the  violation. 


The  first  card  of  each  pxogxtuui  submitted  to  the  computer  must  have  PROS 
recorded  in  colimns  (7-I0).  Coliamns  (II-76)  of  this  card  may  be  used  to  record 
a  problem  number,  persons  name,  title  of  problem,  etc.  The  information  recorded 

vjii  wiiiD  ucii  u.  J.&  ctJ-Wd^iD  xc^JWXu.c;u.  <ao  ua.x  u.  wx  one:  v.^u.ojju.o  ovj  xu.c;iioxx^  one 

subsequent  results  which  follow,  if  any! 


on 


To  identify  the  beginning  of  a  program,  an  arbitrary  name  composed  of  at 
most  6  characters  must  be  recorded  in  columns  (l-6),  i.e.,  opposite  the  first 
statement  to  be  executed.  In  the  example,  we  chose  the  name,  START.  Note 
that  commencing  in  column  11  of  this  card  we  wrote  RKAD(X)Y,  the  instruction 
designed  to  start  the  processing  of  a  data  pair  by  copying  in  internal  storage 
the  numerical  values  recorded  on  the  first  data  card.  Next,  we  wrote  the 
function  to  be  evaluated.  No  doubt  the  beginner  will  note  the  difference  in 


the  conventional  mathematical  notation,  Z  =  i/ic'  +  Y‘"  and  the  cumbersome 

notation  Z  =  SQRT(X**2  +  Y**2).  Unfortunately,  the  former  is  not  acceptable 

to  the  input  media.  In  contrast,  the  latter  uses  a  serial  arrangement  of 

acceptable  symbols  or  their  combinations.  Subscripts,  superscripts,  exponents 

or  other  conventional  symbolism  which  is  not  included  in  the  restricted  set  of 

characters  acceptable  by  input  media  are  denoted  by  special  characters  in  the 

restricted  set  or  by  combinations  of  them.  For  example,  in  place  of  the  symbol 

,  we  use  SORT  (  )  or  (  )**.5  enclosing  the  argument  in  parentheses. 

The  latter  expression  indicates  exponentiation,  as  do  the  forms  X**2  and 

2  2 

Y**2,  which  correspond  respectively  to  X  and  Y  .  The  general  form  of  ex¬ 
ponentiation  is 


(E1)**(E2) 


which  represents  the  expression  El  raised  to  the  E2  power.  The  equality 
symbol,  =,  is  used  in  a  different  sense  from  the  customary  mathematical  meaning 
in  that  even  expressions  such  as 

X=X+1  or  U=Z=Y+  Z**2 


are  valid.  The  special  meaning  of  equality  is:  EVALUATE  the  expression  to  the 
right  of  the  equality  symbol,  using  values  of  the  quantities  that  currently 
exist  in  internal  storage,  then  record  the  resulting  value  in  internal  storage 
as  the  existing  value  of  the  function  or  variable  whose  name  (or  names)  appear 
to  the  left  of  the  equality  symbol.  In  the  example,  U  =  Z  =  Y  +  Z**2  the 
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vjuiiijjater  will  first  square  the  existing  value  of  Z,  add  this  result  t< 


existing  value  of  Y  and  then  record  this  resiilt  in  internal  storage  as  the 
existing  value  of  U  and  Z;  the  previous  values  of  Z  and  U  are  erased  just 
■■^rior  to  the  recordin<^  of  the  new  value. 


On  line  k  of  FIGURE  k  appear  two  distinct  statements  with  the  special 
symbol^  used  to  separate  them.  One  may  write  many  statements  or  symbols 
on  a  given  line  in  the  space  provided.  Each  line  of  the  coding  form  is 
pimched  on  a  single  card.  Three  lines  (or  3  cards)  would  have  sufficed  for 
EXAMPLE  1.  as  illustrated  below: 


START 


PROB  0906  1 

READ(X)Y  ^  Z  =  SQRT(X**2  +  Y**2)^  PRINT(X)Y)Z5t  GOTO(START)  2 
END  GOTO (start)  3 


The  rules  for  the  use  of  parenthesis  will  be  explained  later.  The  last  card 
carries  the  word  END  recorded  in  the  ORDER  TYPE  space.  This  card  serves  two 
purposes : 

1. )  It  signifies  to  the  FCRAST  compiler  the  end  of  the  translation  of 

the  FORAST  program  into  primitive  machine  language.  The  equivalent 
of  this  card  is  not  recorded  in  internal  storage. 

2. )  The  GOTO(  )  statement  recorded  on  this  card  directs  the  computer 

to  the  first  statement  of  the  progreun  to  be  executed.  The  first 
statement  to  be  executed  need  not  follow  the  PROB  card,  it  may  be 
any  statement  in  the  program. 

The  GOTO(START)  statement  on  line  2  is  recorded  in  internal  storage  and  directs 
the  computer  to  process  every  data  pair  following  the  first  pair.  The  GOTO(START) 
on  line  3  directs  the  computer  to  commence  the  program  with  the  reading  of  the 
first  data  pair.  It  is  not  generally  true  that  the  processing  of  subsequent 
cases  begins  with  the  first  instruction  of  every  program. 
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FIGURE  7. 


Illustrated  in  FIGURE  J  are  the  machine  components  used  to  obtain  the 
solution  of  oxir  problem.  Note  that  the  cards  shown  in  FIGURE  6  are  first 
inserted  in  the  CARD  READER ^  ( input -device ) .  The  FORAST  compiler  causes  the 
five  (5)  program  cards  to  be  read  and  the  information  on  them  is  recorded  in 
internal  storage  as  illustrated.  When  the  compiler  recognizes  the  card  with 
END  punched  in  columns  7-10;  (the  5th  card),  it  directs  the  control  xmit  to 
commence  executing  the  instructions  that  have  been  recorded  in  internal  storage. 
Note  also  that  in  addition  to  utilizing  space  for  the  instructions  in  internal 
storage,  space  is  also  reserved  in  internal  storage  for  the  nvunerical  values 
of  X,  y,  and  Z.  Each  time  the  computer  carries  out  the  instruction,  RKAD(X)Y 
the  numerical  values  on  the  card  at  the  read  station  are  transmitted  and  re¬ 
corded  in  internal  storage,  the  card  itself  being  moved  to  the  storage  bin. 

Next,  using  the  existing  values  of  X  and  Y  currently  in  internal  storage,  the 


computer  evaluates  Z  X^  +  Y^  and  records  the  numerical  value  thus  obtained 
in  internal  storage.  To  carry  out  the  next  instruction,  PRINT  (X)y)Z  the 
control  unit  activates  the  CARD  PUNCH  and  records  on  a  card  the  numerical  values 
of  X,  Y,  and  Z  that  currently  exist  in  internal  storage.  This  completes  the 
required  processing  for  a  given  data  pair,  (X,y),  and  then  the  instruction 
GOTO(START)  directs  the  computer  to  READ  the  "next"  data  pair  and  process  it 
in  the  same  manner  as  the  previous  pair.  After  all  of  the  cards  in  the  CARD 
READER  have  been  read  and  processed,  the  computer  stops  and  one  obtains  the 
results  in  the  form  of  the  punched  cards  from  the  bin  of  the  card  punch.  A  few 
of  these  output  cards  are  shown  in  FIGURE  8. 


A  tabulation  of  the  FORAST  program,  the  input  data,  and  the  results 
obtained  are  shown  in  FIGURE  9- 
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A  tabulation  of  the;  FORj^^ST  program,  the  input  data,  and  the  results  are  shown  in  2  3  JULY,  63  PAGE  I 
FIGUFiE  9.  PROS  M.  J.ROMANELLI  45107  EXAMPLE  1.  1 

START  READ(X)Y  2 

Z=SQRr(X**2+Y**2)  3 

PRINUX)  Y)Zi«G0T0(START)  4 

END  GOTO! START)  5 
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FIGURE  9. 


EXAMPLE  2. 


To  illustrate  a  convenient  means  of  identifying  the  "output"  produced  by 
the  computer j  we  consider  an  additional  recjuirement  in  the  prohlem  of  EXAT’IPLE  1. 
In  particular^  we  require  that  the  letter  X  be  printed  (approximately  centered) 
above  the  columns  corresponding  to  the  numerical  values  of  X,  and  similarly  the 
letters  Y  and  Z  centered  over  the  columns  that  correspond  to  their  respective 
numerical  values . 

To  provide  this  facility^  FOPvAST  accepts  a  PRIiJT  statement  wherein  we 
enclose  in  "special  quotes"  the  literal  characters  that  are  to  be  printed.  The 
"special  quotes"  are  the  less  than  and  greater  than  symbols,  <  >  .  The  PRINT 
statement  takes  the  form 

PRINT  <  > 


and  all  characters  within  the  <  and  >  symbols  are  recorded  on  the  output  card. 
For  our  purpose  we  write : 


PRINT  <  X  Y  Z  > 

Since  the  numerical  values  of  X,  Y,  and  Z  reqtiire  12  characters  in  the  standard 
form  of  output  that  we  obtained,  some  spacing  between  the  letters  enclosed 
within  the  <  >  is  necessary  if  we  desire  that  the  letters  be  approximately 
centered  over  their  respective  columns.  Hence  if  we  choose  to  have  the  letter 
X  over  column  6,  Y  over  column  l8,  and  Z  over  column  50,  we  indicate  this  spac¬ 
ing  by  writing 


PRINT  <  (5b)  X  (111^  Y-fliy.Z  > 


and  (11^  indicate  to  a  key  punch 


where  the  circled  quantities,  e)  '  _  ... 

operator  and  hence  to  the  computer  the  number  of  blanks  between  the  literal 
characters.  (This  is  the  only  place  in  the  FORAST  language  where  "blank" 
characters  are  not  ignored!  i.e.,  blank  characters  are  Ignored  everywhere 
except  between  the  <  and  >  symbols  of  PRINT  statements).  Knowing  how  to 
achieve  this  convenient  identification,  our  immediate  problem  now  is  to  deter- 
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mine  hov  to  "fit"  this  requirement  into  the  original  program  of  EXAi»lPLE  1. 
More  specifically,  if  ve  recorded  the  above  PRIM*  statement  on  a  card,  vhat 
is  its  logical  place  in  the  original  program?  If  ve  inserted  it  between 
cards  3  and  k,  i.e.,  in  front  of  the  card  that  instructs  the  computer  to 
PRINT  the  numerical  values  of  X,  Y,  and  Z,  ve  would  obtain  the  printing  of 
the  literal  letters  X,  Y,  and  Z  each  time  the  PRINT  <  X  Y  Z  >  was 
encountere  ci  •  FQRAST  in  structions  are  executed  in  sequence  one  after  another 
unless  a  specific  instruction  directs  the  computer  to  do  otherwise.  Since  we 
require  the  identification  to  be  printed  only  once,  a  logical  place  for  the 
identifying  PRINT  statement  is  at  the  very  beginning  of  the  program.  Con¬ 
sequently  the  following  program  will  produce  the  desired  results. 


ADDREQ 

START 


PROS  C906  M.J.  ROMANELLI  EXAMPLE  2. 
PRINT  <  (2)  X  Y  Z  > 

READ(X)Y 

Z  =  SQRT(X**2  +  Y**2) 
print(x)y)z  ^  GOTO(START) 

END  GOTO(ADDREQ) 


1 

1.1 

2 

3 

4 

5 


Note  that  card  5  differs  from  the  original  in  that  instead  of  directing  the 
computer  to  begin  by  reading  a  card  at  START,  it  directs  the  computer  to  begin 
at  ADDREQ  where  it  is  instructed  to  print  a  card  with  the  desired  X  Y  Z 
heading.  The  above  program  and  the  results  produced  are  listed  in  FIGURE  10. 


29 


23  JULY, 63  PAGE  1 

P’ROB  C906  M.J.  ROMANELLI  EXAMPLE  2.  1 

ADDREQ  PRINK  X  Y  Z>  1.1 

START  READIIX)Y  2 

Z=SQRT(X**2+Y**2)  3 

PRINKX)  Y)ZX  GOTO(START)  A 

END  GOTOII  ADDREQ)  5 

10000000  01  10000000  01 
30000000  02  40000000  02 
-80000000  01  60000000  01 
76000000  00-50000000-01 
-40000000  01-50000000-01 
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PROB 

C906  M.J. 

ROMANELLI 

EXAMPLE  2. 

» 

X 

Y 

Z 

0000001 

10 000000 

1  10000000 

1 

14142136 

1 

0000002 

30000000 

2  40000000 

2 

50000000 

2 

0000003 

-80000000 

1  60000000 

1 

10000000 

2 

0000004 

75000000 

-50000000- 

01 

75166482 

0000005 

-40000000 

1-50000000- 

01 

40003125 

1 

0000006 

FIGURE  10. 


NAMES  AND  OPERATIONS 


EXAMPLE  1  &  2  illustrated  the  substitution,  evaluation,  internal -external 
recording  and  controls  required  to  obtain  a  high-speed  digital  computer  solution 
of  a  simple  problem.  Common  to  each  of  these  processes  is  the  fundamental  re¬ 
quirement  for  identifying  particular  quantities  and  operations  of  interest. 
PORAST  uses  names  constructed  from  a  restricted  set  of  symbolic  characters. 

In  EXAMPLE  1,  we  used  the  symbolic  names  X  and  Y  to  identify  variables,  Z  to 
identify  a  particular  function,  SORT  another  function,  and  **  denoted  parti¬ 
cular  arithmetic  operations,  READ  and  PRINT  as  names  of  recording  operations, 
and  START  to  describe  an  instruction. 


To  avoid  ambiguities  in  the  definition  of  names,  operations,  and  numerical 
values,  the  49  characters  of  the  PORAST  language  are  divided  into  two  classes: 

Class  I  .'ABC*''Z  012‘"’9  58  characters 

Class  II  +-*/()=^^  11  characters 

The  characters  of  Class  II  have  special  meanings  and  must  not  be  used  to  define 
non-indexed  names  of  variables,  functions,  statements  or  other  parameters. 

+  is  used  to  denote  addition 

is  used  to  denote  subtraction 

*  is  used  to  denote  mjultipli cation 

/  is  used  to  denote  division 

’jo  is  used  to  denote  end  of  statement 

=  is  used  to  denote  (equality)  evaluation 

,  is  used  to  denote  an  indexed  (subscripted)  name 

**  is  used  to  denote  exponentiation 

<  is  used  to  denote  "less  than" 

>  is  used  to  denote  "greater  than" 

(  )  are  used  to  denote  miiltipli cation,  enclose  argimients,  indicate 

order  of  operations,  separate  names,  parameters,  numerical  values, 

C  UVJ  • 
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is  used  to  denote  that  the  information  on  a  card  following  the 
5^5^  is  to  be  ignored  by  the  computer. 

Non-indexed  symbolic  names  of  variables,  functions,  statements  and  para¬ 
meters  may  be  constructed  by  combining  characters  of  Class  I  "onder  the  follow¬ 
ing  restrictions: 

1 .  The  leading  character  must  not  be  zero  and  at  least  one  character 
must  be  other  than  a  decimal  digit. 

2.  Non-indexed  names  must  not  exceed  6  characters  in  length  unless 
those  after  the  first  6  are  not  required  for  unique  identification. 

5.  The  special  names  such  as  SELF,  SIN,  COS,  READ,  PRINT,  etc.,  have 

been  reserved  and  shoiild  not  be  used  as  arbitrary  names.  (A  complete 
list  of  reserved  names  is  given  on  page  59  of  [_!•]  • 

4.  Names  of  indices  should  not  exceed  three  characters  in  length,  4  are 
permitted  only  if  the  rightmost  2  are  decimal  digits. 

Despite  the  above  restrictions,  considerable  freedom  in  the  choice  of  names  is 
permitted  as  illustrated  in  the  examples  given  below: 


X 

X« 

BOX  1. 

X3VEL. 

LAST 

Y8 

XDOT 

IBOX 

Z” 

SINA 

17A 

XBAR 

YACC 

DZ" 

TANX 

IX 

START 

VELZ 

8.9 

BETA2 

XI 

SAM 

FOURTH 

10.2.4 

9GAMMA 

X15A 

WORK 

Y» 

A.  'B' 

DELTAX 

F300 

N.Y. 

L.A. 

FLA. 

TEMPO 

EVALYS 

EPS 

RHO 

ATPRES 

DENSTY 

FOFX 

FOFY 

FOFZ 

FOFZ" 

FINIS 

Many  problems  require  the  processing  of  groups  or  arrays  of  numerical 
data,  particularly  discrete  functions,  vectors,  matrices,  etc.  Conventional 
mathematical  notation  provides  identification  by  subscripts  or  indices.  For 
example ,  ,  (^i  =l,2,...,n^  j  =0,c:,4,...,m^  j  ^3'’  "  "  ’ 


Z 


k* 
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the  FORAST  language j  we  use  the  special  character ^  comma,  to  denote  an 
indexed  name ;  i . e . , 

is  represented  as  X,I 

Y.  is  represented  as  Y.J 

J 

Z,  is  represented  as  Z,K 

K 

As  is  conventional  in  mathematical  notation,  the  indexed  names  take  on  specific 
meanings  which  depend  on  the  name  before  the  comma  and  the  existing  value  of 
index  whose  symbolic  name  appears  after  the  comma.  Flexible  and  convenient 
means  for  setting  and  manipulating  index  values  is  provided;  however,  it  is 
important  to  recognize  that  the  internal  form  of  representation  for  index 
values  differs  from  the  floating-point  representation  discussed  previously. 
Index  values  are  restricted  to  a  limited  range  of  integers  and  consequently 
combining  index  values  and  floating-point  values  is  not  permissible.  The 
language  does  provide  for  transforming  from  one  form  to  another.  The  con¬ 
venient  use  of  index  names  and  means  for  establishing  and  generating  desired 
values  for  them  will  be  illustrated  in  later  examples.  The  index  names  were 
introduced  here  briefly  to  indicate  the  general  means  of  identification  by 
symbolic  name  permitted  in  the  FORAST  language. 

We  introduce  next  the  means  of  denoting  the  fundamental  arithmetic 
operations  of  addition,  subtraction,  multiplication  and  division,  which  are 
denoted  by  the  special  characters,  +  -  *  and  /  respectively.  Hence,  the  sum, 
difference,  product  and  quotient  of  two  quantities,  A  and  B,  are  expressed  as 
follows : 

A  +  B 
A  -  B 
A  *  B 
A  /  B 
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The  special  symbol^  is  used  to  denote  multiplication  to  distinguish  the 
product  A  times  B  from  the  single  quantity,  AB.  Other  means  of  denoting  a 
product  are 


(a)(b) 

or 

A)(B) 

or 

a)b 

or 

A(B 

or  even  the  redundant  symbol,  *,  is  permitted  in  the  above  forms.  In  general, 
parentheses  may  be  used  to  denote  multiplication  or  to  group  operations  in  any 
desired  sequence.  In  the  absence  of  parentheses,  particularly  in  ambiguous 
expressions,  priority  rules  govern  the  sequence  of  operations. 

When  the  order  of  operations  is  not  specified,  the  established  priorities 
are  as  follows: 

1. )  Single-valued  functions  of  one  argument; 

2. )  Exponentiation 

3. )  Multiplication  and  Division 

U.)  Addition  and  Subtraction 

Next,  we  introduce  the  means  to  denote  evaluation  of  the  convenient 
elementary  functions  such  as  sine,  cosine,  square  root,  etc.  These  functions 
are  denoted  by  special  names  as  is  conventional  in  mathematics.  For  example: 

sine  is  denoted  by  SIN 
cosine  is  denoted  by  COS 
is  denoted  by  SORT  . 

We  classify  these  and  similar  functions  as  single -valued  functions  of  a  single - 
argument.  (This  is  not  to  be  interpreted  as  a  single-valued  function  of  one 
variable,  on  the  contrary,  the  argument  may  be  a  function  of  many  variables  and 
indeed  many  other  functions.  We  emphasize  the  classification  of  this  set  of 
convenient  functions  to  distinguish  them  from  another  set  which  includes 
functions  of  more  than  one  argument). 
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X  =  y  +  Z  *  W  implies  x  =  y  +  zw 

since  miiltiplication  has  priority  over  addition. 

X  =  Y  -  Z/W  implies  x  =  y  -  — 


:e  division  has  priority 


X  =  Y  +  Z  **  2  implies  x  =  y  + 

since  exponentiation  has  priority  over  addition; 

X  =  SIN(Y+Z)**3  implie  s  X  =  sin'^(y+z) 

since  single-valued  functions  have  priority  over 
exponentiation . 
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EXAMt'LE  3. 


This  example  is  given  to  illustrate  the  single -valued  functions  of  single - 
arguments  that  are  available  in  the  FOPtAST  programming  language . 

GIVEN; 


Numerical  values  of  two  variables,  x  and  y.  Assimie  that  they 
3iVQ  A*0cord.6d.  in  si.s.nd.s.nd.  f*lo9,*tin§*poin"b  fomn  on  2.  pvincli0d.  c227ci« 
For  simplicity,  we  will  let  x  =  1  and  y  =  2. 

Compute  the  following  functions;  identify  and  record  each 
function  on  punched  cards. 


A  =  yy 
B  =  sin  X 
C  =  sin  xy 
D  =  sin  4x 


E  =  sin  5x 
F  =  sin  (-x) 
G  =  cos  x 


H  =  cos  (x  +  rt) 
I  =fn(y  + 


L  =  arctan  (x  -  y) 
M  =  arctan  (y) 

N  =  arccot  (-x) 

0  =  arccot  (y) 

P  =  arccos  (x/y) 


Q  =  arccos  (-x/y) 


R  =  arc sin  x/5 

S  =  arcsin  (-x) 

T  =  tan  xy 

U  =  tan  (-x) 

V  =  cot  (y  -  x) 

W  =  cot  (-y) 

FX  =  sec  X 

Fi  =  CSC  y 

FI  =  sign  of  X 

F2  =  sisn  of 
_  —  - - \ 

F3  =  sign  of  (y  -  y) 

F4  =  sinh  x 

F5  =  cosh  X 

F6  =  tanh  X 

Ff  =  integer  part  of  A 

f8  =  fractional  part  of  A 

2  2 

F9  =  sin  y+  cos  y 
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The  program,  input  and  output  for  this  example  are  listed  in  i?xGURE  12. 


BEGIN 


FLOW  CHART  FOR  EXAMPLE  5. 
FIGURE  11. 
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I 

1 


COMM  SINGLE-VALUED  FUNCTIONS  OF  ONE  ARGUMENT  2 

BEGIN  KEAO(X)Y  3 

PRINT<  X  =  >X<  Y  =  >Y  4 

A=SQRT(Y)«  PRINT<  SORT(Y)  =  >A  5 

B  =  SIN(X)-«  PHINT<  SIN(X)  =  >8  6 

C=SIN(X*Y)«  PRINT<  SIN(XY)  =  >C  7 

D  =  SIN(4*X)?;  PRINT<  SIN(4X)  =  >U  8 

E=SlN(X*b)I  PRINT<  SIN(5X)  =  >E  9 

F  =  SIN(-X)*  PRINT<  SliM(-X)  =  >F  10 

G  =  COS(X)-<;  PRINT<  COS(X)  =  >G  11 

H  =  COS  (  X4-3.  14159)  ^  PKINT<  COSIX  +  PI)  =  >H  12 

I =LOG ( Y  +  SQRT ( X) ?  PR  I  NT <LN ( Y  +  SQRT ( X )  =  >I  13 

J=LOG10(Y)%  PRINT<  LOGIO(Y)  =  >J  14 

K  =  EXP(X<-Y)S  PRINT<  EXP(X+Y)  =  >K  15 

L=ARCTAN( X-Y) «  PRINT<  ARCTAN(X-Y)  =  >L  16 

M  =  ARCTANIY)?  PRINT<  ARCTAN(Y)  =  >M,  17 

N  =  Ai<CCOT  (-X)  *  PRINT<  ARCCOT(-X)  =  >N  18 

0=AKCC0r(Y)^  PRINT<  ARCCOT(Y)  =  >0  19 

P  =  AKCCOS(  X/Y)  Pf(INT<  ARCCOS(X/Y)  =  >P  20 

g=ARCCOS( -X/Y)?  PRINT<ARCCOS(-X/Y )  =  >Q  21 

R  =  AKCSIN(  X/5)  PRINT<  ARCSIN{X/5)  =  >R  22 

S=ARCSIN(-X) «  PR1NT<  ARCSIN(-X)  =  >S  23 

r  =  TAN(X*Y)?;  PRI(nT<  rAN(XY)  =  >T  24 

U=TA.m(-X)%  PRINT<  TAN(-X)  =  >U  25 

V=COr(Y-X)«  PRInT<  COT(Y-X)  ^  >V  26 

w  =  CUT(-Y):«  PRINT<  COTl-Y)  =  >W  27 

FX=SEC(X)%  PRINT<  SEC(X)  =  >FX  28 

FY=CSC(Y)?  PRINr<  CSC(Y)  =  >FY  29 

Fl  =  SIGN(X)ii  PRI;\iir<  SIGN(X)  =  >F1  30 

F2=SIGN(-X)4  PRINT<  SIGNl-X)  =  >F2  31 

F3  =  SIGNl  Y-Y)  as  PRINT<  SIGN(O)  =  >F3  32 

F4=S1NH(X)X  PRINT<  SINH(X)  =  >F4  33 

F5  =  C0SH(X)S;  PRINT<  COSH(X)  =  >F5  34 

F6=rANH(X)S;  PRINT<  TANHIX)  =  >F6  35 

F7=WH0LE  (  A )  a;  PRINT<  WHOLE(A)  =  >F7  36 

F8=FRACT(A)^  PKINT<  FrACT(A)  =  >F8  37 

F9= ( SIN( Y ) *»2+C0S( Y)**2) **.5«  PRINT<  CHECK  =  >F9  38 

GOTO(N.PKOB)  39 

END  G0TO(BEGlN)i  40 

lOOOUOOO  01  20000000  01 
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X  =  10000000  1  Y  = 

20000000  1 

0000001 

SQRT( Y) 

zi 

14142136 

1 

0000002 

S  I  N  (  X  ) 

= 

84147098 

0000003 

SIN( XY) 

s 

90929743 

0000004 

SIN(4X) 

= 

-75680250 

0000005 

SINJ6X) 

s 

-95892427 

0000006 

SIN(-X) 

-84147098 

0000007 

cost  X) 

= 

54030231 

0000008 

COS ( X+PI ) 

-54030454 

0000009 

LN(Y+SCRT(X) 

10986123 

1 

0000010 

LOGIO(Y) 

30103000 

0000011 

EXP  (  X4-Y) 

= 

20085537 

2 

0000012 

ARC  TAiNt  X-Y) 

-78539816 

0000013 

ARCTANt  Y) 

1 1071487 

1 

0000014 

ARCCOT (-X) 

-78539816 

0000015 

ARCCOT(Y) 

=. 

46364761 

0000016 

ARCCOSt  X/Y) 

s 

10471976 

1 

0000017 

AKCCOSt-X/Y) 

5 

20943951 

1 

0000018 

ARCSINtX/S) 

= 

20135792 

0000019 

ARCSliMt-X) 

S 

-15707963 

1 

0000020 

TANt  XY ) 

S 

-21850399 

1 

0000021 

TAW(-X) 

-15574077 

1 

0000022 

CO  rt Y-X) 

S 

64209262 

0000023 

COT(-Y) 

= 

45765755 

0000024 

SCC(X) 

•= 

18508157 

1 

0000025 

CSCt  Y) 

10997502 

1 

0000026 

SlGNt  X) 

= 

10000000 

1 

0000027 

SIGWt-X) 

■:= 

-10000000 

1 

0000028 

SIGNtOI 

OOOOOOOO 

0000029 

SiNHt  X) 

= 

11752012 

1 

0000030 

CUSHt  X) 

15430806 

1 

0000031 

FAHHIX) 

76159416 

0000032 

WHOLE (A) 

10000000 

1 

0000033 

FKACT ( A) 

= 

41421356 

0000034 

CHECK 

10000000 

1 

0000035 

FIGURE  12. 


Restrictions  on  the  single -valued  functions  of  one  argument  illustrated 
in  EXAMPLE  5  are: 

1. )  Every  argument  is  enclosed  in  parentheses,  the  closing  right 

parenthesis  is  optional. 

2.  )  Every  argument  and  every  result  is  in  floating-point  form. 

5.)  The  arguments  of  trigonometric  functions  and  results  of  inverse 
trigonometric  functions  are  in  radians . 

4. )  The  results  of  ARCTAN  and  ARCSIN  are  in  the  interval 

-  ^  <  Result  <  ^ 

5. )  The  results  of  ARCCOT  and  ARCCOS  are  in  the  interval 

0  <  Resillt  <  n 

6. )  The  arguments  may  be  functions  of  many  variables  and  fiuictions; 

however,  only  the  sign  of  the  resulting  argument  is  considered  in 
the  evaluation  of  the  particular  function. 

Although  not  illustrated,  names  of  variables,  arguments  and  results  may 
be  indexed.  Violations  of  restricted  bounds  of  arguments  will  be  detected  and 
recorded  during  the  execution  of  the  program. 
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EXAMPLE  4. 


EXAMPLE  5.  illustrated  the  available  functions  of  one  argument.  Another 
set  of  convenient  functions  are  those  which  are  functions  of  more  than  one 
argument  and  often  produce  more  than  one  result.  Included  in  this  set  are 
methods  for  interpolation,  qiiadrature ,  matrix  operations,  solution  of  ordinary- 
differential  evj_uations^  etc.  To  instruct  the  computers  to  evaluate  functions 
in  this  set  -we  -write  a  statement  of  the  form. 


L.'  r\ir  I k  L^j  t  TVT ^ 


where  the  word  ENTER  followed  by  a  particular  name  enclosed  in  parentheses 
identifies  the  particular  function  or  method  of  interest.  The  N1,W2, ..... 
represent  the  names  or  values  of  the  arguments,  parameters  and  the  defining 
names  of  the  results .  Many  of  the  available  functions  in  this  set  will  be 
illustrated  in  separate  examples. 

To  introduce  one  of  these  functions,  we  will  consider  an  extension  of 
EXAMPLE  1.  Assume  that  the  given  pairs  (x,y)  represent  the  rectangular 
coordinates  of  a  point  in  a  plane  and  that  we  want  to  produce  a  table  which 
includes  the  corresponding  polar  coordinates  of  the  point,  i.e.. 


GIVEN: 


REQUIRED: 


}  i  =  . . . 


Corresponding  to  each  point,  (x,y)  detexiiiine, 
list  in  a  table  the  polar  coordinates  and 
degrees ) . 


■;  ^  4- j 

±u.tz;m^±iy  aiiu 

0^.  (list  &  in 
1 


o . 
"1 


S'.  =  arctan  (y./x.) 

1  ■'i'  1 


kl 


Before  constructing  the  flow-chart  and  writing  the  program  which  will 

produce  the  desired  solution  we  emphasize  the  distinction  between  the  available 

arc tan  function  of  one -argument  illustrated  in  EXAMPLE  and  the  arc tan 

function  of  two -arguments  to  be  used  and  illustrated  in  this  example.  The 

distinction  is  that  the  arctan  function  of  one -argument  produces  a  result  which 

lies  in  quadrants  I  or  IV.  (Note  4.  on  page  38^-21  <  Result  <  £  )  •  Since 

2  2 

the  arctan  function  of  two  arguments  take  into  account  the  signs  of  both 
arguments,  it  produces  a  res\ilt  which  may  be  in  quadrants  II,  III  or  IV,  i.e., 

-n  <  Result  <  n  . 


START  READ  PT 


FLOW  CHART  FOR  EXAMPLE  4. 
FIGURE  13. 
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I’he  program,  input  and  output  for  this  example  are  listed  in  FIGURE  l4. 

PK08  C906  M.J.  KOMANELLI  EXAMPLE  4. 

COMM  TKANSFOKMATION  OF  RECTANGULAR  TO  POLAR  COORDINATES 
START  PRINT<  X  Y  RHO  THETA> 

REAOPT  READlXlYi  RHO-SCRT ( X»X+Y*Y ) 

ENT£R(ARTAN  )Y)X)TH 
irH==57.29S70O*rH 


■p- 


END 

10000000 

30000000 

-80000000 

75000000 

-40000000 


PRINKX)  Y)Rh0)  TH'-S  GOTO(READPT) 
GOTO (START) 

Ul  10000000  01 
02  40000000  02 
01  60000000  01 
00-50000000-01 
01-50000000-01 
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X 

10000000 

30000000 

-80000000 

75000000 

-40000000 


Y 

1  10000000 
2  40000000 
I  60000000 
-50000000- 
1-50000000- 


RHO 

14142136  I 
50000000  2 

10000000  2 
75166482 
40003125  I 


THETA 

45000000  2 

53130103  2 

14313010  3 

-38140749  1 

-17928384  3 


1 

2 

1 

01 

01 


1 

1 

2 

3 

4 


5 


6 

7 

8 


0000001 

0000002 

0000003 

0000004 

0000005 

0000006 


FIGURE  14 


EXAMi^LE  5 . 


In  the  previous  examples ^  the  output  results  were  recorded  in  standard 
floating-point  form.  In  this  example  we  will  illustrate  a  means  of  obtaining 
the  tabular  results  in  a  form  with  the  decimal  point  actually  punched  (and 
printed)  instead  of  the  exponent.  For  example,  instead  of  obtaining  the 
printed  value  for  J2  as  lhlh21^6  1,  we  will  obtain  the  printed  value  for 
as  1.4l4  .  This  example  will  illustrate  a  special  case  of  the  general 
PRIUT  statement  where  a  departirre  from  standard  form  is  desired  and  must  be 
specified  accordingly. 


To  designate  a  non-standard  form,  the  PRINT  statement  is  written  in  the 


form 


PRINT -FORi'iAT  (Name  of  specifying  form)  -  (  ))  ••=  )  ^ 


where  the  "Name  of  specifying  form",  enclosed  in  parenthesis  after  the  word 
FORMAT,  identifies  a  quantity  which  explicitly  defines  forms  and  horizontal 
spacings  desired.  We  refer  to  this  specifying  quantity  as  a  "format  word". 
To  provide  for  arbitrary  spacing  of  the  output  forms  and  for  transformations 
from  internal  to  external  forms  the  format  word  is  composed  of  parenthesized 
expressions  of  the  form 

(T  -  S  -  L) 


where  T,  S,  and  L  are  decimal  digits  with  defined  meanings: 


T 

S 

L 


generally  specifies  a  type  of  transformation,  a  repetition  of  a 
previous  format,  spacing  of  quantities,  end  of  format,  etc.  ; 


generally  specifies  a  scale  or  relative  location  of  a  decimal  point; 


generally  specifies  the  length  of  a  field,  i.e., 
columns  or  characters  to  be  used  in  representing 


the  number  of 
a  given  form. 


.) 


The  parenthesized  expressions  may  contain  all  three  of  the  above  descriptors^ 


(T  -  S  -  L) 

or  only  two  descriptors,  (T  -  L) 
or  only  one  descriptor  ,  (T) 

Similar  expressions  may  be  constructed  for  use  with  the  general  READ  statement. 
The  permissible  values  of  T  and  their  corresponding  meanings  are  given  in  ^l.J  . 
Only  a  few  of  the  most  frequently  used  T*s  will  be  illustrated  here. 

In  EXAMPLE  4.  the  entries  in  the  output  table  were  expressed  in  standard 
floating-point  form;  i.e.,  coefficients  with  corresponding  exponents.  To 
depart  from  the  standard  form  and  obtain  numerical  values  with  decimal  points 
printed  instead  of  exponents,  we  will  specify  (using  a  format  word)  the  parti¬ 
cular  transformations  desired,  the  relative  location  of  the  decimal -point,  the 
number  of  columns  to  be  used  for  each  quantity  and  the  spacing  between  indivi¬ 
dual  quantities.  In  the  standard  form  we  used  12  columns  for  each  quantity  and 
provided  for  no  spacing  between  them. 

Assume  that  we  want  the  following  format  for  the  table  corresponding  to 
the  output  of  EXAMPLE  4. 


+  DDD.DDD 


DDD.DDD 


DDDD.DD 


THETA 
+  DDD.DDD 


where  the  D’s  represent  decimal  digits.  That  is,  we  will  specify  a  total  of 
8  characters  for  each  quantity  and  4  blardc  spaces  between  each  quantity.  Note 
that  in  the  total  length,  8,  we  include  the  algebraic  sign  and  decimal  point 
in  the  count.  (We  chose  four  spaces  between  quantities  since  the  column  head¬ 
ings,  X,  y,  PRO  and  THETA  were  approximately  centered  over  12  column  entries; 


henee,  with  the  choice  of  8  characters  for  each  quantity  and  4  spaces  between 
quantities,  we  will  not  have  to  alter  the  spacing  of  the  characters  of  the 
column  headings). 

To  specify  this  format  we  construct  a  format  word  which  we  arbitrarily 
label  SEEl. 


LOG  O.T. 

1  6  T  : 

SEEl .  FORM 


rORMULAS,  STATEI''uiji'jTS ,  COiwIENTS  IDEN. 

11  TT-80 

(5-2)12-3-8)5-4)12-5-3)5-4)12-4-8)5-4)12-5-8)2  7-1 


Note  that  the  arbitrary  name  chosen  for  the  format  word  is  written  in  the 
LOCATION  columns  and  the  word  FORM  in  col'amns  7-10  of  the  coding  sheet, 
parenthesized  expressions  are  written  in  columns  11-70. 


The  PRINT  Statement  that  refers  to  this  format  w'ord  is  written  in  the 
following  form: 


PRINT -FCeMAT( SEEl.)  -  (X)y)RH0)T  ^ 


7’ 


The  expressions  (5-L)  in  the  format  word,  SEEl.,  correspond  to  T  =  5  which 
denotes  that  a  spacing  of  L  columns  is  desired.  (Note  that  we  begin  with  2 
spaces  before  the  first  quantity,  X,  and  denoted  spacings  of  4  thereafter). 
The  expression  (12-5-8)  corresponds  to  T  =  12 

S  =  5 
L  =  8 

and  represents  the  form  for  the  quantity,  X. 

T  =  12  denotes  that  the  quantity  exists  internally  as  a  floating-point 
number  and  the  output  form  desired  is  a  "fixed"  form  with  the 
"fixed"  decimal-point  location  specified  by  the  corresponding  3. 

S  =  3  denotes  that  the  "fixed"  decimal -point  is  to  be  recorded  after 
the  5  digits  that  follow  the  algebraic  sign. 

L  =  8  denotes  that  a  total  of  8  columns  (or  characters)  are  to  be  used 
for  the  quantity. 
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Note  the  identical  forms  for  the  quantities  X,  'i,  and  THETA  and  the  decimal- 
point  after  h  digits  in  the  representation  for  RHO.  The  2  after  the  last 
parenthesized  expression  corresponds  to  T  =  2  and  denotes  the  end  of  the 
format  word. 

To  obtain  the  desired  results,  we  will  replace  Card  7  of  EXAMPLE  U.  with 
the  card  identified  as  7’  above,  (i.e.,  the  PRINT  Statement  that  refers  to  the 
specified  format),  and  add  the , card  identified  as  7*1  that  defines  the  speci¬ 
fied  format .  The  program,  input  and  output  are  listed  in  FIGURE  15 . 


23  JULY, 63  PAGE  1 

PROB  C906  M. J.ROMANELLI  46107  EXAMPLE  6  1 

COMM  TRAtMSIFORMAriON  OF  RECTANGULAR  TO  POLAR  COORDINATES  2 

START  PRINT<  X  Y  RHO  THETA>  3 

REAUPT  READ(X)YX  RUO= S QRT U ♦X +Y *Y )  4 

ENTER  TARTAN  ) Y ) X ) TH  5 

TH=57.295780*TH  6 

PRINT-FOKMAI (SEEl.)-(X)Y)RHO)TH)?  GOTO ( READPT ) *  7* 

SEEl .  FORfU  3-2 ) 12-3-8 ) 3-4 ) 1 2-3-8 ) 3-4 ) 1 2-4-8 ) 3-4 ) 1 2- 3-8 ) 2  7.1 

END  GOTO (START)  8 


lOOOOOOO  (J1  10000000  01 
30000000  02  40000000  02 
-80000000  01  60000000  01 
75000000  00-50000000-01 
-40000000  01-50000000-01 


.23,63 

PKOB 
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EXAMPLE  5 

« 

X 

Y 

RHO 

THETA 

0000001 

1.000 

l.OOO 

1.41 

45.000 

30.000 

40.000 

50.00 

53. 130 

8.000 

6.000 

10.00 

143.130 

,  750 

.050 

.  75 

-  3.814 

4.000 

.050 

4.00 

-179.284 

FIGURE  15. 


In  addition  to  accepting  FORMITLA  statements  which  involve  addition, 
subtraction,  multiplication,  division,  exponentiation  and  many  s ingle -value d 
functions  of  one-argument,  FORAST  also  accepts  eleven  English  word  statements 
and  twenty "two  special  English  words.  The  latter  are  called  '*pseudo  order- 
types”  and  are  labelled  "ORDER  TYPE"  over  columns  7  through  10  of  the  standard 
coding  form.  Only  11  of  the  22  permissible  order  types  will  be  illustrated  in 
this  report.  Included  in  the  eleven  English  word  statements  is  the  general 
EiWER  statement  wiiich  provides  for  evaluations  of  approximately  forty  functions 
of  more  than  one  argument. 

We  have  in  the  previous  examples  illustrated  special  cases  from  each  of 
the  above  categories.  Each  example  included  at  least  one  English  word  state¬ 
ment  and  at  least  one  single -valued  function  of  one  argument.  The  READ,  PRINT, 
GOTO  and  ENTER  statements  illustrated  4  of  the  11  English  word  statements. 

PROB  and  END,  (which  are  required  in  every  problem),  illustrated  two  of  the 
special  O.T.  English  words.  Listed  below  are  the  11  English  word  statements 
and  the  11  order-types  to  be  illustrated  in  this  report. 


ENGLISH  WORD  STATEMENTS 


SPECIAL  ENGLISH  WORDS  (O.T.) 


GOTO 

SET 

SETEIA 

INC 

rtr\T  Tnim 

IF 

CLEAR 

MOVE 

ENTER 

READ,  PRINT  (PUNCH) 
HALT 


PROB 

(Problem 

El® 

DATE 

COMM 

( Comment ) 

nrMVTm 

( Continue ) 

FORM 

(Format) 

LIST 

(Listing) 

BLOC 

SYN 

(Synonym) 

DEC 

(Decimal) 

DEC  = 

(Decimal  Equality) 

X  VI  — 


1  TTnrrl-Jo'h  a+Q  +  ^iTnon+o  Q-n-noar^  omrurViovo  n  q  ■nT’/^rrr'QTn 

V.'XXV-.  xjXX^_l.XOii  W^X  k_r  s/\.A  XX  <-*xxjr  n  xx\.«x.  u.xx  <-«• 


and  are  transformed  to  computer  instructions  that  will  be  executed  during  the 


k9 


XT _ 


running  oi  one  program.  The  special  English  words  are  not  transformed  into 
computer  instructions^  the  words  themselves  may  only  appear  in  columns  J 
through  10.  (Note  that  none  of  these  names  exceeds  four  characters).  Rather 


'♦"ViQ'n  o  e  "Ko  oq/-»Vi  ■? 


••  V  Vw  X 


n  their  full  generality,  we  will 


as  in  the  previous  examples  to  introduce  a  few  in  each  example,  with  vari¬ 
ations  illustrated  in  succeeding  examples. 


Many  problems  req_uire  definitions  of  functions,  or  more  generally,  many 
problems  require  control  of  processes  ^rtiich  are  based  on  conditional  relation¬ 
ships.  For  example,  it  may  be  necessary  to  evaluate 


f(X)  =  A  +  EX  +  CX^  if  X  <  1.5 
f(X)  =  D  +  EX  if  X  >  1.5 


We  assume  that  X  may  take  on  values  in  either  range  but  as  stated  above 
f(X)  is  defined  accordingly.  Schematically,  the  above  may  be  illustrated  as 
follows: 


We  have  arbitrarily  labelled  one  of  the  definitions  DEF2  and  indicated  that 
after  the  appropriate  definition  is  applied  (i.e.,  after  the  appropriate 
function  is  evaluated)  the  computer  is  to  continue  operations  at  a  place  called 
WORK.  To  instruct  the  computer  in  the  FORAST  language  to  carry  out  the  above 
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plan,  "we  write: 


LOG.  O.T.  FORMULAS,  STATEMEOTS,  COMMEITOS 

IF(X  >  1.5)G0T0(DEF2)^  FX  =  a  +  X(B  +  C  *  X)^  GOTO(WORK) 
DEF2  IX  =  D  +  E  *  X 

WORK 

The  above  is  an  illustration  of  a  simple  conditional  statement.  It  has  the 
form 

IF  (a  specified  condition  is  satisfied)  GOTO( someplace 

As  indicated  in  the  form,  if  the  specified  condition  is  satisfied  the  computer 
is  directed  to  GOTO  the  someplace  denoted  for  subsequent  instructions,  other¬ 
wise  the  computer  will  execute  the  statement  immediately  following  the  con¬ 
ditional  statement.  It  is  to  be  emphasized  that  the  above  is  an  illustration 
of  a  simple  conditional  statement,  more  general  compound  forms  are  permissible. 
The  conditional  expressions  may  contain  <  ,  >  ,  and  equality  relations  whose 
terms  involve  the  arithmetic  operations  and  single -valued  functions  of  one 
argument.  Indeed  a  single  conditional  statement  may  contain  many  conditional 
expressions  separated  by  the  logical  operations  AND  or  OR  and  each  expression 
may  be  prefixed  with  ABS,  NOT,  INT,  etc.,  where:  ABS  denotes  absolute  value; 

NOT  denotes  negation  of  the  parenthesized  relation,  INT  denotes  that  the 
quantities  Involved  in  the  relation  are  integers  in  integer  form,  (not  floating¬ 
point  form).  In  the  following  example,  we  will  illustrate  the  simple  form  of 
a  conditional  statement. 
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EXAMPLE  6. 


GIVEN:  k,  B  and  N  recorded  in  standard  floating-point  form  in  the 

first  three  fields  of  a  piinched  card;  (Assume  several  such 
cards . ) 

f(X)  =  sin  X  ;  H  =  B  -  A 

N 


REQUIRED; 


Use  the  trapezoidal  rule  to  obtain  an  approximation  for  the 
definite  integral 


I 


f(X)dx 

4 


Trapezoidal  rule: 


I 


I  =  H 


f(A)  +  f (A+H)  +  f(A+2H)  + 
2 


+  f(B-H)  4  h|l 


Print  and  identify;  A;  B;  N  and  I  . 


Note  that  the  trapezoidal  rule  requires  that  the  function 
(integrand)  be  evaluated  at  (N+l)  discrete  values  of  X, 
the  end  values  (A  and  b)  are  weighted  l/2,  the  interior 
values  have  unit  weight. 

To  construct  a  flow  chart  that  outlines  a  plan  to  obtain  the  desired 
solution  we  begin  by  denoting  that  we  want  the  computer  to  read  the  card  that 
contains  the  pertinent  values  k,  B  and  N.  Recall  that  once  they  have  been 
read,  they  are  recorded  Internally  for  as  many  future  references  as  desired. 
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Hence  to  indicate  that  we  want  to  "begin  by  having  the  computer  read  a  card, 
we  write 


•RPrtTM 


N. 

READ; 

V 

> 

A,B,N 

Next,  we  observe  by  studying  the  definition  of  the  trapezoidal  rule  that  we 
need  the  value  H  for  two  purposes: 

1. )  it  is  needed  as  a  factor  to  obtain  the  final  result  after  the 

weighted  integrand  values  are  slimmed; 

2. )  it  is  needed  to  construct  the  discrete  interior  values  of  X,  i.e., 

A  +  H,  A  +  2H,  A  +  _3H,  etc. 

Hence,  since  we  need  H  for  the  two  purposes  given  above  we  append  to  the  chart 
a  box  indicating  the  desired  definition  and  evaluation  of  H. 


BEGIN 


READ: 

H  =  B  -  A 

7 

A;B,N 

7 

N 

Next,  noting  that  the  definition  requires  weights  of  l/2  on  the  "end"  terms, 
we  can  dispose  of  these  by  instructing  the  computer  to  evaluate  the  following 
function 

I  =  l/2(sin  A  +  sin  b) 

We  indicate  this  on  the  chart  accordingly 


csn  "bs 


Now  note  that  the  remaining  required  arguments,  A  4  K,  A  4  2H,  etc., 
obtained  in  general  by  adding  H  to  the  previous  argument.  If  ve  let  X 
be  an  existing  argument,  the  next  argument  can  be  obtained  by  adding  H  to  X. 
Symbolically  we  can  write  this  as 


V*  V*  •  TT 

=  .A.  n 

where  it  is  understood  that  the  X  on  the  right  of  the  equality  represents 
an  existing  value  which  will  be  used  to  generate  a  "new"  existing  value. 

Hence  we  have  generalized,  in  that  anytime  we  want  the  computer  to  generate 
a  new  argument,  we  need  only  direct  it  to  the  generalized  expression.  Note 
also  that  each  of  these  arguments  will  undergo  "similar  treatment",  i.e., 
each  is  the  argument  for  the  function,  sine.  Hence  we  can  generalize  a  two 
step  process  and  write  it  in  the  form 

X=X4H^  I=l4  SIN(X) 

If  we  let  X  take  on  the  initial  value  of  A  and  direct  the  computer  to 
the  above  two  step  process,  it  (the  computer)  would  in  the  first  step  evaluate 
the  function  for  existing  values  of  X  and  H  which  result  in  a  value  of 
X  =  A  +  H.  In  the  next  step,  the  computer  would  add  the  sin(A  4  h)  to  the 
existing  value  of  I.  Hence,  after  the  first  execution  of  the  two  step  process, 

I  would  then  represent  the  sum  of  three  terms  of  the  bracket  [  ]  ,  namely  the 

first,  second  and  last  terms.  To  schematically  illustrate  the  above,  we  augment 
the  flow  chart  as  follows: 


BEGIN 


W(s  havG  Si  uitiSiily  XauGllGu.  ohe  geriGraXizGu.  definition  of  arguinGnoS  and  "fciiG 
summing  of  the  genoraX  torm,  NXTEF^M.  We  wiXX  v/ant  to  make  reference  to  this 
process  and  direct  the  machine  to  carry  it  out  as  many  times  as  required. 

At  this  point  we  can  appXy  a  simpXe  conditionaX  statement. 


If  one  refers  to  the  definition  of  the  Trapezoidal  Rule,  it  wiXX  be 
noted  that  the  "last  term"  we  want  the  computer  to  add  to  the  existing  I  is 
f(B  -  H),  i.e.,  sin  (B  -  H),  (since  we  have  already  accounted  for  the  term 

f ^B^  ■ 

'  ).  Hence  we  need  only  instruct  the  computer  to  determine,  after  adding 

a  term  to  the  existing  I,  whether  the  corresponding  argument  was  equal  to 
B  -  H.  If  indeed  X  =  B  -  H,  all  of  the  required  terms  have  been  evaluated 
and  summed;  otherwise,  we  can  instruct  the  computer  to  "go  back"  to  evaluate 
and  add  the  next  term  to  the  existing  I.  Schematically  this  condition  can  be 
illustrated  as  follows: 


HXTERM  ^ 

1 

X  =  X  +  H 

I  =  I  +  sin  X 

' 

f 

X  =  B  -  H  ? 

No  1 

TiF 


The  corresponding  conditional  statement  in  the  FORAST  language  would  be  written 
as 


IF(X  =  B  -  H) GOTO ( FINISH GOTO(NXTERM) 

however,  since  the  computers  use  a  finite  number  of  places  in  the  calculations 
the  conditionaX  statement  written  above  may  not  be  adequate  for  all  values  of 
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B  and.  N. 


To  illustrate  this,  suppose: 


A  =  0 

B  =  1 

N  =  3  , 


'Z.  • 


then,  H  =  l/3  and  the  equivalent  internal  representation  =  .333' 
similarly,  B  -  H  =  2/3  and  the  equivalent  internal  representation  =  .666""-7 
(The  internal  representations  result  from  conventional  rounding  and  truncation 
procedures).  The  X  argument,  "corresponding"  to  B  -  H,  generated  by  the  speci¬ 
fied  definition  would  be  A  +  H  +  H  =  .666* • *6  !  Note  the  "6"  in  the  least 
significant  digit.  Hence  this  value  is  not  equal  to  the  .666* ••7,  i.e., 

.666' • ‘6  /  .666' • 'J .  In  conditional  expressions,  equality  relations  are  not 
satisfied  unless  the  resulting  numerical  values  on  both  sides  of  the  equality 
have  identical  representations.  Hence,  we  should  relax  the  stringent  equality 
condition  and  ask  if  X  differs  from  B  -  H  by  a  "tolerable"  amount,  say  e,  i-e., 


|X  -  (B-H)  I  <  € 

If  we  choose  e  =  ,  we  will  provide  for  a  maximum  tolerable  error.  We 

express  this  conditional  statement  in  the  following  form: 


if-abs(x-(b-h)  <  h/2)goto( finish)^  GOTO(NXTERM)^ 

When  ABS  precedes  a  conditional  expression,  the  computers  first  evaluate  both 
sides  of  the  conditional  expression,  then  take  the  absolute  values  of  both 
sides  before  checking  to  see  if  the  relation  is  satisfied.  In  essence,  the 
above  statement  corresponds  mathematically  to 

if  jX-B-rilj  <  go  to  finish,  otherwise  go  to  the  computation  of 

the  next  term.  At  the  location  called  FINISH,  we  want  the  computer  to  multiply 
the  existing  value  of  I  by  H  to  complete  the  definition  of  the  trapezoidal  rule. 
The  complete  flow-chart  is  shown  in  FIGURE  l6. 
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•RFnXTJ 
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FLOW  CHART  FOR  EXAMPLE  6,  (TRAPEZOIDAL  RULE). 
FIGURE  l6. 


The  program,  input  and  output  for  EXAMPLE  6  are  listed  in  FIGURE  17 . 
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21  JULY ,63  PAGE  I 

PROS  C906  M.J.  ROMANELLI  A5107  EXAMPLE  6.  1 


BEGIN 

REAU ( A ) H 1 N2 

H-= 

(B-A)/N«  I 

=  .5(S1N(A)+SIN|B)  ).«  X=A 

2 

ivJXTERM 

x  = 

X>H^,  l  =  l+SliN{  X)^  lF-ALJS(X-B  +  H<H/2)G0ir0{FINlSH)«  GO  TO  1  NX  TERM  )i 

3 

FINISH 

1  = 

H»n:  PKlNir<A 

=  >A<  6 

=  >8<  N  =  >N<  I  =  >i«  GOTO  (BEGIN!) 

4 

END 

GOTO{OEGIN) 

5 

OOOOOOUO 

00 

16707930 

01 

25000000 

02 

UOOOOOOO 

00 

16707930 

01 

50000000 

02 

0(DOOOOOO 

00 

16707930 

01 

10000000 

03 

00000000 

00 

16707930 

01 

20000000 

03 

00000000 

00 

15707930 

01 

30000000 

03 

00000000 

00 

15707930 

01 

40000000 

03 

00000000 

00 

16707930 

01 

60000000 

03 

MAY. 23,63  BKLLSC  FORAST  F62 


PROS  C906 

M. 

J. 

ROMANELLI 

45107 

EXAMPLE 

6. 

* 

A 

= 

00000000 

8 

= 

15707930 

1 

N  = 

25000000 

2 

I 

= 

99966767 

0000001 

A 

= 

00000000 

B 

= 

15707930 

1 

N  = 

50000000 

2 

I 

= 

99991443 

0000002 

A 

00000000 

6 

= 

15707930 

1 

N  = 

10000000 

3 

I 

= 

99997611 

0000003 

A 

= 

00000000 

6 

15707930 

1 

N  = 

20000000 

3 

1 

= 

99999153 

0000004 

A 

00000000 

D 

= 

15707930 

1 

N  = 

30000000 

3 

I 

■= 

99999439 

0000005 

A 

00000000 

B 

= 

15707930 

1 

N  = 

40000000 

3 

I 

= 

99999539 

0000006 

A 

= 

00000000 

B 

= 

15707930 

1 

N  = 

50000000 

3 

I 

99999585 

0000007 

FIGURE  17. 


To  iCj.Gn.'tiTy  "tliG  Ou.'tpu.'t  Q^u.s.n'ui'biGS^  ws  i13.vg  lisod.  3.  vRi'iR'tion  of  "tii© 
general  PRINT  statement.  The  statement 


PRINT  <  Ah='b>  A  <  hhBb=h>  B  <  bbNb=h>  N  <  bblb=b>  I 


instructs  the  computer  to  PRINT(PUircH)  the  name  of  the  quantity,  an  equality 
symbol,  followed  by  the  corresponding  numerical  value  of  the  quantity.  (The 
lower  case  b’s  merelv  serve  to  indicate  blanks  and  are  used  to  obtain  con¬ 


venient  bpaciiig  between  symbols,  numerical  values,  etc.), 
specified  between  <  and  >  in  PRINT  statements  are  printed 
ment  provides  for  output  in  the  form 


Since  all  characters 
,  the  above  state - 


A  - -  B  = -  N  = 


Observe  in  FIGURE  I7  that  the  numerical  values  are  represented  in 
floating-point  form  since  we  did  not  specify  any  format  to  depart 
standard  representation. 


n  +  Q  P 

from  the 
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iijCAMPLE  7- 


To  illustrate  a  concept  that  is  applicable  in  many  problems  we  will  write 
another  program  to  obtain  a  solution  to  the  problem  of  EDCAMPLE  6.  The  concept 
to  be  illustrated  is  one  in  which  we  will  "vary”  the  "exit”  of  a  general 
function  evaluation  so  that  we  will  direct  the  computer  to  various  destinations 

^  r~\  T  T  r~\  T  ^  ^  ^  r-\  r-\  Trt  -i-  ^  T  ^  n  i-i  \  -t-  In  m  m  im  T  iim  ^  i  m  ^  •?  m  v-\  n  m  in  n  mm  m  m  m 

xxi  ui  U.CJ.  ou  wcx^iio  ux  xii  gjCriicx  cxx  u.oxx.xzic:  y  uiic  cva.jLua.ocix  x  uiix  oxuii  xii  acxuxxaiicc 

with  given  definitions.  Recall  that  the  definition  of  the  trapezoidal  rule 


I  =  H 


+  f(A+H)  +  f(A42H) 


+ 


+ 


V.  -^  / 
2 


requires  the  end  terms  be  weighted  by  l/2,  and  the  interior  terms  require 
weights  of  unity.  In  EXAMPLE  6,  we  disposed  of  the  end  terms  by  treating  them 
separately;  l.e.,  we  explicitly  defined  them  and  constructed  general  definitions 
for  the  arguments  and  functional  values  corresponding  to  the  interior  terms. 

In  this  example  we  will  include  all  arguments  and  terms  in  general  definitions. 
Although  the  generalizations  are  relatively  simple  for  this  example,  the  con¬ 
cept  is  most  advantageous  when  general  definitions  are  complicated  and  mere 
repetitions  are  laborious. 

Consider  the  following  generalization: 


EVRM 


f(X)  -  . 

GO  TO 


*  C?  T 

u'vy'iiic  jjj-o.'w.c; 


I 


“TJ.  •  ^  _ _ J  0-1 _ 0-  _ f _ o-_  _ _ X _ j _ TmmivT  ( _ t _ j _ or* _ to _ ^  V  i _ x-,i- _ 

j.t  lis  fciSsuxiieu.  utifciu  prxui’  uu  eriuex'xug  Veva.xua.ufci  luneuxuu^^  a  iiab  uSAea 

on  a  prescribed  value.  The  function  is  then  evaluated  and  recorded  as  per 
definition  and  the  computer  then  goes  to  "someplace"  for  subsequent  instruct- 


Tv%  /-»  -1  -Pt  r  ■ 

Xii  CI.U.U.X  uxVi^ii  ov./  X  J.  ^  . 


t  i-r  +  V»  Y”  -f*/— t.tVv  -1  <-*  Vi  W  -f*!  f  ro  T  1  1  Ci 

1^  L/ii*:::  ox  x  wx  wiix’v-ii  v.»ii»:;  x  ucii^  oxv-ziicix  vcxxu.^ 


required,  we  can  specify  the  "someplace”  for  subsequent  instructions’.  That 
is,  we  can  let  "someplace"  take  on  a  specific  name  prior  to  entering  EVEN. 


An 


For  example^  suppose: 


Let  X  =  A 

Let  "someplace"  =  Place  1 
Go  to  evaluate  function 


First: 


This  is  appended  to  the  above  generalization  and  illustrated  schematically 
below 


When  the  computer  "arrives"  at  Place  1^  (since  someplace  has  existing  value 
Place  l)j  f(X)  has  been  evaluated  for  X  =  A  and  hence  is  available  for  what¬ 
ever  treatment  is  desired.  Suppose  then  at 

PLACE  1  :  Form  I  =  f(X)/2 

Let  "Someplace"  =  Place  2 
Go  to  generate  next  X 

where  at 

NEXTX  :  Form  X  =  X  +  H 

Go  to  evaluate  function 

The  flow-chart  is  augmented  to  reflect  the  above  as  follows: 


6l 


^  v*  ■?  r\v*  +  CiTro  Tno+T-nrr  +Viq  "Pi  tti  r»  +  ■> 

X  X  ^'w'X  O'w'  V  -1- 00.1.1.^  oxxv«  X  . 


ion  for  X  =  X  +  H,  "som^ 


,...c:iT\nor*c:i*  Vicac?  cavnc+T-nrr  t  ro  1  t»  ca 

xxxv«  XX^C?  C'X.XX^ 


Place  2,  hence  at  Place  2,  f(X)  corresponding  to  f(X+H)  is  available  for  sub¬ 
sequent  treatment.  Now,  at  Place  2,  "subsequent"  treatment  depends  on  whether 
f(X)  has  just  been  evaluated  for  X  =  B,  (for  if  it  has  we  want  to  apply  weight 
1/2,  otherwise  weight  unity;).  Hence  at 


PLACE  2. 


.M 


If  X  -  B  < 


otherwise,  form 


go  to  finish  ; 

1=1+  f(x) 


Go  to  next  X  . 


(Note  that  "someplace"  retains  the  undisturbed  existing  value,  (Place  2). 

Wlien  the  computer  "arrives"  at  finish,  the  existing  value  of  I  represents  the 
sum  of  the  first  N  terms.  Now  we  need  only  add  l/2  of  the  existing  f(x)  =  f(B) 
to  I  and  multiply  by  H  to  obtain  the  desired  solution.  Hence  at 

FINISH:  I  =  H  (I  +  f(X)/2) 

Print - etc. 


The  complete  flow-chart  which  includes  the  above  is  shown  in  FIGURE  l8,  the 
corresponding  program  and  results  are  shown  in  FIGURE  19* 
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FLOW  CHART  FOR  EXAMPLE  J. 
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We  have  in  our  previous  examples  made  direct  references  hy  identifying 
quantities  by  names,  i.e.,  we  constructed  symbolic  names  of  variables,  functions, 
statements,  etc.  It  Is  Important  to  distinguish  between  "names”  of  quantities 
and  the  quantities  themselves.  If  for  example, 


Y  =  1  21 

•*».  ~ 

Y  =  72 


wc  wx  A.  uc: 


Z  =  X  +  Y 


Z,  X  and  Y  are  names,  whereas  75=5^  1=5  and  72  are  the  quantities  or  one  may 
say  that  75*5^  1-3-  and  72  are  the  numerical  values  of  Z,  X,  and  Y. 

In  arithmetic  expressions  and  formulas  such  as  the  above  it  is  understood 
that  the  operations  are  to  be  performed  on  the  numerical  values  and  not  their 
names.  As  we  shall  illustrate  shortly,  in  some  statements  it  is  the  name  and 
not  the  numerical  value  that  is  inferred.  A  trivial  exanple  is  the  statement, 
GOTO(BEGIN).  In  general,  one  need  not  associate  a  numerical  value  with  the 
name  BEGIN,  (we  shall  leave  such  Interpretations  to  the  professional  pro¬ 
grammer!  ). 

In  this  example  we  introduce  the  SET  statement  which  illustrates  a  case 
where  the  name  and  not  the  value  is  inferred.  SET  statements  are  generally  of 

uii6  J.  ux  lii 

set(a  =  b) 

where: 

A  is  generally  the  symbolic  name  of  an  index; 

B  is  either  a  symbolic  name  or  an  explicit  integer. 


If  B  is  a  symbolic  name. 


f  V.  ^■4-  +  V»  ^ 

liw  U  UllC 


vct-nxc:  wx  XJ  y  . 


then  A  takes  on  a  value  that  is  the  name  B, 


If  B  is  an  explicit  integer,  then  A  takes  on  the  value  that  is  the 
explicit  integer. 
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K0nG0  ^  "tiiG  'v3,lu.0  'fciiS'b  A  "bSilcGs  on  (isj)snccs  on  whether  B  is  3.  sy iiiholi c 
name  or  whether  B  is  an  explicit  integer. 


EXAMPLES: 

SET(A  =  SAM)  ,  then  the  existing  value  of  A  is  SAM. 
SET(A  —  ^  then  the  existing  value  of  A  is 


Recall  that  in  the  previous  discussion  of  general  names,  a  name  that  contains 
the  special  character  comma  is  an  *indexeu.  name  and  indexed  names  depend  on 
the  existing  "value'*  of  the  index.  For  example,  XI, I  is  an  indexed  name. 

We  call  XI  the  base  name  and  I  an  index  name. 

If  I  has  existing  value  2,  then  XI, I  represents  the  name  X5. 

Similarly,  the  name  ,J  is  an  indexed  name. 

If  J  has  existing  value  PETE,  then  ,J  represents  the  name  PETE  or  we  may  say 
that  the  "value"  of  ,J  is  PETE. 


Hence,  indices  may  take  on  two  types  of  values: 

a. )  the  value  of  an  index  may  be  a  name  *  ; 

b. )  the  value  of  an  index  may  be  an  integer. 

(in  EXAMPLE  8  we  will  illustrate  an  example  where  the  value  of  an  index  is  an 
integer) . 


In  the  generalization 


EVEN 


i 

I  V-^/  - . 

GO  TO 
"Someplace" 


to  let  "someplace"  =  Place  1  prior  to  entering  the  evaluation  of  the  function 


*  In  general,  an  index  can  only  take  on  integer  values  since  all  names  are 
automatically  transformed  to  unique  Integers  dvuring  the  transformation  from 
the  FORAST  language  to  the  primitive  machine  language . 
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SET(SMPL 


PLACE  l) . 


with  X  =  we  write 


Then  following  the  definition  of  f(X)  we  write 

G0T0(,SMPL) 

Hence,  after  evaluation  of  f(A).  the  existing  value  of 

SMPL  is  PLACE  1 

conseq_uently  the  eomputer  goes  to  PLACE  1  for  subseq^uent  instructions.  Simi- 
larly  at  PLACE  1^  we  defined 


I  =  f(X)/2  and  followed  this  definition  with 
SET ( SMPL  =  PLACE  2) 

thereby  establishing  the  existing  value  of  SMPL  as  PLACE  2.  Throughout  the 
remainder  of  the  computation,  for  the  existing  A,  B  and  N,  SMPL  retains  its 
value  of  PLACE  2  so  that  all  exits  from  the  evaluation  of  the  function  are 
directed  to  PLACE  2. 
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EXAMPLE  8. 


This  example  illustrates  a  convenient  use  of  an  index  which  takes  on 
prescribed  integer  values.  We  will  also  introduce  a  COUNT  statement,  a 
BLOC  statement,  anothf^r  form  of  the  general  READ  statement  and  a  GOTO  (next 
problem)  statement. 

GIVEN: 

recorded  on  cards  in  standard  form,  six  per  card,  i  =  1,2, ...,100 

REQUIRED: 


100 


+  "P  v»tti 

u.  1  L/’j.  ill  • 


START  PARSUM 


FLOW  CHART  FOR  EXAMPLE  8 
FIGURE  20. 
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On  the  flow  chart  of  FIGURE  20,  we  have  indicated  that  we  want  the 
computer  to  START  by  reading  and  recording  the  100  niimerical  values  correspond 
ing  to  the  given  .  We  assimie  that  they  are  recorded  on  cards  in  standard 
form,  6  numerical  values  per  card.  To  instruct  the  computer  to  do  this  we 
write 

read(ioo)nos.at(xi)^ 

Ths  In‘t0g0r  0nclQS0ci  in  p3.r0n'th0S0S  i2i]ni0d.i3.'tely  READ  sp0Qifi0s  1^110 

niimber  of  values  that  are  to  be  read  and  recorded.  (Actually,  one  may  write 
the  explicit  integer  as  above,  or  the  name  of  an  existing  integer).  The 
name  enclosed  in  parentheses  immediately  after  AT  specifies  where  they  are 
to  be  recorded  internally. 

To  insure  that  there  will  be  no  conflict  in  identifying  the  100  X's,  we 
reserve  internal  space  by  writing  a  BLOC  statement 

BL0C(X1  -  XIOO)^ 

The  BLOC  statement  not  only  reserves  100  consecutive  spaces  for  the  X's,  it 
also  establishes  X2,  X3,  ...,  X99  as  valid  names  associated  with  the  X's. 
(This  may  seem  trivial,  however,  it  is  necessary  since  names  such  as  XA,  X' , 
X.,  etc.  are  also  permissible  names'.).  Since  BLOC  is  an  order  type  (O.T.) 
word  it  may  only  appear  in  columns  7-10;  the  bloc  definitions  themselves  are 
recorded  commencing  in  column  11.  Also,  since  BLOC  is  an  order  type  (o.T.) 
word,  it  is  not  transformed  into  machine  instructions,  hence  it  is  generally 
placed  after  the  PROB  card  and  before  the  START  card. 

To  obtain  the  sum  of  the  X' s  we  could  write  a  statement 

S  =  XI  +  X2  +  X5  +  etc . , 

however;  i  may  be  large  (several  thousand)  and  consequently  we  generalize  so 
that  we  do  not  have  to  write  and  identify  each  and  every  term.  We  will  genera 
lize  by  computing  "partial  sums",  generally  denoted  by  S^.  .  We  first  SET(l=0) 
and  define  the  corresponding  partial  sum,  S^  =  0.  The  general  partial  s\im  is 
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defined  in  terms  uf  the  previous  partial  sum;  i.e.. 


"l.l 


=  S. 
1 


+  ‘X.  , 
i+1 


When  our  ultimate  goal  is  to  obtain  the  final  sum^  we  need  not  even  distinguish 
between  the  (i+l)st  and  the  ith  partial  sum,  hence  we  may  write 

s  =  s  . 


because  the  existing  value  of  S  is  used  on  the  right  of  the  equality  before 
the  resulting  value  is  assigned  to  the  name  on  the  left  of  the  equality.  We 
have  arbitrarily  labelled  the  general  partial  sum  definition,  PARSUM.  The 
program,  corresponding  to  the  portion  of  the  flow  chart  which  has  been  dis¬ 
cussed  thus  far,  is  as  follows; 


START 

PARSUM 


PROS 

BLOC 


C9O6  M.J.  ROMANELLI  L5IO7  EXAMPLE  8 
(XI  -  XlOO) 

read(ioo)nos.at(xi)^  SET(I=0)^  S  =  0 
S  =  S  +  XI, I 


When  the  computer  encounters  the  PARSUM  statement  for  the  first  time,  the 
existing  value  of  S  is  zero  and  1=0,  hence  after  executing  the  PARSUM  state¬ 
ment  for  the  first  time,  the  computer  obtains 

S  =  0  +  XI  =  XI  . 


Since  we  have  generalized,  we  could  obtain  o\ir  desired  solution  if  we  could 
instruct  the  computer  to  execute  the  PARSUM  statement  exactly  100  times ,  cor¬ 
responding  to  I  =  0,1,2, ... ,99.  The  general  COUNT  statement  was  designed 
specifically  for  such  situations'.  The  COUNT  statement  is  generally  written  in 
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Uiic  ±  lii 


COIMT  (A)  IN  (b)  goto  (C)  ^ 


where: 

A  is  the  name  of  an  explicit  integer  or  name  of  an  integer; 
B  is  the  name  of  an  integer; 

C  is  the  name  of  a  place. 


Firsts  the  index,  whose  name  appears  in  the  parenthesis  after  IN,  ’is  auto¬ 
matically  increased  by  unity. 

Next,  the  existing  value  of  the  B  index  is  compared  with  the  integer  indicated 
in  parenthesis  after  COUNT,  (or  with  the  existing  value  of  the  index  whose  name 
appears  in  parenthesis  after  COUNT),  i.e.  the  existing  value  of  B  is  compared 
with  the  existing  value  of  A. 

Finally,  if  the  existing  value  of  B  is  less  than  the  existing  value  of  A,  then 
the  computer  goes  to  C  for  subsequent  instructions;  otherwise,  i.e.,  if  the 
existing  value  of  B  is  >  or  =  the  existing  value  of  A,  the  computer  goes  to 
the  next  statement,  i.e.,  the  statement  immediately  following  this  COUNT  state¬ 
ment  . 


Again,  to  generalize,  if  the  initial  existing  value  of  B  is  the  integer  b  and 
the  existing  value  of  A  is  the  integer  a,  then  the  instructions  at  C  will  be 
executed  exactly 

(a  -  b)  times  . 


For  our  example,  we  write 

count(ioo)in(i)goto(parsum)  5^ 

Since,  we  initially  set  I  =  0  so  that  XI, I  initially  represents  the  first  X, 
we  have  a  =  100  and  b  =  0,  hence  the  PARSUM  statement  will  be  executed 
(lOO-O)  times.  Of  equal  importance  is  the  fact  that  the  index  i  will  auto¬ 
matically  take  on  the  desired  values  1,2,5^ =■■ >100*  have  emphasized  the 
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fact  on  the  flow  chart  hy  writing  i  =  i  +  1  in  the  box  just  prior  to  the 
determination  as  to  whether  i  < 100  ?  Note  also  that  even  though  i  will 
take  on  the  value  100,  the  PARSUM  statement  will  not  be  executed  corresponding 
to  i  =  100,  since  as  indicated,  i  is  advanced  by  unity  inmediately  after  the 
execution  of  the  PARSUM  statement  and  just  prioi’  to  the  determination  of 
i  <100  ?  Hence,  the  single  COUNT  statement  corresponds  to  the  two  boxes  of 
the  flow  chart  indicated  below: 


Hence,  immediately  after  the  PARSUM  statement,  we  write  the  COUNT  statement. 

The  COUNT  statement  is  followed  by  the  PRINT  statement  and  then  the  "new" 
statement 

goto(n.prob) 

is  written.  We  use  this  latter  statement  to  signal  to  the  computer  the  end 
of  the  execution  of  our  problem,  and  it  directs  the  computer  to  the  reading  of 
the  next  problem  residing  in  the  input  device. 

The  program,  input  and  output  corresponding  to  EXAMPLE  8  is  listed  in  FIGURE  21. 
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PK08  C906  M.J.  ROMANELLI  45107  EXAMPLE  8 
BLOC(XI-XIOO) 

START  READ{  100)N0S.AT(X1)2  S=Ot  SETU=0) 

PARSUM  S=S+X1,I 

COUNT ( 100) IN ( I ) GOTO (PARSUM) 

PRINT<S  =  >S%  GOTO{N.PROB) 


END 

10000000 

GOTO! START ) 

01  20000000 

01 

30000000 

01 

40000000 

01 

50000000 

01 

60000000 

01 

70000000 

01 

80000000 

01 

90000000 

01 

10000000 

02 

1  1000000 

02 

12000000 

02 

13000000 

02 

14000000 

02 

15000000 

02 

16000000 

02 

17000000 

02 

18000000 

02 

19000000 

02 

20000000 

02 

21000000 

02 

22000000 

02 

23000000 

02 

24000000 

02 

25000000 

02 

26000000 

02 

27000000 

02 

28000000 

02 

29000000 

02 

30000000 

02 

31000000 

02 

32000000 

02 

33000000 

02 

34000000 

02 

35000000 

02 

36000000 

02 

37000000 

02 

38000000 

02 

39000000 

02 

40000000 

02 

41000000 

02 

42000000 

02 

43000000 

02 

44000000 

02 

45000000 

02 

46000000 

02 

47000000 

02 

48000000 

02 

49000000 

02 

50000000 

02 

51000000 

02 

52000000 

02 

53000000 

02 

54000000 

02 

55000000 

02 

56000000 

02 

57000000 

02 

58000000 

02 

59000000 

02 

60000000 

02 

61000000 

02 

62000000 

02 

63000000 

02 

64000000 

02 

65000000 

02 

66000000 

02 

67000000 

02 

68000000 

02 

69000000 

02 

70000000 

02 

71000000 

02 

72000000 

02 

73000000 

02 

74000000 

02 

75000000 

02 

76000000 

02 

77000000 

02 

78000000 

02 

79000000 

02 

80000000 

02 

81000000 

02 

82000000 

02 

83000000 

02 

84000000 

02 

85000000 

02 

86000000 

02 

87000000 

02 

88000000 

02 

89000000 

02 

90000000 

02 

91000000 

02 

92000000 

02 

93000000 

02 

94000000 

02 

95000000 

02 

96000000 

02 

97000000 

02 

98000000 

02 

99000000 

02 

10000000 

03 

1 

1 

2 

3 

4 

5 

6 
7 
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S 


50500000  4 


FIGUEE  21. 


0000001 


EXAMPLE  9. 


This  example  illustrates  the  MOVE  and  INC  (increment)  statements,  and 
the  SYN  ( sjTiQnym)  order  type . 


GIVEN: 


and  f(X^)  recorded  on  cards  in  standard  form,  i  =  1, 2, 3; • • • ; 20. 
Ap;c;uTri(=>  t.Vi(=>  Y  nn  f QUT  COnseCUtive  CS'cds  .  ■'^lues  nn  each  nf  the 

i 

first  three  cards,  and  two  values  on  the  fourth  card.  Assixme 
the  corresponding  f(X^)  on  the  next  four  cards,  again  6  values 
on  each  of  the  first  three  cards  and  two  values  on  the  last  card. 

Assume  that  the  X.  are  not  in  monotone  sequence. 

1  - 


REQUIRED :  Obtain  a  "re-arranged"  table  of  the  X_.  and  corresponding  f(X^) 


X.  <  for  all  i. 

1  =  1+1 


Print  both  tables;  i.e.,  the  original  and  re-arranged  tables 
with  the  numerical  entries  in  standard  form.  In  particular, 
print  a  table  in  the  following  form: 


ORIGINAL 


REARRANGED 


XI  FI 


XI  FI 


Since  we  will  have  to  instruct  the  computer  to  determine  if  each 

we  will  first  provide  for  getting  the  given  table  of  numerical  values  recorded 

internally.  This  can  be  accomplished  with  the  following  READ  statements: 


READ ( 20 ) NOS . AT ( XI 
READ  ( 20  )N0S  .  AT  ( FI )  5^ 
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Alternatively,  we  can  provide  for  the  internal  recording  of  both  the  and 
corresponding  f(X^)  with  the  following  READ  statement: 

READ(U4)N0S.AT(X1)  ^ 


With  this  single  READ  statement,  we  are  instructing  the  computers  to  include 
in  the  44  numerical  values,  the  foiir  "extraneous"  X^,  (i  =  21,22,25,24),  of 
the  4th  card.  (Recall  that  only  two  X_. ,  (i  =  19, 20)  are  recorded  on  the  4th 
card).  In  the  absence  of  a  format,  the  computers  assimie  6  standard  values 
per  card  until  the  READ  (or  PRINT)  statement  is  fulfilled. 

To  reserve  space  in  Internal  storage  for  the  X^  and  f(X^)  we  write  the 
BLOC  statement 


BLOC (XI  -  X44) 


For  future  references  to  the  f(X^)  we  write  a  SYN  statement 

SYT^Fl  =  X25) 

In  SYR  statements,  the  equated  names  are  assigned  identical  internal  names 
(identical  integers)  and  hence  references  to  any  of  the  equated  names  refer 
to  the  same  quantity.  Again,  recall  that  the  BLOC  statement  not  only  reserves 
space  in  internal  storage,  it  also  establishes  X1,X2,X5, . • . ,X44  as  valid  names. 
The  above  SYIi  statement  establishes  FI  and  X25  as  synonomous  names.  Without 
the  SYN  statement  we  could  refer  to  the  f(X^)  in  many  ways,  for  example. 


or 


X24,I  where  I  =  1,2, ... ,20 
XI, J  where  J  =  24,25, ... ,45- 


However,  with  the  SYN  statement  we  can  refer  to  the  f  (X^.  )  as 


FI, I 


thereby  retaining  the  symbolic  identification  of  the  given  data.  Further,  if 
we  refer  to  the  X^  as  XI, I  and  if  we  refer  to  the  f(X^)  as  FI, I  we  need  only 
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"manipulate"  the  numerical  value  for  I  to  make  references  to  and  the 
corresponding  f(X^.  ). 

Our  problem  requires  that  both  the  original  and  re-arranged  forms  of  the 
table  are  to  be  printed.  To  have  both  forms  available  in  internal  storage 
for  printing,  we  can  have  the  given  table  recorded  in  two  areas.  In  one  area 
we  can  retain  the  original  form,  in  another  area  we  can  interchange  any  values 
required  to  obtain  the  desired  re-arranged  form.  To  reserve  space  for  a 
second  area  we  can  au^ent  the  BLOC  statement;  i.e.,  we  simply  add  another 
parenthesized  expression  which  defines  the  second  area.  Let  us  arbitrarily 
denote  the  second  area  as  (X'l  -  X'hU).  The  BLOC  statement  which  reserves 
space  for  both  areas  is  then  of  the  form 

BLOG (XI  -  X44)(X‘l  -  X»44)^ 

(Alternatively,  one  could  define  a  BLOC  (Xl  -  X88)  and  an  accompanying 
SYN  (X'l  =  X45)  to  reserve  space  for  both  tables  and  identification  of  the 
second  area  begirming  at  X45). 

Now,  the  single  READ  statement  will  provide  for  the  recording  of  the 
given  table  in  the  XI  -  X44  area.  To  obtain  a  copy  of  the  given  table  in  the 
X'l  -  X'44  area,  we  use  the  convenient  MOVE  statement.  The  MOVE  statement  is 
generally  of  the  form 

move(a)nos.frcm(b)to(c)^ 

where : 

A  is  either  the  name  of  an  integer  or  an  explicit  integer  which 
defines  the  number  of  quantities  that  are  to  be  moved; 

B  is  the  name  of  the  first  quantity  in  a  source  area; 

C  is  the  name  of  the  first  quantity  in  a  destined  area. 

The  quantities  are  retained  in  the  source  area  and  hence  after  the  MOVE  state¬ 
ment  has  been  executed  the  quantities  exist  in  both  areas.  In  the  above  MOVE 
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statement  it  is  assumed  that  the  (Quantities  in  both  areas  will  be  uniJ 


UJ.  lii-LJ 


spaced, one  unit  apart.  Uniform  spacings  other  than  one  unit  may  be  specified 
for  the  source  and(or)  destined  area  by  folloving  the  symbolic  name  with  a  / 
and  the  integer  defining  the  spacing.  For  example, 


move(n)nos.from(a/5)to(b/5)^ 


indicates  that  the  N  quantities  in  the  source  area  are  three  units  apart, 
the  N  quantities  in  the  destined  area  are  to  be  spaced  five  units  apart. 

For  oiiT  example  we  write 

MOVE  (  UU  )  NOS .  FRQt4(Xl  )T0  ( X  ‘  1 )  ^ 

Ordinarily  we  plan  a  solution  by  drawing  a  flow  chart.  In  this  example, 
we  delayed  the  drawing  of  the  chart  until  the  capabilities  of  the  new  state¬ 
ments  were  discussed  so  that  the  novice  obtains  some  insight  as  to  "how"  the 
flow  chart  solution  can  easily  be  expressed  in  the  FORAST  language.  The 
complete  flow  chart  for  a  solution  is  given  in  FIGURE  22. 
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START 


IffiAD: 

X!  =  X. 

1  1 

X.  f(X.) 

f'(Xj_)  =  f(X.) 

i 

1 

1  =  1,2, ...,20 

i  =  1,2,. ..,20 

Yes 


’'jn  <  1 

9 

i+2  • 

■ 

INTERCHANGE; 


j+2 


j+1 

^j+2 


f 

f 


j+1 

j+2 


=  t 


t 


TVrr\4*£i  4*Vic>4*  +Vic»  -Pt  ■k»c+  Vif^v^c*  A'f*  "fl  /^"P  "fTTriTTRTT  OO 

XV  w  ^liu.  ^  ^liv,^  x^x>3«.<>  r*^  U  o  Wx  OllC  X  ^  w  r?  Lril'-A>x  v/  x  s^x  vxj  r  f 

indicate  the  recording  of  the  given  data  in  two  areas  of  internal  storage, 
the  primed  and  imprimed  areas.  The  third  box  indicates  that  the  index  i 
IS  to  talce  on  valne  sero.  The  next  box  labelled  TSTFWD  (test  fon-.'ard) 
indicates  the  monotone  test.  If  the  monotone  condition  is  satisfied,  the 
index  i  is  advanced  by  one  and  thence  a  determination  is  made  as  to  whether 
all  of  the  X's  in  the  table  have  been  tested.  If  all  have  been  tested,  the 
existing  value  of  i  is  19  and  hence  we  indicate  the  printing  of  the  tables. 
If  all  X's  have  not  been  tested,  (i  <  19);  hence  we  want  the  computer  to  test 

+  >n^'V+  -rxQ  *1  v» 

If  the  monotone  condition  is  not  satisfied,  we  want  to  interchange 
X^^^  and  ^^^2’  corresponding  functional  values.  This  interchange  will 

guarantee  that  this  pair  now  satisfies  the  monotone  condition,  but  once  an 
interchange  is  made,  we  have  no  guarantee  that  the  monotone  condition  prevails 
for  all  X's  that  have  previously  been  tested.  Hence,  we  must  test  backwards 
until  one  pair  of  X's  does  satisfy  the  monotone  condition.  To  perform  this 
backward  test,  note  that  we  define  another  index  j  which  initially  takes  on 
the  existing  value  of  i.  Next,  if  j  =  0  no  further  backward  testing  is 
necessary,  hence  we  indicate  that  forward  testing  is  to  be  resumed.  If  J 
is  not  zero,  the  existing  value  of  j  is  diminished  by  one  and  the  backward 
test  continued.  Note  that  an  interchange  is  indicated  any  time  the  monotone 
condition  is  not  satisfied.  As  soon  as  the  condition  is  satisfied  in  the 
backward  testing,  the  forward  testing  is  resumed. 

The  program,  input  and  output  for  this  example  are  listed  in  FIGURE  25. 
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BL0C{X1-X44)X' 1-X*441  2 

SYN  (E1=X25)  3 

START  KEAD{44)N0S.AT{Xl)  4 

MOVE{44)NOS.FROM(X1)TO{X*1)  5 

SET(I=0)  6 

TSTFWD  1F( Xl,I<=X2,I  )GOTOCCI  )  7 

T=X1,I2  Xl,I=X2,I%  X2,I=T  8 

T  =  Fl,U  FI,  l  =  Fl,  (  I+l  J*  F1,(I+1)  =  T*  J=I  9 

IF-INT( J=0)G0T0(TSTFWO)«  INC{J=J-1)  10 

1F(X1,J<=X2,J)G0T0{TSTFWD)*  11 

T=X1,J«  X1,J=X2,J*  X2,J=T  12 

T  =  F1,J!?  F1,J=F1,(  J  +  l)%  F1,(J  +  1)=T?  GOTO(CJ)  13 

Cl  C0U.NT(19)IN{I  )GOTOiTSTFWO)*  14 

PHINT<  ORIGINAL  REARRANGED>  15 

CX3 

^  PRINT<  XI  FI  XI  FI>%  SET! 1=0)  16 

WORK  PRINT ( X' 1 , I )X'25, I) XI , I )F1, I )?  COUNT  I  2  0 ) IN ( I ) GOTO [ WORK )  17 

GOTO(N.PROB)  18 

END  GOTO! START)  19 


93000000  00  95000000  00  98000000  00  10000000  01  12000000  01  12500000  01 

17000000  01  17500000  01  18000000  01  18500000  01  19000000  01  19300000  01 

13000000  01  13500000  01  14000000  01  15000000  01  16000000  01  16500000  01 

19800000  01  20000000  01 

59783000  00  58168000  00  55702000  00  54030000  00  36236000  00  31532000  00 

-12884000  00-17825000  00-22720000  00-27559000  00-32329000  00-35153000  00 
26750000  00  21901000  00  16997000  00  70740000-01-29200000-01-79120000-01 
-39788000  00-41615000  00 
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original 

REARRANGED 

0000001 

XI 

FI 

XI 

FI 

0000002 

93000000 

59783000 

93000000 

59783000 

0000003 

95000000 

58168000 

95000000 

58168000 

0000004 

98000000 

55702000 

98000000 

55702000 

0000005 

10000000 

1  54030000 

10000000 

1  54030000 

0000006 

12000000 

1  36236000 

12000000 

1  36236000 

0000007 

12500000 

1  31532000 

12500000 

1  31532000 

0000008 

17000000 

1-12884000 

13000000 

1  26750000 

0000009 

17500000 

1-17825000 

13500000 

1  21901000 

0000010 

18000000 

1-22720000 

14000000 

1  16997000 

0000011 

18500000 

1-27559000 

15000000 

1  70740000- 

■01 

0000012 

19000000 

1-32329000 

16000000 

1-29200000- 

01 

0000013 

19300000 

1-35153000 

16500000 

1-79120000- 

01 

0000014 

13000000 

1  26750000 

17000000 

1-12884000 

0000015 

13500000 

1  21901000 

17500000 

1-17825000 

0000016 

14000000 

1  16997000 

18000000 

1-22720000 

0000017 

15000000 

1  70740000-01 

18500000 

1-27559000 

0000018 

16000000 

1-29200000-01 

19000000 

1-32329000 

0000019 

16500000 

1-79120000-01 

19300000 

1-35153000 

0000020 

19800000 

1-39788000 

19800000 

1-39788000 

0000021 

20000000 

1-41615000 

20000000 

1-41615000 

0000022 

FIGURE  25. 


The  program  contains  two  new  concepts  which  have  not  been  discussed. 
Note  first  the  name 

FI, (I  +  1) 


on  line  9-  This  is  used  as  a  reference  for  f(X^_^2)-  above  is  an 

illustration  of  an  indexed  name  of  the  form 


A.C'R  +  n'l 


where  A  is  a  defined  base  name,  B  is  the  name  of  an  index 
and  n  is  an  integer.  The  resulting  name  is  obtained  by  increasing  (or 
decreasing)  the  existing  value  of  B  by  n  and  adding  this  result  to  the  base 
name  A.  That  is,  the  parenthesized  expression  is  evaluated  and  the  resuJLting 
integer  is  added  to  A  to  obtain  the  resulting  name.  If  for  example,  the 
existing  value  of  I  is  7,  then 


Fl,(l  +  l)  represents  the  name  F9 
Fl,(l  -  3)  represents  the  name  F^. 


Next,  on  line  10  is  an  illustration  of  the  INC  statement  of  the  form 

INC(A  =  a  +  n) 

where  A  is  the  name  of  an  index  and  n  is  an  explicit  integer.  This  statement 
is  generally  used  to  Increase  (or  decrease)  index  values  by  constant  amounts. 
(More  general  means  for  operating  on  integer  values  will  be  illustrated  in 
later  examples)  For  our  example  we  wrote 

INC(J  =  J  -  l) 


which  told  the  compiler  to  dimish  the  value  of  J  by  one. 


SXAI-IPLE  10. 


This  example  illustrates  the  LIST  order-type  word  which  is  used  for 
two  purposes: 

a. )  It  is  used  by  the  programmer  to  obtain  a  listing  which  shows  the 
one  to  one  correspondence  between  the  symbolic  names  in  his  program  and  unique 
integers.  The  unique  integers  represent  the  absolute  machine  names  of  the 
internal  storage  units.  Hence  this  listing  shows  the  programmer  the  explicit 
storage  assignment  by  the  FORAST  compiler.  Since  the  compiler  does  not  check 
for  all  possible  violations  and  conflicts,  a  study  of  the  listing  can  reveal 
conflicts  of  storage  assignment  or  other  programming  errors. 

b.  )  It  is  used  by  the  professional  programmer  to  obtain  the  absolute 
machine  language  corresponding  to  his  symbolic  FORAST  language.  We  will 
indicate  how  this  is  obtained  and  list  the  absolute  machine  language  of 
EXAMPLE  10;  however,  we  leave  the  interpretation  to  the  professional  programmer. 


During  the  transformation  from  the  symbolic  FORAST  language  to  the  absolute 
machine  language,  indications  of  errors  will  be  printed  to  inform  the  pro¬ 
grammer  of  the  errors.  Further,  the  listing  (dictionary)  referred  to  under 
a.)  above  will  be  produced  and  no  attempt  will  be  made  to  execute  the  program. 
The  novice  will  find  on  many  occassions  that  even  though  no  errors  were  detect¬ 
ed  during  compilation,  his  program  does  not  run  correctly.  In  such  a  case,  no 
dictionary  is  obtained  and  to  begin  to  determine  a  source  of  error  (or  errors) 


it  is  advisable  to  obtain  the  dictionary.  To  obtain  the  dictionary,  we  simply 
record  LIST  in  columns  7  through  10  of  a  card  and  place  this  card  in  front  of 
the  END  card.  This  indicates  that  the  dictionary  is  desired  and  will  be  pro¬ 
duced  whether  or  not  any  errors  are  detected  during  compilation.  (When  the 
absolute  machine  language  is  desired  in  addition  to  the  dictionary,  S.CODE  is 
recorded  beginning  in  column  11  of  the  LIST  card  referred  to  above ) . 


We  have  inserted  a  LIST  card  in  the  program  of  EXAI»IPLE  10. 
are  listed  in  FIGURE  24. 


The  results 
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PROB  C906  M.J.  ROMANELLI  45107 
BLOC(  XI-X44)  X'  ]L-X'44) 

SYN  (F1=X25) 


EXAMPLE  10 


23  JULY, 63  PAGE 


1 
1 

2 

3 

START  READ(44)N0S.AT(X1)  4 

M0VE(44)N0S.FR0M(X1)T0(X*1)  5 

SETII=0)  6 

TSTFWO  IF{X1 ,I<  =  X2,I )GOTO<CI  )  7 

T=X1,I%  X1,I=X2,1X  X2,I=T  8 

T  =  F1,U  F  1  ,  I  =  F1  ,  ll  +  l  )  «  F1,(H-1)  =  T*  J=I  9 

CJ  IF-If\IT(  J=0)GCIT0(  TSTFWO)^  INC(J=J-1)  10 

1F{X1, J<=X2,J)GOTO(TSTFWO)*  11 

T  =  X1,J?  X1,J=:X2,J*  X2,J=T  12 

T  =  F1,J«  F1,J=:F1,(  J  +  l)%  F1,(J<-1)  =  T%:  GOTO(CJ)  13 

Cl  COUNT( 19) IN(I )GOIO{TSTFWO)?  14 

PRINT<  ORIGINAL  REARRANGED>  15 

PRINT<  XI  FI  XI  FI>3;  SET{I  =  0)  16 

WORK  PRINT(X'1,I)X‘25„I)X1,I)F1,I)*  COUNT  I  20 ) I N ( I ) GOTO ( WORK )  17 

GOTO (N. PROB)  18 

LIST  S.COOE  18.1 

END  GOTO( START)  19 


93000000  00  95000000  00  98000000  00  10000000  01  12000000  01  12500000  01 

17000000  01  17500000  01  18000000  01  18500000  01  19000000  01  19300000  01 

13000000  01  13500000  01  14000000  01  15000000  01  16000000  01  16500000  01 

19800000  U1  20000000  01 

59783000  00  58168000  00  55702000  00  54030000  00  36236000  00  31532000  00 

-12884000  00-17825000  00-22720000  00-27559000  00-32329000  00-35153000  00 
26750000  00  21901000  00  16997000  00  70740000-0 1-29200000-0 1-79 1 20000-0 1 


-39788000  00-41615000  00 


CD 

VJl 


MAY. 23*63  BRLESC  FORAST  F62 


PROB 

C906 

M.J. 

ROMANELLI 

45107 

EXAMPLE 

10 

Cl 

118  L 

N.PROB 

N70  F 

XI 

12S 

BSM 

CJ 

lOF  L 

START 

100  L 

X44 

001 

B 

FI 

143 

SM 

T 

157  M 

X*  1 

158 

B  M 

1 

OOK 

I 

TSTFWD 

106  L 

X'44 

001 

B 

J 

00  S 

I 

WORK 

126  L 

?INDEX 

OON 

MAY. 23*63  BRLESC  FORAST  F62 

PROB  C906  M.J.  ROMANELLI  45107  EXAMPLE  10 

IOC)  05  6K  0  70  OF  12S  2N  1  01  8000  NOOOIOOOO 

104  13  8000  2N  103  0128000  0  0  622812N28L2S  118 

108  KK  02812N2812S  KK,  0  1572812N  KK  028143  157 

ION  KK  0  15728144  KK,  0  K  S  S6  S  50  106 

110  622N12N2N12S  106  KK  02N12S  157  KK  02N12N2N12S 

114  KK  02N143  157  KK  02N1442N143  KK  0  1.572N144 

118  1328000  13  106  05  6N  0  71  IF  0  0  9K9 

IIN  IF  0  0  0  IF  29551K694655J5  1F53J00  0  0 

120  IF  0  37640000  IF  0  016640  IF  0  0  37 

124  IF  1667J0000  0  0128000  0  0  05  6N  0  71 

128  1028143  0  0  1328000  14  126  04  0  0  N70 

MAY. 23, 63  BRLESC  FORAST  F62 

PROB  C906  M.J.  ROMANELLI  45107  EXAMPLE  10 

ORIGINAL  REARRANGED 

XI  FI  XI  FI 

93000000  59783000  93000000  59783000 

95000000  58168000  95000000  58168000 

98000000  55702000  98000000  55702000 

lOOOOOOO  1  54030000  10000000  1  54030000 

12000000  1  36236000  12000000  1  36236000 

12500000  1  31532000  12500000  1  31532000 

17000000  1-12884000  13000000  1  26750000 

17500000  1-17825000  13500000  1  21901000 

18000000  1-22720000  14000000  1  16997000 

18500000  1-27559000  15000000  1  70740000-01 

19000000  1-32329000  16000000  1-29200000-01 

19300000  1-35153000  16500000  1-79120000-01 

13000000  1  26750000  17000000  1-12884000 

13500000  1  21901000  17500000  1-17825000 

14000000  1  16997000  18000000  1-22720000 

15000000  1  70740000-01  1B500000  1-27559000 

16000000  1-29200000-01  19000000  1-32329000 

16500000  1-79120000-01  19300000  1-35153000 

19800000  1-39788000  19800000  1-39788000 

20000000  1-41615000  20000000  1-41615000 


*NOS.  OSO 
*SUBS.  N70 


KN  812S  0  0158 
KK  02812S  157 
KK  02814428143 
022LLLL  0  0 
KK  0  1572M12N 
04  0  0  lOF 
1F6576654630  0 
05  6N  0  71 
1F64000  0  0 
1O281582017O2812S 


0000001 

0000002 

0000003 

0000004 

0000005 

0000006 

0000007 

0000008 

0000009 

0000010 

0000011 

0000012 

0000013 

0000014 

0000015 

0000016 

0000017 

0000018 

0000019 

0000020 

0000021 

0000022 


FIGURE  2k. 


Note  first  that  all  of  the  symholic  names  used  in  the  program 
in  "alphabet ical"  order.  (Numerals  in  general  are  "less  than"  any  alphabetic 
character  and  appear  in  the  "alphabetical  sequence"  before  any  alphabetic 


a-r*  . 


Immediately  to  the  right  of  a  symbolic  name  is  the  unique  integer 


assigned  by  the  compiler.  This  integer  is  recorded  in  the  sexadecimal  number 
system^  (base  l6).  The  characters  in  this  system  which  correspond  to  the 
decim.al  ten^  eleven^ ...  ^f  if  teen  are  respectively  K^S.N.J.F  and  L.  Hence  the 
integer  listed  corresponding  to  the  symbolic  name 


XI 


IS 


1^ 


which  represents  the  decimal  integer 

1  X  (l6)^  +  2x(i6)^  +  llx(l6)°  =  256  4  52  +  11  =  299- 

The  integer  (absolute  machine  name)  is  generally  followed  by  one  or  more  of 
the  alphabetic  characters,  A,B,F,I,L,M,R,S  or  U.  Only  a  few  of  the  meanings 
are  listed  below,  the  meanings  of  each  are  given  in  j^l^j  , 

B  indicates  a  BLOC  name  and  opposite  the  end  name  of  the  BLOC 
is  an  integer  which  represents  the  uniform  spacing  between 
the  elements  of  the  BLOC; 

F  indicates  the  nsime  of  a  function; 

I  indicates  that  the  name  is  an  index  name; 

L  indicates  a  name  that  was  recorded  in  the  LOCATION  columns, 
usually  a  statement  name; 

S  indicates  that  the  name  was  defined  in  a  SYN  statement. 

U  indicates  that  the  name  appeared  only  once  in  the  program 
and  hence  this  single  reference  may  indicate  a  programming 
or  punching  error. 


o 


Ob 


As  an  illustration  of  an  undetected  error,  suppose  on  line  J  of  EXAMPLE  10 
one  wrote  GOTO(Cl)  and  on  line  l4  in  the  LOCATION  column  one  wrote  Cl.  Note 
that 

Cl  ^  Cl 

and  both  Cl  and  Cl  are  distinct  and  valid  names.  Both  would  appear  in  the 
dictionary  followed  by  the  letter  U  to  indicate  single  references .  The 
obvious  error  is  that  Cl  and  Cl  were  not  designed  as  distinct  references, 
and  either  the  1  should  be  changed  to  I  or  the  I  to  1.  The  last  three  entries 
in  the  listing 

^  INDEX 
i  NOS. 
i  SUBS. 

correspond  to  the  extent  of  the  storage  used  for  indices,  constants  and 
functions.  The  integer  opposite  ^  INDEX  corresponds  to  the  integer  that 
would  be  assigned  to  the  next  index  encountered.  Similarly,  the  integer 
opposite  io  NOS.  corresponds  to  the  integer  that  would  be  assigned  to  the  next 
constant  encountered.  Finally,  the  functions  (subroutines)  are  allocated 
space  at  the  end  of  internal  storage  (largest  integers)  and  hence  the  next 
function  encountered  would  be  assigned  space  just  above  the  integer  opposite 
i  SUBS. 
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11. 


This  example  introduces  the  CLEAE  and  SETEA  statements  and  illustrates 
convenient  variations  of  the  COUKT,  READ  and  BLOC  statements. 

GIVEN:  X.  and  f(X. )  recorded  on  cards  in  standard  form,  i  =  1.2...., 20 

Assiame  the  X^.  are  in  monotone  sequence. 

REQUIRED :  Compute  and  print  a  table  of  forward  differences. 


..1 

a: 

1 


where 


j  =  1.  2,  5,  ^ 
i  =  1,  2,  5,  • • •  >  20 
A°  =  f 

and  all  undefined  A' s  =  0,  i.e., 

A^  =  0  for  j  =  1,  2,  5, 

Ag  =  0  for  j  =2,  ^ etc. 

Print  the  table  in  the  following  form: 


X,I 

F,I 

IDEL, I 

2DEL,I 

5DEL,I 

Ldel,i 

X. 

a_ 

f(x^) 

0 

0 

0 

0 

'^2 

f(X2) 

1 

^2 

0 

0 

0 

S 

^(x^) 

r\ 

KJ 

0 

f(\) 

.1 

% 

0 

f(x^) 

1 

A 

), 
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There  are  many  ways  that  one  can  obtain  the  solution  for  this  example. 

We  chose  the  method  indicated  in  the  flow-chart  of  FIGURE  25  to  illustrate 
a  variation  of  the  BLOC  statement  which  permits  us  to  conveniently  index  and 
reference  elements  in  a  two  dimiensional  array.  Illustrated  below  is  a  symbolic 
form  of  the  desired  two  dimensional  table. 


E 

1 

2 

3 

4 

3 

6 

1 

X. 

X 

X 

0 

0 

0 

0 

P 

x^ 

f  - 

0 

0 

0 

2 

2 

2 

Z 

J 

='3 

4 

0 

0 

2 

.3 

4 

0 

P 

5 

4 

5 

« 

^5 

"5 

« 

20 

^20 

"6 

4 

4 

4 

4 

We  have  labelled  the  rows  of  the  table  1  throu^  20  and  the  columns 
1  through  6.  Note  that  the  entries  in  the  first  two  columns  correspond  to  the 
given  data,  the  and  f(X^). 

To  reserve  internal  storage  space  for  this  array  we  write  the  BLOC  state¬ 
ment 

BL0C(E1,1  -  E20,6) 


The  first  name  in  the  enclosed  parenthesis,  El,l  defines  the  name  of  the 
first  element  in  the  array;  in  particular  it  defines  the  name  of  the  element 
in  the  first  row  and  first  column.  The  last  name  in  the  enclosed  parentheses, 
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E20,6  defines  "the  Iasi  element  in  the  array;  in  particular  EaO^o  defines  the 
element  in  the  20th  row  and  6th  column.  Hence ^  space  is  reserved  for  the  120 
elements  in  the  20  by  6  array.  Having  defined  the  array  with  the  BLOC  state- 


on'K c* 

xi.  (..f  ^  V.X  1 1  it;' 


ferences  to  any  element  in  the  array  can  be  made  by  writing 


E  followed  by  the  explicit  integers  which  define  the  row  and  column,  the  row 
and  column  integers  must  be  separated  by  a  comma.  For  example,  we  can  refer 
to  f(X^)  by  writing  E5,3  and  similarly  we  can  refer  to  by  writing  E5,l. 
Further,  any  of  these  references  may  be  used  as  a  base  name  in  a  general  index 
name.  For  example, 

E1,1,J 


is  a  valid  name  and  the  particular  element  referenced  depends  on  the  existing 
value  of  J . 

If  J  =  1,  E1,1,J  is  equivalent  to  the  name  El, 2  and  refers  to  element  f^. 
If  J  =7,  E1,1,J  is  equivalent  to  the  name  E2, 2  and  refers  to  element  fg. 
If  J  =  h,  E2,3;<I  is  equivalent  to  the  name  E3;l  and  refers  to  element  x^. 


The  programmer  must  bear  in  mind  that  successive  elements  in  a  row  are  spaced 
one  unit  apart  whereas  successive  elements  in  the  columns  are  spaced  six  unit 
apart.  Hence,  to  advance  references  from  one  element  to  another  the  index 
values  must  be  manipulated  accordingly.  The  flow-chart  for  this  example  is 
shown  in  FIGURE  25 • 
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START 


FLOW  CHART  FOR  EXAMPLE  11. 

FIGURE  25. 

12  5  U 

Note  first  the  indication  that  the  elements  A  ,  A  ,  A  and  A  are  to  be 

o  o’  o  o 

defined  as  zero  since  they  correspond  to  undefined  differences .  To  instruct 
the  computer  to  do  this  we  use  the  convenient  CLEAR  statement  which  is  of  the 
form 

clear(a)nos.at(b)^ 

where  A  is  the  name  of  an  index  or  an  explicit  integer,  and  B  is  the  name  of 
the  first  of  the  consecutive  quantities  that  are  to  be  cleared,  i.e.,  their 
numerical  values  are  set  to  zero.  If  A  is  the  name  of  an  index,  then  the 
current  value  of  that  index  specifies  the  number  of  quantities  that  are  to  be 
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cleared.  If  A  is  an  integer  then  the  integer  defines  the  number  of  quantities 
that  are  to  be  cleared.  It  is  assumed  that  the  quantities  are  uniformly  spaced 
in  internal  storage  and  in  the  absence  of  a  specified  spacing  it  is  taken  to 
be  unity.  If  a  spacing  other  than  unity  is  desired  it  is  specified  by  follow¬ 
ing  the  name  B  with  a  slashj  j ,  followed  by  the  Integer  which  defines  the 
spacing.  For  example,  if  we  want  to  clear  all  of  the  elements  in  the  4th 
column  we  would  wnrite 


CLEAR ( 20 ) NOS . AT ( El , 4/6 


For  our  example  we  write 


CLEAR ( 4 )N0S . AT(E1 , 5 


Next,  to  read  and  record  the  given  and  f(X^)  in  th  first  two  col’jmns 
of  the  reserved  array,  we  use  a  variation  of  the  READ  statement  where  we  can 
specify  that  the  uniform  spacing  of  the  quantities  to  be  recorded  is  not  unity. 
In  particular,  the  unifom  spacing  is  six,  hence  we  write 

READ(20)N0S.AT(E1,1/6)5^ 


which  provides  for  the  reading  and  recording  of  the  twenty  in  the  reserved 

area.  Similarly,  we  write 

READ(20)N0S.AT(E1,2/6)^ 

which  provides  for  the  reading  and  recording  of  the  twenty  f(X^*s)  in  the 
reserved  area. 


Now,  assuming  that  the  given  data  have  been  read  and  recorded  in  the 
reserved  area,  we  plan  to  generalize  the  difference  definition  so  that  by 
generating  prescribed  values  of  an  index  we  will  instruct  the  computer  to 
compute  and  record  the  differences  for  a  given  line.  Then  we  will  instruct  the 
computer  to  repeat  the  process  for  subsequent  lines. 


To  generalize,  note  that  each  difference  is  defined  as  the  difference  of 
two  neighboring  quantities  in  the  array,  namely  the  preceding  row  element  minus 
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its  preceding  coliomn  element.  In  particular. 


E2,3  =  E2,2  -  El, 2 

E5,5  =  E5,2  -  E2,2 

E4,3  =  E4,2  -  E5,2 

and  in  general  we  can  define 

E2,3,J  =  E2,2,j  -  E1,2,J 

to  obtain  all  of  the  desired  differences  provided  that  J  takes  on  prescribed 
values . 

First,  if  J  took  on  values  0,1,2,  and  3  we  would  obtain  results  for  the 
second  rowj  however,  only  the  first  difference  is  defined. 

Next,  if  J  took  on  values  6,7,8  and  9  we  would  obtain  results  for  the 
third  row;  however,  only  the  first  two  differences  of .this  row  are  defined. 

Similarly,  if  J  took  on  values  12,13,14  and  I5  we  would  obtain  results 
for  the  4th  row;  here  the  first  three  differences  are  defined. 

Finally,  if  J  took  on  values  18,19,20  and  21  we  would  obtain  the  differ¬ 
ences  for  the  514  row.  In  row  5  and  subsequent  rows  all  differences  are  defined! 

For  convenience,  we  will  permit  the  computer  to  compute  and  record  the 
incorrect  differences  in  the  2nd,  3rd  and  4th  rows  to  maintain  generality. 

Prior  to  printing  the  desired  table,  we  will  instruct  the  computer  to  replace 
the  few  "incorrect"  differences  with  their  correct  zero  values, 

To  achieve  the  generality,  note  that  the  index  j  is  set  to  zero,  and  for 
j  =  0,1,2  and  3  the  general  definition 

E2,3,J  =  E2,2,J  -  E1,2,J 

will  provide  the  four  entries  for  the  second  row.  Prior  to  entering  the 
general  definition,  we  set  another  index,  to  have  value  J  +  4  so  that  we 
can  use  a  COUNT  statement  to  automatically  advance  j  by  unity  and  terminate 
the  computations  after  the  four  differences  for  a  given  row  have  been  computed 
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and  recorded.  The  COUJNT  statement  -will  have  the  form 


count(jo)in(j)goto(nxtlin) 

After  the  differences  for  a  given  row  have  been  computed  and  recorded, 
we  instruct  the  computer  to  determine  if  all  rows  have  been  completed.  Again 

we  can  use  a  COUNT  statement  for  this  purpose  since  we  know  the  precise  number 

of  rows  that  we  want  processed.  Further,  if  all  rows  have  not  been  processed, 
the  existing  value  of  J  should  be  advanced  by  2  to  commence  the  computations 

for  the  next  row.  Again,  the  COUNT  statement  is  convenient  for  both  purposes; 

the  determination  as  to  whether  all  rows  have  been  processed;  and  the  auto¬ 


matic  increase  of  the  index  j  by  2. 
both  is  written  in  the  form 


The  COUNT  statement  which  will  achieve 


CQLTIT  ( ll4 /2)IN(  J  )G0T0(  iRCTLIN 

This  variation  of  the  COUNT  statement  specifies  that  the  current  value  of 
the  index  J  is  to  be  advanced  by  2,  then  determine  if  this  new  value  of  J  <  ll4 
If  J  <ll4,  go  to  next  line,  otherwise  continue  with  the  statement  immediately 
following  this  COUNT  statement.  Note  that  the  specification  of  the  amount  by 
which  the  index  is  to  be  advanced  is  specified  after  the  slash  following  the 
maximum  count,  114- .  (The  student  should  verify  that  the  maximum  count,  ll4, 
will  indeed  provide  for  the  processing  of  the  20th  row  and  terminate  the  compu¬ 
tation  of  the  general  differences.  He  should  also  verify  that  maximum  counts 
Of  109,  110^  111^  112  or  115  would  achieve  the  same  result). 

Note  that  after  all  rows  have  been  processed  we  instruct  the  computer  to 
replace  the  few  incorrect  undefined  differences  with  their  correct  zero  values 
and  then  print  the  desired  table. 


We  stated  previously  that  at  NXTLIN  we  would  instruct  the  computer  to 
let  another  index,  j^,take  on  a  value  equal  to  j  +  4.  Here  we  use  the  con¬ 
venient  SETEA  statement  which  has  the  general  form 

SETEA(A  =  B,C  +  n) 
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whex'e : 

A  and  C  are  index  names: 

B  is  any  name; 
n  is  an  explicit  integer. 

The  resTilt  of  the  general  SETEA  statement  is  that 

Index  A  takes  on  a  value  eQ.ual  to  the  "name"  B  +  the  existing  "value" 
of  C  n. 

(Note^  not  the  "value"  of  B  and  not  the  "name"  C). 

For  our  example  we  wrote 

SETEA(J0  =  J  +  U) 

A  more  general  means  for  operating  on  integers  is  provided  by  prefixing 
formulas  with  the  name  INT  .  Prefixing  formulas  with  lOT  indicates  a  departure 
from  standard  floating-point  operations  and  in  particular  indicates  that  in¬ 
teger  arithmetic  is  to  be  used  to  obtain  the  result.  For  example, 


INT(A  =  B  +  5  -  C  *  D/E) 

is  valid  provided  that  B,C,D  and  E  are  names  of  existing  integers  and  inter¬ 
mediate  operations  do  not  yield  results  which  exceed  the  restricted  range  of 

permissible  integer  representation.  (For  compatibility  on  both  ORDVAC  and 

19 

BRLESC,  index  integer  values  should  not  exceed  2  in  magnitude,  non-index 
integers  are  permitted  and  for  compatibility  on  both  ORDVAC  and  BRLESC  should 
not  exceed  in  magnitude).  A  further  note  of  caution  to  the  novice,  inter¬ 
nal  representation  of  integers  necessitates  special  rules  for  integer  operat¬ 
ions  which  do  not  yield  exact  integer  results.  If,  for  example,  in  the  above 
INT  statement,  existing  values  of  D  and  E  are  1  and  3  respectively,  then 
D/E  =  l/5  and  the  computers  will  yield  1  as  the  result  for  D/E.  In  general, 
if  the  result  of  integer  division  is  a  positive  mixed  number,  (i.e.,  an  in¬ 
teger  part  and  a  fractional  part),  the  integer  result  assigned  is  one  greater 
than  the  integer  part  of  the  mixed  number.  If  the  result  of  integer  division 
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is  a  negative  mixed  number,  the  result  assigned  is  the  signed  integer  part 
of  the  mixed  nimiber. 

The  program,  input  and  output  for  EXAMPLE  11  are  listed  in  FIGURE  26. 
Note  on  line  8  of  the  program  the  statement 

Tn/*^  ).  f—  r~  ^7  ^  /i/ 

HiSjH  =  ~  =  u  % 

This  is  an  exam.ple  of  assigning  a  rramerical  result  to  many  distinct  q\ 

A  result  may  be  assigned  to  as  many  as  15  distinct  quantities  where  the  name 
of  each  is  separated  by  the  equality  symbol.  We  covild  not  conveniently  use 
the  CLFIAR  statement  in  this  case  since  all  of  the  elements  were  not  uniformly 
spaced  in  internal  storage . 
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PKOB 

C9C)6  M.J. 

ROMANELLI  45107  EXAMPLE 

11 

• 

X,I 

F,  I 

IDEL,  I 

2DEL, I 

30EL,  I 

4DEL,  I 

0000001 

93000000 

59783000 

00000000 

00000000 

00000000 

00000000 

0000002 

93000000 

58168000 

-16150000-01 

00000000 

00000000 

00000000 

0000003 

98000000 

55702000 

-24660000-01 

-85100000-02 

00000000 

00000000 

0000004 

lOOOOOOO 

1  54030000 

-16720000-01 

79400000-02 

16450000-01 

00000000 

0000005 

12000000 

1  36236000 

-17794000 

-16122000 

16916000 

•18561000 

0000006 

12300000 

1  31532000 

-47040000-01 

13090000 

29212000 

46128000 

0000007 

13000000 

1  26750000 

-47820000-01 

-78000000-03- 

13168000 

■42380000 

0000008 

13500000 

1  21901000 

-48490000-01 

-67000000-03 

1  1000000-03 

13179000 

0000009 

14000000 

1  16997000 

-49040000-01 

-55000000-03 

12000000-03 

10000000-04 

0000010 

15000000 

1  70740000- 

01-99230000-01 

-50190000-01- 

49640000-01- 

•49760000-01 

0000011 

16000000 

1-29200000- 

01-99940000-01 

-71000000-03 

49480000-01 

99120000-01 

0000012 

16500000 

1-79120000- 

01-49920000-01 

50020000-01 

50730000-01 

12500000-02 

0000013 

17000000 

1-12884000 

-49720000-01 

20000000-03- 

49820000-01- 

■10055000 

0000014 

17500000 

1-17825000 

-49410000-01 

31000000-03 

11000000-03 

49930000-01 

0000015 

18000000 

1-22720000 

-48950000-01 

46000000-03 

15000000-03 

40000000-04 

0000016 

18500000 

1-27559000 

-48390000-01 

56000000-03 

10000000-03- 

•50000000-04 

0000017 

19000000 

1-32329000 

-47700000-01 

69000000-03 

13000000-03 

30000000-04 

0000018 

19300000 

1-35153000 

-28240000-01 

19460000-01 

18770000-01 

18640000-01 

0000019 

19800000 

1-39788000 

-46350000-01- 

-18110000-01- 

37570000-01- 

56340000-01 

0000020 

20000000 

1-41615000 

-18270000-01 

28080000-01 

46190000-01 

83760000-01 

0000021 

FIGURE  26. 


EXAMPLE  12. 


Specific  indications  of  errors  detected  diiring  compilation  of  a  program 
are  recorded  on  an  output  device  to  inform  the  programmer  of  the  particular 
errors  detected.  To  illustrate  the  particular  information  obtained  pertain¬ 
ing  to  detected  errors,  we  wrote  a  short  program  which  contains  several  pro¬ 
gramming  errors.  We  subiuitted  the  program  to  the  computer^  obtained  and 
listed  the  results  so  that  one  can  correlate  the  information  produced  with 
the  particular  statements  of  the  program  submitted. 


The  program  and  output  obtained  are  listed  in  FIGURE  27  below. 
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PROB  C9()6  M.J.  ROMANELLI  45107  EXAMPLE  12 
START  X=Y+2it  U=V-W*  A=B*C  0=S1N(X/W) 

E=X, I+G1,2,K 
F=785,.3<9'99 
I F-A8S(X=Y) GOTO (CHINA) 

IF-INC(J=0) GOTO (ALASKA) 

F(X)=X**2 

C0UNT(50) IN( J)GOTO{FRANC£) 
MOVE(N)NOS.FROM(A, I ) TO ( B , J ) GOTO ( S I  AM ) 

END  GOTO(BEGIN) 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


ERROR 

45 

ILL..  ~ 

3 

2 

=SIN{X/W 

X=Y+Z*  U=V 

PROB 

C906 

M.J. 

ROMANELL 

ERROR 

03 

2  COMMAS 

2 

3 

,K 

E=X, I+Gl,2 

PROB 

C906 

M.J. 

ROMANELL 

ERROR 

35 

ILL,<  OR  > 

1 

4 

99 

F=785.3<99 

PROB 

C906 

M.J. 

ROMANELL 

ERROR 

08 

ABS  IN  IF= 

1 

5 

Y 

IF-ABS(X*Y 

PROB 

C906 

M.J. 

ROMANELL 

ERROR 

11 

ILL  PAR  IF 

1 

6 

(  J=0 

IF-INC( J=0 

PROB 

C906 

M.J. 

ROMANELL 

ERROR 

01 

ILL,  O.T. 

1 

7 

(X)=X**2 

F( X)=X**2 

PROB 

C906 

M.J. 

ROMANELL 

ERROR 

14 

CL  MV  NO  % 

4 

9 

SIAM) 

MOVE(N)NOS 

PROB 

C906 

M.J. 

ROMANELL 

ERROR 

23 

END  SYMB. 

2: 

10 

) 

GOTO! BEGIN 
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A 

105 

MU 

I 

OOS 

lU 

W 

ION 

MU 

INOS. 

OSO 

B 

106 

MU 

J 

OOK 

I 

X 

lOJ 

MU 

?SUBS. 

N7  0 

BEGIN 

107 

MU 

START 

100  L 

u 

Y 

lOF 

MU 

FRANCE 

108 

MU 

U 

lOK 

MU 

Z 

lOL 

MU 

GIL 

109 

MU 

V 

lOS 

MU 

X  INDEX 

OON 

FIGURE  27. 


Programs  in  general  contain  many  statements  and  consequently  m^any  cards. 
Since  errors  are  detected  during  compilation^  of  immediate  concern  to  the 
programmer  is  the  recognition  of  the  errors.  In  particular,  which  card  or 
cards  contain  the  errors?  Identifying  the  cards  which  contain  the  errors 
presents  no  prohlern  to  the  compiler;  however,  it  is  the  responsibility  of 
the  programmer  to  use  the  means  provided  for  this  purpose.  This  requires  the 
programmer  to  uniquely  identify  each  card  of  his  program  so  that  the  compiler 


will  provide  the  unique  identification.  In  particular,  the  programmer  should 
use  the  identification  columns  provided  on  the  standard  programming  sheets 
and  as  suggested  previously,  simple  numerical  ordering  will  suffice  to  quickly 
identify  the  card  or  cards  vrhich  contain  errors.  In  general,  wiien  the  com¬ 
piler  detects  an  error,  it  will  record  a  line  of  information  on  the  output 
device,  (cards  or  teletype).  A  portion  of  the  information  recorded  is  the 
identification  of  the  card  that  contains  the  error.  Hence  if  the  programmer 
uses  the  simple  numerical  Ordering  scheme,  the  compiler  will  record  on  the 
output  device  the  identifying  number  (or  that  number  +  l)  that  appeared  on 
the  card  containing  the  error.  Other  relevant  information  is  recorded  per¬ 
taining  to  the  error;  however,  identification  of  "the"  card  is  most  important. 
The  results  shown  in  FIGURE  27  correspond  to  the  compilation  errors  detected 
on  BRLESC  ,  Similar  but  less  detailed  resiilts  are  obtained  from  ORDVAC.  Tables 
which  classify  the  tj^'pes  of  errors  are  given  on  pages  78  through  82  of  j^lj  . 
The  essential  difference  between  the  ORDVAC  and  BRLESC  results  is  the  amount 
of  detail  information  recorded  pertaining  to  the  error.  ORDVAC  results  in- 


dics-'tG  "bliG  cl3.s s  1  f'lcs.'biQn  "typG  of  error  ^  "the  card  and  field  of  "blie  card  ■thaf 
contains  the  error.  BRLESC  prints  the  word  ERROR,  the  first  and  subsequent 
characters  of  the  field  containing  the  error,  the  characters  recorded  in 
columns  11  through  20  of  the  card  in  addition  to  the  classification,  card  and 
field  identification.  BRLESC  also  prints  the  first  50  characters  that  were 
recorded  on  the  PROB  card. 


102 


V?hen  an  error  on  a  card  is  detected,  the  remainder  of  that  card  is 
ignored.  Further,  if  the  error  was  detected  near  the  end  of  a  card  the  next 
card  may  also  be  ignored,  hence,  all  errors  may  not  be  detected  in  a  single 
compilation.  The  programmer  generally  uses  a  conventional  process  of 
elimination  of  errors;  i.e.,  he  submits  a  program,  corrects  the  detected 
errors,  re-submits  his  program,  corrects  and  re-submits  until  no  further 
compilation  errors  are  detected.  Depending  on  the  nature  and  complexity  of 
the  program,  several  compilations  may  be  necessary  to  detect  all  of  the  errors. 

Note  in  FIGURE  indication  of  errors  on  cards  2, d,U,5,6,7j9  ^nd  10. 

The  student  should  study,  identify  and  correct  each  of  the  errors  indicated. 
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EXAMPLE  15. 


The  fact  that  the  compiler  does  not  detect  any  errors  in  a  given  program 
does  not  infer  that  there  are  not  errors  in  the  program.  Another  type  of  error 

TatK Tr*V»  mo'ir  r»Q>n  r\r»r»nv»  o vo r*i ■?  r\n  r\'P  +  Vi ~I  ■rN-r»/-\ rr-r»cj m 

ri  j.  j  \A.  v^v_<.i.x  '-4.  MX  N_  S_/X  L/XXS_r  ^  WllX^X.-!.^  VU.  dill  • 

For  example^  one  may  write  a  valid  statement 


7  _  Gm5fT^/V^ 

XJ  — 


and  through  an  error  in  input  or  otherwise  the  existing  value  of  X  may  become 
less  than  zero.  Similarly  one  may  write 

z  ^  arcstnCx') 

-  - - -  \  ""/ 

and  at  some  point  in  the  computation  X  may  attain  an  existing  value  much  greater 
than  one  in  ma<?nitude  due  to  an  error  in  innut  or  otherwise . 


Errors  which  occur  during  execution  are  not  as  easily  identified  as  those 
detected  during  compilation.  Here  again,  the  responsibility  of  correlating 
the  error  with  the  program  lies  with  the  programmer.  Many  of  the  function  sub 
routines  (the  functions  of  single  and  multiple  arguments)  include  checks  for 
errors.  If  errors  are  detected  during  execution,  information  is  recorded  on 
an  output  device  and  the  computer  will  stop  unless  the  programmer  provides  for 
continuation.  If  for  example,  the  programmer  wrote  twenty  distinct  statements 
involving  the  SQRT  function,  relevant  information  is  recorded  on  the  output 


u.GvJ.'Jc:  pcrx  uclx.111.11^  uu  uiit::  t^xxv-'x  Xii  one;  »-^'qy.\x  ctx  ^oujidi  iiv.-'We;  v  cx  ^  xuidi  ox  i  x  L.ci  o  X»^ii 

as  to  which  of  the  twenty  caused  the  error  is  the  responsibility  of  the  pro¬ 
grammer  . 


To  provide  for  continuation  after  an  error  is  detected  and  relevant  in¬ 
formation  printed,  the  programmer  must  Include  in  his  program  the  location 
name,  EEIROR.  At  the  place  called  ERROR  in  his  program,  the  programmer  may  pro¬ 
vide  for  whatever  action  he  desires.  He  may  for  example,  at  EtffiOR,  instruct 
the  computer  to  print  the  twenty  arguments  of  the  twenty  square  roots  so  that 


he  can  then  ascertain  which  of  the  twenty  caused  the  error.  At  EREOR  one 
could  write  instructions  which  test  each  of  the  twenty  arg’nments  to 
and  identify  which  of  the  twenty  was  less  than  zero  and  provide  for  corrective 
action  and  continuation  without  stopping  the  execution  of  the  program. 


In  the  example  which  follows  we  will  illustrate  the  use  of  ERROR  so  that 
the  computer  will  not  stop  after  the  first  error  encountered.  In  particular. 


at  ERROR,  we  will  instruct  the  computer  to  print  GOOF  AGAIN  and  conti 
the  next  statement  in  the  program.  In  this  way  we  can  illustrate  the 
information  recorded  corresponding  to  several  types  of  errors. 


vn  1  o  T.rn  +  Vi 

ri  ^ 

pertinent 


The  program  and  output  are  listed  in  FIGURE  28. 
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PROB  C906  M. J.ROMANELLI  46107  EXAMPLE  13  1 

START  X=-.000000000i:«  Y=  1000000000 i  P I  =  3 . 1 4 1 5926536  2 

SET(SPD=2.)*  Z=SQRT(X)  3 

2.  StT(SPD=3.)^  W=LOG(X)  4 

3.  SET(SPD=4.)X  Z=Y*Y*  A=EXP(Z)  5 

4.  SET(SP0  =  5.)*  B  =  ARCSIN( 1-1000*X )  6 

5.  SIET(  SIPU  =  6.  )«  C  =  C0S(Y*»2)  7 

6.  SET(SPD=7.)i  D=.3*Y?  ENTER ( POWER ) X ) D )  8 

7.  SET(SP0=8.)«  E=rAN(PI/2)?  PRINK  E  >E  9 

8.  GOTOIN.PROB)  10 

ERROR  PRINT<G00F  AGAIN  >*  GOTO(,SPO)  11 

END  GOTO(STAKT)  12 
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« 


RUN  ERROR  LOG  X  NEG.  MAY. 23,63  C906 
GOOF  AGAIN 

RUN  ERROR  EXP  BIG  X  MAY, 23, 63  C906 


GOOF  AGAIN 

RUN  ERROR  ARCSIN  1+  MAY. 23, 63  C906 
GOOF  AGAIN 

RUN  ERROR  SINCOS  N  S  MAY. 23, 63  C906 
GOOF  AGAIN 

RUN  ERROR  LOG  X  NEG.  MAY. 23, 63  C906 
GOOF  AGAIN 


E  -19595930  12 

RUN  ERROR  NEG.  SORT  MAY. 23, 63  C906 


M. J.ROMANELLI 
M. J.ROMANELLI 
M. J. ROMANELLI 
M.  J.ROMANELLI 
M,J. ROMANELLI 

M.J. ROMANELLI 


45107  EX 
45107  EX 
45107  EX 
45107  EX 
45107  EX  < 


262-10000000000-09 

0000001 
265  14426950409  19 

0000002 
269  10000001000  1 

0000003 
272  15915494309  18 

0000004 

2905-10000000000-09 

0000005 

0000006 


45107  EX 


FIGURE  28. 
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In  EXAMPLE  11,  the  program  vas  designed  to  compute  and  print  a  table  of 

differences  corresponding  to  the  given  and  f(X^)  where  i  =  1,2,..., 20; 

,  !! 


1  .e . 


it  was  designed  for  a  table  with  precisely  twenty  "rows' .  Itirther, 


internal  storage  space  was  reserved  for  the  complete  table  of  differences. 

In  this  example,  we  will  design  another  program  for  a  solution  to  the  problem 
of  EXAMPLE  11,  to  introduce  two  new  concepts:  the  program  will  be  designed  for 
a  table  of  'Variable "  length,  i  <  1000;  internal  storage  space  will  be  reserved 
for  the  given  data  and  only  two  rov7s  of  differences . 

To  provide  for  reading  and  recording  a  "variable"  amount  of  data,  the 
READ  statement  specifies  the  maximum  number  of  quantities  to  be  read  and  re¬ 
corded  and  the  variable  amount  of  data  is  terminated  with  a  blank  card.  The 
computers  consider  a  READ  statement  fulfilled  whenever  the  specified  number  of 
quantities  have  been  read  and  recorded  or  whenever  a  blank  card  is  encoimtered. 
(if  a  blank  card  is  encountered  as  the  first  card,  it  is  ignored).  Correspond¬ 
ing  to  any  READ  statement,  the  precise  number  of  quantities  and  the  precise 
number  of  cards  that  were  read  and  recorded  is  made  available  to  the  programmer. 
In  particular,  after  a  READ  statement  has  been  executed,  the  number  of  quanti¬ 
ties  that  were  read  and  recorded  is  available  in  integer  form  in  index  09 , 
Similarly,  the  number  of  cards  that  were  read  and  recorded  is  available  in 
integer  form  in  index  08.  Each  subsequent  READ  statem.ent  erases  the  previous 
integers  in  08  and  09  and  consequently  the  existing  integers  in  08  and  09  at 
any  given  time  correspond  to  the  latest  READ  statement  that  was  executed. 


(Blank 


not  included  in  the  integer  counts ) . 


To  instruct  the  computer  to  read  and  record  the  given  X^.  ,  i  =  1,2,  ...,NX 
where  NX  <  1000,  'we  write 

read(iooo)nos.at(xi)^  NX  =  09 

The  above  assumes  that  the  X^.  are  in  standard  form,  six  per  card.  If  NX  <996, 

then  tVip  X  are  terminated  with  a  blank  card.  If  for  example,  there  were  997  X*s, 
1  "  '  ■ 
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X  Q  would  be  recorded  in  the  first  field  of  the  167th  card^  and  even  though 

I 

the  remainder  of  this  card  were  blank  the  computer  would  substitute  zeros  for 
the  values  of  ^999  ^1000  ‘  integers  in  09  and  08  would  be 

respectively  iOOO  and  IbT-  The  statement 


ArV  _  AO 
inA  —  wy 


tells  the  computer  to  assign  a  "value"  to  NX  equal  to  the  existing  integer  in 
09"  Since  "09"  is  an  absolute  machine  name,  the  integer  9  is  not  the  value 
assigned  to  NX,  rather  the  "value"  in  index  09  is  assigned  as  the  existing 
value  of  NX.  We  wrote  the  NX  =  09  statement  immediately  following  the  READ 


statement  so  that  we  will  have  in  NX  the  number  of  X's  read.  After  the  f(X. ) 

1 

are  read  and  recorded,  we  can  instruct  the  computer  to  check  if  the  number  of 

X's  and  f's  read  are  identical.  To  read  and  record  the  f (X. )  we  write 

1 


READ(1000)N0S.AT(F1)  i  NF  =  09 


To  include  the  program  check  we  write 

IF-INT(NX  =  nf)goto(work)^ 

This  conditional  statement  is  immediately  followed  by  a  PRINT  statement  which 
identifies  and  prints  the  number  of  X's  and  the  number  of  f's  in  the  event 
there  is  a  discrepancy. 

To  determine  the  internal  storage  space  required  for  the  four  differences 
of  a  given  row,  we  illustrate  schematically  below  the  quantities  required  for 
the  computation  and  printing.  We  plan  to  compute  a  row  of  differences  and  print 
that  row  immediately. 


4-1 

4-1 

A^ 

1-1 

1.  _  __  j 

X. 

1 

f(X.) 

.3 

AT 

1 

p 

AT 

1 

I> 

U 

A. 

1 

D1  D2  D3  D4 
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If  we  let  D1,D2,D3  and  Bk  be  the  names  of  the  four  differences  in  the  ith 
row  and  P1,P2,P5  the  names  of  the  required  differences  in  the  (i-l)st  row, 
we  can  write 

D1  =  F,I  -  F,(l-l) 

D2  =  D1  -  PI 
D3  =  D2  -  P2 
D4  =  D5  -  P3 

where  we  prepare  in  advance  for  the  computations  in  the  next  row  by  writing 

PI  =  D1 
P2  =  D2 
P3  =  D3 

so  that  these  quantities  are  available  for  the  computations  of  the  next  row  of 
differences.  (We  will  make  special  provisions  for  printing  rows  with  \in- 
defined  differences)*  We  have  used  the  above  notation  in  the  flow-chart 
Illustrated  in  FIGURE  29. 
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23  JULYt63  PAGE  1 

PROS  C90o  M. J.ROMANELLI  45107  EXAMPLE  14.  1 

BLOCC X-XIOOO ) F-F 1000 ) D1-D4  2 

START  KEAO(  iOOO)NCS.AT(Xl)^  NX  =  09%  READ(  1000 )  NOS .  AT  {  Fi  )  *  iNF  =  09  3 

IF-I NT {NX=NF )  GOTO(NORK)  4 

PK1NT-IF0RMAT(  SEEl.  )-<nUM8ER  OF  X,S  =  >NX  5 

PRINT-FORMAT ( SEEl. ) -<NUMBER  OF  F,S  =  >NF^  GOTOIN.PROEn  6 

WORK  PRINT<  I  X,I  F,1  IDELI  2DELI  3DELI>  7 

C.ONT<  40ELI>«  ENTER!  PRINTS)*  SET!  1  =  1)*  8 

PRliVr-FORMAT  !  S£ E  1 .  )  -  !  1)  X  ,  1  )  F  ,  1  9 

AATl  C0UNT!NX  +  1)  INI  I 1GOTO!N.LINE)*  GOTOIN.PROB)  10 

N.LINE  D1  =  F  ,  I~F  ,  !  I-l  )  *  02  =  li1-PU  D3=D2~P2*  D4  =  D3-P3  11 

P1  =  U1*  P2  =  l)2*  P3=D3%  IF-1NT!I>4)  G0T0!GENPR)  12 

PRINT-FORMAT ! SEEl. )-! I )X, 1 ) F, 1 )( I-l  )N0S.AT!D1)*  G0TO!AATI)  13 

GENPR  print-format ! SEEl. )-! I )Xt 1 )F, I )  ! 4  I  NOS. A T ! D1 ) *  G0T0!AAT1)  14 

SEEl.  F0RM14-5 ) 3-2 ) 1-1 ) 12-1-7) 3-2) 1-6)2  15 

END  GOTO! START)  16 

93000000  00  95000000  00  98000000  00  lOOOOOOO  01  12000000  01  12500000  01 

13000000  (Jl  1  3500000  01  14000000  01  15000000  01  16000000  01  16500000  01 

17000000  01  17500000  01  18000000  01  18500000  01  19000000  01  19300000  01 

19800000  01  20000UOO  01 

59783000  00  58168000  00  5!)702OU0  00  54030000  00  36236000  00  31532000  00 

26750000  00  21901000  00  16997000  00  70740000-0 1-29200000-01-79 120000-01 
-12384000  00-17825G00  00-22720000  00-27559000  00-32329000  00-35153000  00 
-39788000  00-‘ti615000  00 
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I 

X,  I 

F,I 

IDELI 

2DELI 

3DELI 

40EL  I 

1 

.9300 

.5978 

2 

.9500 

.5817 

-- 

.0  161 

3 

.9800 

.5570 

- 

.0247 

-  .0085 

4 

1.0000 

.540  3 

-- 

.0167 

.0079 

.0165 

5 

1.2000 

.  3624 

- 

.  1779 

-  .1612 

-  .1692 

-  .1856 

6 

1.2500 

.3153 

-- 

.0470 

.  1309 

.2921 

.4613 

7 

1.3000 

.2675 

- 

.0478 

-  .0008 

-  .1317 

-  .4238 

8 

1. 3500 

.2190 

~ 

.0485 

-  .0007 

.0001 

.1318 

9 

1.4000 

.1700 

-- 

.0490 

-  .0005 

.0001 

.0000 

10 

1.5000 

.0  707 

- 

.0992 

-  .0502 

-  .0496 

--  .0498 

11 

1.6000 

- 

.0292 

“ 

.0999 

-  .0007 

.0495 

.0991 

12 

1 .6500 

- 

.0791 

- 

.0499 

.0500 

.0507 

.0013 

13 

1 .7  000 

- 

.1288 

- 

.0497 

.0002 

-  .0498 

-  .1006 

14 

1.7500 

- 

.1783 

~ 

.0494 

.0003 

.0001 

.0  499 

15 

1 . 8000 

- 

.2272 

- 

.0489 

.000  5 

.0002 

.0000 

16 

1.8500 

- 

.2  756 

- 

.  0484 

.000  6 

.0001 

--  .0001 

17 

1 . 9000 

- 

.3233 

- 

.0477 

.  0007 

.0001 

.0000 

18 

1.9300 

- 

.3515 

- 

.0282 

.0195 

.0188 

.0186 

19 

1.9800 

- 

.3979 

- 

.0464 

-  .0181 

-  .0376 

-  .0563 

20 

2.0000 

- 

.41u2 

-- 

.0183 

.0281 

.0462 

.0838 

21 

.  0000 

.0000 

.4161 

.4344 

.406  3 

.3602 

22 

.  0000 

.  0000 

.0000 

-  .4162 

-  .8506 

-1.2569 

23 

.  0000 

.0000 

.0000 

.0000 

.4161 

1.2667 

24 

.  0000 

.0000 

.0000 

.0000 

.0000 

-  .4162 

FIGURE  30. 


0000001 


correspond  to  I  =  21,22,23  and  24.  To  eliminate  these,  we  could  have  used 
a  second  terminating  condition,  namely: 


IF  (XI, I  =  FI, I  =  0)G0T0(N,PR0B) 


This  should  he  inserted  after  the  PRINT  statement  on  line  9  oi’  on  a  card 


j.uj.Lin:;\A  J.CX  jr  _pj.  ^  vj.  j.  vu.  u-\-/ • 


il4 


HiXAMrliE  15 . 


This  example  illustrates  several  of  the  convenient  ENTER  statements. 

They  provide  a  means  of  evaluating  functions  which  may  req^uire  more  than  one 
argument  and  for  which  more  than  one  result  may  be  produced.  As  in  the  class 
of  single  valued  fimctions  of  one  argument,  most  of  these  functions  require 
floating-point  arguments  and  yield  floating-point  results.  Exceptions  will 
be  explicitly  stated  in  the  descriptions  which  follow.  Arguments  for  the  tri¬ 
gonometric  functions  and  results  for  the  inverse  trigonometric  functions  are 


GIVEN: 


A  = 
B  = 
C  = 
D  - 


-  2 

3 

A. 5 

-  6.5 


REQUIRED; 

Compute ;  ; 

print  : 

SA  = 

sin(A) 

WC 

CA  = 

cos(A) 

FC 

BC  = 

B^ 

WI 

©'  = 

arcsin(SA) 

FI 

CliiU. 

arccos(CA) 

m. 

W  =  i 

arctan(B/C) 

CI 

WC  =  whole  part  of  C 
FC  =  fractional  part  of  C 


integer  form  of  C 
integer  form  of  FC 
integei 


floating-pt  representations 


integer  representations 


-  floating-point  of  M  V  floating-pt  representa.tion 


Tne  flow  chart  for  this  example  is  shown  in  FIGURE  31- 
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FLOW 


Print  and  Identify  : 


CA;  TriETA;  BCj.W;  WC;  FC;  ul 
FI;  NI 


_ ! _ 

Print  one  blank  card. 

Set  output  card  counter  =  0 


Print  (standard  form) 

CA;  BC;  THETA;  W;  WC;  FC;  Cl 


ix  FOR  i^.XAI'lPLE  15 
TGURE  31. 


16 


Itie  statement 


ENTER  ( SINGQS  )  A  )  SA  )  CA  ) 

tells  the  computer  to  compute  both  the  sine  and  cosine  of  A  and  to  assign 
the  results  to  SA  and  CA  respectively. 

The  statement 

enter(power)b)c)bc 

tells  the  computer  to  raise  the  quantity  B  to  the  C  power  and  to  assign  the 
result  to  BC.  Since  the  logarithm  function  is  used  to  obtain  this  result, 

B  >  0. 


The  statement 


enter  (  ARCSC  )  SA  )  CA  j'i'tlETA 


tells  the  computer  to  determine  the  angle  THETA  whose  arcsin  and  arccos  are 
respectively  SA  and  CA. 

(-IT  <  THETA  <  3t)  . 


The  statement 


enter(artan)b)c)w 

tells  the  computer  to  determine  the  angle  W  whose  arctan  is  B/C.  (-:t<  W<  fl). 


The  statement 

ENTER ( WH . FRA ) C )WC )FC 


tells  the  computer  to  separate  the  floating-point  quantity  C  into  its  whole  and 
fractional  parts,  assigning  the  whole  part  to  WC  and  the  fractional  part  to  FC. 
Both  parts  are  recorded  in  floating-point  form. 


The  statement 


enter(cvftoi)c)wi 


tells  the  computer  to  convert  the  floating-point  quantity  C  to  integer  form, 
assigning  the  result  to  WI .  (For  rounding  purposes,  the  quantity  10  is  first 
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rs 

CXU.U.CU.  UU 


n 


ucx  ux  c: 


l:uii vex  JDXUii  xa 


ex  X  ee  ueu 


^ 

j  • 


The  statement 


ENTER ( CVITOF )NI )CI 


tells  the  computer  to  convert  the  integer  quantity  Nl  to  floating-point  form, 
assigning  the  result  to  Cl. 

The  statement 

ENTER (PRINT  B) 

tells  the  computer  to  print  one  blank  card. 

The  statement 

ENTER  (ZEROCC) 


+.<=k11c  +."h<=k  pniTiTMl+.^r'  pniin  +  o-r'  omi«1  +  '7or*<^  IJV^on  ■?  n  + -?  n  rr 

SMi  k./'ww  wi.i'w  s_  Vw  w  A  waw  J.  w  -4-  ^s^Xa-/*  rfxxa_.xx  L/^xXf^ 

standard  card  formats,  the  successive  cards  produced  contain  serial  identifi¬ 
cation  in  columns  77  through  80,  (output  card  count).  The  serial  count 
commences  at  one  with  the  first  output  card  produced  and  is  advanced  by  one 
with  each  successive  card  produced  until  either  the  count  reaches  9999  or  when 
restored  to  zero  by  the  above  ENTER  statement. 


Non-standard  formats  may  also  include  the  output  card  count,  however,  it 
must  be  specified  in  the  format.  We  used  a  special  format  for  the  integer 
quantities  that  are  to  be  printed.  The  format 

K  FORM  (U-U)U-U)U-U)2 


is  of  the  form 


^ m  T  T  T  ^ 

\  X— uyx— xj yx— xjy 


where  in  each  case  the  length  is  U  and  the  type  is  U.  The  U  type  tells  the 
computer  that  the  quantity  exists  internally  as  an  integer  and  is  to  be  re¬ 
presented  in  integer  form  in  the  four  columns  allocated,  the  units  position 
on  the  card  corresponding  to  the  fourth  column.  The  characters  specified  with 
in  "t-iio  ^  snd.  ^  oi  "t/lio  PRINT  s’t-s.’tcnicn'ts  Svipcr*c0d.s  "tho  column  s-HoGsliiion  spcci 
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fled  in  "the  foxuiaoj  i.e.^ 


<*»  »^  •?  try  rr  /~v  ^  To  p«  ^  Q  ^  ^  wi^  *o  ^ 

l;  i-'x  X  c:  Djj;vjiiu.liig  oO  oii^  o  uo,  c/«=:iin3ii  o 


PRINT -FORMAT(k)  -<WIb=b>WI<bbFI>-  — 

WI  is  first  recorded  in  columns  1  and  2,  then  a  blank  in  column  3;  the 
equality  in  coliunn  4^  a  blank  in  column  5  then  the  integer  in  columns 
6,7,8  and  9  ■with  insignificant  zeroes  suppressed,  (i.e.,  replaced  by  blanks). 
Next,  blanks  are  recorded  in  columns  10  and  11,  FI  recorded  in  12  and'  13,  etc. 

The  program,  input  and  output  for  this  example  are  listed  in  FIGURE  32. 
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START 


K 


PROS 


FORM 

END 


C906  Ml.  J.ROMANELLI  A5107  EXAMPLE  15 
A=-2*  B=3X  C==4„535  D=-6.3 

ENTER  (SINCOS)A)SA)CA)*  ENTER  (  POWER )  Ei )  C  )  BC 
ENTER  (ARCSC]ISA)(:A)TIHETA%  ENTER  (ARTAN)B)C)W 
ENTER  (WH.FRAEOWC  )FC*  ENTER  (  C  VFTO  I  )  C  )  W  I 
ENTER  (CVFTOI  )  FCJF  I  )*  ENTER  (  CVFTO  I  ]l  D )  N I 
ENTER  (CVITOF)NI 1CI« 


PRINT<SA  =  >SA< 

CA  =  >CA< 

THETA 

=  >rHETA 

PRINT<BC  =  >BC< 

w  =  >w  < 

WC 

=  >WC 

PRINT<FC  =  >f-C< 

Cl  =  >CI 

PRINT-FORMATCK)- 

<WI  =  >WI< 

FI  = 

>FI<  NI  =  >NI 

ENTEREPRINT  6)X  ENTE R ( ZEROCC ) * 

PRINT(SA)CA)THE:TA)BC)W)WC)FC)CIX  GOTO  (N.PROB) 

(4-4) 4-4) 4-4) 2 

GOTO  (START) 
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1 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


SA  := 

-90929743 

CA  = 

-41614684 

THETA  =  -20000000  1 

0000001 

BC  = 

14029612  3 

W  - 

58800260 

WC  =  40000000  1 

0000002 

FC  = 

50000000 

Cl  == 

-60000000  1 

[ 

0000003 

Ml  = 

4  FI  - 

Nl 

1  =  -  6 

90929743 

-41614684 

-20000000 

1  14029612  3  58800260 

40000000  1  0000001 

50000000 

-60000000 

1 

0000002 

FIGURE  52. 
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In  many  problems,  the  given  functions  are  not  always  known  analytically 
and  in  lieu  of  explicit  analytic  definitions  discrete  tabular  data  is  given. 
Similarly,  as  illustrated  in  the  previous  examples,  many  solutions  obtained 
from  the  computers  are  expressed  in  tabiilar  form.  Hence,  in  many  problems 
conventional  interpolation  is  often  desired.  This  example  is  designed  to 
illustrate  how  the  ENTER  statement  provides  for  interpolation.  The  statement 


has  the  general  form 


EN’TER(D.B.IK)x)FX)Xl)Fl)tpt)n)ix)if)^ 


1231^5678 

where  C.D.IN  is  an  abbreviation  for  divided  difference  interpolation.  As  many 
as  8  parameters  may  be  specified  to  describe  the  particular  interpolation 
desired.  It  is  assumed  that  a  table  of  discrete  values,  X^  and  f(X^.  ),  is 
recorded  in  internal  storage,  where  i  =  1,2, . . . ,tpt-and  the  X^  are  in  monotone 
increasing  or  monotone  decreasing  sequence.  *  It  is  further  assumed  that  the 
tabular  values  are  uniformly  spaced  in  internal  storage.  (This  does  not  mean 
that  uniformity  refers  to  the  space  allocation  in  internal 

storage).  Further,  the  -uniform  spacing  of  the  f(X^)  may  indeed  be  different 
from  the  uniform  spacing  of  the  independent  -variable,  X^.  In  general,  if  we 
■want  the  computer  to  perform  particular  interpolations,  we  must  specify  the 
pertinent  information.  The  parameters  which  follow  the  ENTER(D.D. IN)  tell 
the  computer  the  specific  information  that  is  required  for  the  interpolation. 
The  meanings  for  the  parameters  are  as  follows; 


X  is  the  name  of  the  argument,  i.e.,  we  -want  the  interpolated  value 
for  f(X); 

FX  is  the  name  of  the  result,  i.e.,  the  result  obtained  is  assigned 
as  the  existing  value  of  FX  and  hence  fu-ture  references  to  this 


*  Since  divided  differences  are  used. 


X.  ^  X. ,  for  all  i. 
1  '  1+1 
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r*dciil+  "hir  r»d'Por»r»n  +rt  ■|■'ho  namo  TTY* 

^  s_^  ii_>  ‘  ■*  ■  ^  s_^  'xcs_^  ^  A  ^'_4aii.v_,  j.  ^  k  ^ 


XI  is  the  name  of  the  first  entry  in  the  table  of  the  independent 


•\ra  "K*  n  a  V\  1  o  • 

V  (^X  0.(^1^ U-V..  y 


Fl  is  the  name  of  the  first  fionctional  value  in  the  table; 

tpt  is  an  integer  which  specifies  the  number  of  points  in  the  table ^ 
i.e.,  the  maximum  value  of  i: 

n  is  an  integer  which  specifies  the  number  of  points  to  be  used  in 
the  interpolation j  i.e.,  if  n  =  2,  linear  interpolation  is 
desired,  if  n  =  5^  parabolic  Interpolation  is  desired,  etc., 
n  <  18} 

jlx  is  an  integer  which  specifies  the  uniform  spacing  of  the  in 
internal  storage; 

if  is  an  integer  which  specifies  the  uniform  spacing  of  the  f(X^)  in 
internal  storage . 

Note  that  the  dependent  and  independent  tabular  data  need  not  have  the  same 
uniform  spacing.  (The  last  three  parameters,  n,  ix  and  if  may  be  omitted  only 
if  ix  =  if  =  1  and  n  =  5)*  In  general,  the  lower  case  letters  listed  in 
the  string  of  parameters  of  EOT'ER  statements  refer  to  integer  values.  As  such, 
the  integers  may  be  expressed  explicitly  or  may  be  expressed  as  the  existing 
values  of  indices  provided  the  index  name  is  preceded  by  the  mandatory  connna. 

If  for  example,  the  total  number  of  points  in  a  table  is  92,  one  may  literally 
write  92  in  the  parenthesis  corresponding  to  tpt.  As  an  alternative,  if  one 
used  the  name  J  as  an  index  with  existing  value  92,  for  example,  SET(J  =  92), 

^  m  «*a  T  r  T  t  >0  /-r 

maj  wx  X  yfj  xii  viic:  c»  x  c»  j^'-'ii'xx  ^  • 


GIVEN: 


and  f(x.),  (j  =  1,2, ...,10).  y_. ,  i  =  l,2,5,!l,5' 


X,  <  y,.  <x,^  for  1  =  1, 2,5,1+,  y^  >  x, , 

X  X  XW  P  X' 
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REQUIRED: 


For  each  v. ,  determine  f (v. )  using  2,5.^  and 
'■i'  "i'  “  '  ' 

o  4-  ^ 

Jj'w'-LO.  OX'w'Jli* 

Print  and  Identify:  x;  f(x)  2  pt;  f(x)  3  pt; 
f(x)  5  pt.  where  x  corresponds  to  the  y^. 

The  flow-chart  for  this  example  is  shown  in  FIGURE  35 • 


5  point  inter - 
f(x)  4pt;  and 


FLOW  CHART  FOR  EXAMPLE  l6. 


FIGURE  55- 


One  other  general  remark  concerning  the  interpolation  function  is  that, 
extrapolation  will  be  performed  provided  that  the  point  is  not  "too  far" 
outside  the  given  range.  Indeed  "too  far"  means  that  the  point  must  lie 


within  "one  interval"  at  either  end  of  the  table.  Explicitly  this 


extrapolation  for  the  point  X  will  be  performed  provided 


(X 


tpt 


-  X 


+.r%+. 
‘'Jr'  ■' 


.n) 


or  at  the  other  end, 

X^  -  (X^  -  X^)  <  X  <  XI 


for  monotone  increasing  X^ .  Similar  conditions  hold  for  extrapolations  for 
monotone  decreasing  X^. 

The  program,  input  and  output  for  this  example  is  listed  in  FIGURE  ’^h. 
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JULY, 63  PAGE  1 


PKOB  C906  ROMANELLI  C.L.  45107  EXAMPLE  16  I 

BLOC  ( YI-Y5) IF2-IF5) X1-X20  2 

SYN  (F1=X2)  3 

BEGIN  Yl=.973!«  Y2  =  1.000ri  Y3=1.085*  Y4=U125?  Y5  =  1.32  4 

PRINT<  DIVIDED  DIFFERENCE  I  NT  ERPOIL  A  TI  0N>  5 

ENTERIPRINT  B)  6 

PRINT<  LINEAR  PARABOLIC  CUBIC  QUARTIO  7 

ENTERtPRINT  B)  8 

PRINT<  X  F(X)  2PT  FIX)  3PT  FIX)  4PT  FIX)  5PT>  9 

ENTERIPRINir  8)*  SETII=0)  10 

RSJ  SETIJ=0)  11 

INT6RP  ENTERI D.U. IN) Yl, I ) IF2, J)X1 )F1 ) 10 ) , J+2 )2 )2 )  12 

C0UNTI4) INI J)G0T0I INTERP)  13 

PRINT-FORMATIH)-I Yl,  I  )  14 ) NOS. ATI  I F2 )  14 

COUNT  15) INC  I )GOTOI RSJ)*  GOTOIN.PROB)  15 

H  FORMl 12-1-7)3-2)1-1)12-1-9)3-3)1-4)2  16 

XI  DEC  I .95 ) .81342) .96) .81919) .98 ) .83050)  1 ) . 84 14 7 ) 1 . 02 ) . 8 52  1 1  17 

CONTI  1.05) .86742) 1.06) . 8 72 3b )  1 . 09 ) . 8866 3 ) 1 . 1 1 ) . 89570 ) 1.  13). 90441  18 

ERROR  PRINr<  X  IS  PROBABLY  OUTSIDE  THE  RANGE  >%  GOTOIN.PROB)*  19 


20 


END  GOTOIBEGIN) 


MAV.23»63  BKLESC  FORAST  F62 

PROB  C906  ROMANtLLI  C.L.  45107  EXAMPLE  16  * 

DIVIDED  DIFFERENCE  INTERPOLATION  0000001 

LINEAR  PARABOLIC  CUBIC  QUARTIC  0000002 

X  FIX)  2PT  FIX)  3PT  FIX)  4PT  FIX)  5PT  0000003 

.9730  .826541  .826576  .826578  .826579 

1.0001  .841523  .841524  .841524  .841524 

ro  1.0850  .884252  .884309  .884308  .884309 

1.1250  .902233  .902266  .902266  .902267 

RUN  ERROR  O.D.IN  MAY. 25,63  C906  ROMANELLI  C.L.  45107  EXA  286  13000000000  1 

+UJ04+USUJV34UK3VMV43K+45NV4TLU/  ♦NS VTMT2MUT5MVU5VN0NNV  ED-05 

X  IS  PROBABLY  OUTSIDE  THE  RANGE  0000004 


IT  CURE  ^k. 


A  s't'iid.jy  of  i!/li0  program  f  02r  fhio  o^c^jiiplo  will  illus  Irals  som0  row  coricopls « 
First;  we  could  have  used  the  general  READ  statement  to  record  the  given 
tabular  data  in  internal  storage.  However,  to  illustrate  another  means  of 
recording  data  in  internal  storage  we  have  used  the 


DEC  order -type  word. 


Note  you  will  find  on-  line  IJ ,  XI  in  the  LOCATION  coliimns  and  DEC  in  the  order- 
type  columns.  Beginning  in  column  11  we  have  written  the  numerical  values  that 
o  o  4-  ^  vn  -C  f  vn  ^  V  o  f  v  o  ^  To  v.  ^  q  H  4ro  n 

O  Ox  jTObpVJiiU.  OVJ  ./V_L  y  J.  ^  J  y  ±  J  y  J.  V  (-*.x 


S  c*  c*  ciT^O  v»o  +  cs  "K-rr 
. .  C-V  uj 


L4.V.;  j.  k 


a  parenthesis.  These  numerical  quantities  are  recorded  in  consecutive  spaces 
in  internal  storage.  The  first  n\merical  value,  -95}  is  associated  with  the 


name  XI  recorded  in  the  location  columns.  This  fact  alone  does  not  establish 


X2,  X3,  etc.  as  valid  names.  To  associate  these  quantities  with  particular 
names  for  future  references  we  have  used  BLOC  and  SYN  statements.  In  the  BLOC 
statement  we  wrote  (Xl  -  X20)  to  reserve  space  and  establish  X  names  for  these 
twenty  quantities.  In  the  SYN  statement  we  wrote  (FI  =  X2)  indicating  that 
the  second  numerical  quantity,  .815^2  was  to  be  associated  with  either  name, 

X2  or  FI.  Note,  also,  that  under  this  ordering  and  arrangement,  the  X’s  are 
two  units  apart  in  internal  storage.  Likewise  the  corresponding  functional 
values  are  two  units  apart  in  internal  storage.  These  concepts  permit  us  to 
establish  the  data  and  names  required  for  the  interpolation.  We  will  make 
reference  to  XI,  and  FI  and  specify  that  the  tabular  data  is  two  units  apart 
both  for  the  X  values  and  F  values. 


■Wc  used  names  IF2  -  IF5  for  the  interpolated  functional  values.  We  have 
also  used  a  non-standard  format 

H  F0RM(12-1 -7)5-2)1-1)12-1-9)5-5 )l-^)2 

for  the  entries  in  the  output  table.  The  form  (12-1-T)  specifies  that  the 
nurriGi’ics-l  vs-Iug  of  Yl^I  Gxisfs  in  in'tGx'ncil  s'toi'a.gG  in  flosuing-poin't  fomn^  'tliG 
GxtGrnal  form  is  to  bo  printGd  -with  a  decimal  point  after  one  digit  and  7 
columns  are  to  be  used  for  entire  external  representation.  The  form  (5-2) 
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indicates  that  2  blank  spaces  are  to  precede  the  next  qmntity  printed. 

The  form  (l-l)  specifies  that  the  preceding  forms  are  to  be  used  just  1  time. 
This  was  inserted  since  we  will  indicate  repetitions  of  another  form  and 
repetitions  commence  with  the  form  immediately  following  the  previous  re¬ 
petition  form.  In  the  above  example^  the  repetition  from  (l-l)  applies  to  the 
(12-1-7)5-2  whereas  the  (l-4)  repetition  form  applies  to  the  (12-1-9)5-5) • 

The  form  (12-1-9)  is  similar  to  the  Initial  form  (12-1-7),  the  only  difference 
is  in  the  total  number  of  columns  to  be  allocated  for  the  quantity. 

The  form  (5-5)  Indicates  5  spaces. 

The  form  (l-4)  indicates  that  the  total  form  (12-1-9)5-5)  is  to  be  used  for 
4  consecutive  quantities,  namely  the  four  interpolated  functional  values, 

IF2,  IF5,  IF4  and  IF5. 

The  form  )2  terminates  the  format  definition. 


TTVAWOTT?  1 'T 
JLl/VrU*!!  i-UJ  I  • 


This  example  illustrates  how  the  ENTER  statement  is  used  to  solve  a 


4-  cam 

3>  j  o 


II 


h  *1  C' 


J.11  11 


is  expressed  as 


AX  ^  B 


a  ^  a  a  •  •  •  a 

(x  1 

b  \ 

11  12  13  In 

1 

^21  ^22  ^23  '  "  ^2n 

"^2 

^nl  ^n2  ^n3  ^nn 

X 

n 

I'd 

The  problem  in  general  is  to  determine  the  which  satisfy  the 

above  system  when  the  elements  of  A  and  B  are  given.  To  use  the  convenient 
ENTER  statement  to  solve  such  a  system,  the  given  elements  must  be  recorded 
in  consecutive  spaces  of  internal  storage  by  rows .  That  is,  the  prescribed 
order  is  as  follows: 


^11^^12^^15^  •  •  •  ^^ln^^l^^21^^22^  ’  ’  ’  ^®-2n^^2" 


nn 


b  . 

n 


One  should  reserve  spaces  for  the  n  x  (n+l)  consecutive  given  elements,  the 
A  and  B.  bv  usine  the  BLOC  statement.  Here,  one  mav  define  a  linear  BLOC,  sav 


ni-j.^ 


Al  -  A  |n  X  (n+l)j  or  a  two  dimensional  array,  Al,l  -  A-j^n 
example,  for  a  system  of  four  equations  in  four  unknowns  one  may  choose  the 
linear  bloc.  Al  -  A20 .  or  the  two  dimensional  arrav  Al.l  -  AU.'^.  In  either 


_ _  _ L  T _ J L 1.  -  _ _  -?-l.^ J.T TTlm  I  I  l.»l  )  _ J.  -  J  -."U  — 

ca.be  ^  urie  mua  u  ue  curia xaoeiiu  wiieri  wi*xu±ii^  one  rjnxrin  aua.oeuieriu  wiixcii  x*exex’a 
the  given  elements.  The  ENTER  statement  has  the  general  form 


j.  ^ 

uu 


I  In Ti riT  w  X  / 


\  A  \ 


niiNTiiin )i\)n)y^ )u 


1  ■zn 


where: 


S.N.E.  is  simply  an  abbreviation  for  Solve  Normal  Equations; 

A  is  the  name  of  the  first  given  element^  ^^11^'' 

n  is  an  integer  (or  name  of  an  integer)  which  specifies  the  number  of 
unknowns ; 

C  is  optional  and  if  specified  is  the  name  where  the  result  of  the  first 
unknown  is  to  be  recorded,  i.e.,  generally  the  name  XI  where  the 
solution  is  to  be  recorded.  The  successive  results  are  recorded  in 
consecutive  spaces  in  internal  storage.  The  identical  results  are 
available  elsewhere,  these  are  explained  below. 

D  Is  the  name  of  the  determinant  of  A  This  name  is  optional,  and  if 
specified,  the  value  of  the  determinant  of  A  is  calculated  and  recorded 
at  the  place  called  D. 

In  the  process  of  obtaining  the  solution,  the  computer  destroys  the  given 
A  and  B.  In  particular,  the  original  A  is  replaced  by  the  inverse  of  A,  and 
B  is  replaced  by  the  solution,  X,  which  satisfies  the  given  system.  Hence  one 
need  not  specify  C  since  one  can  make  references  to  the  solution  by  referring 
to  the  names  of  the  original  B*s.  However,  if  one  wants  to  specify  D  and  not 
C  then  the  ENTER  statement  must  be  written  in  the  form 

EI«’ER(S.N.E.  )A)n))D 

To  illustrate,  we  chose  the  simple  system  of  four  equations  in  four  unknowns 
with  solution  XI  =  1,  X2  =  2,  Xk  =  U. 

GIVEN: 


XI 

4 

X2 

4 

X5 

4 

X4 

=  10 

XI 

4 

2X2 

4 

X5 

4 

x4 

=  12 

XI 

4 

X2 

4 

5x5 

4 

x4 

=  16 

XI 

4 

X2 

4 

X5 

4 

4x4 

=  22 
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REQlixRiixi: 


Determine^  identify  and  print  the  X1,X2,X3  and  X4  which  satisfy 
the  above  equations . 


Note  that  the  given  elements  in  the  prescribed  order  by  rows  are; 


1)1)1)1)10)1)2)1)1)12)1)1)3)1)16)1)1)1)4)22 


To  reserve  space  in  internal  storage  for  the  given  elements  and  the 
solution,  we  write  a  BLOC  statement 

BL0C(A1,1  -  A4,5)X1  -  X4) 

To  record  the  given  elements  in  the  reseinred  space  we  use  the  convenient  DEC 
order-typo  word  and  identify  the  first  element  by  writing  Al,l  in  the  LOCATION 
space . 

Al,l  DEC  (1)1)1)1)10)1)2)1)1)12)1)1)3)1)16)1)1)1)4)22 

The  BLOC  statement  above  provides  for  reserving  the  consecutive  space  for  the 
given  elements  and  solution.  The  '*A1,1  DEC”  provides  for  recording  the  given 
system  in  internal  storage  in  the  reserved  space.  The  flow  chart  for  this 
example  is  given  below. 


1 


T\T  ^  ^  V»  o  ^  ^  Q  A  T3  ^  ^  tvi  /^  /^  ^  V»  n  ^  +  Vk  iTn  rr  •?  »  r/^  v^  A  v^  *D  »~«  v»  /-\  v»  /-\  /^  /~» ^  4- 

iiuoc:  oiict  L»  G.  o  uxrvL\x  x  o  xo  ctoouiucu.  0110.0  one  ^xveii  n  oiiu.  u  o±  e  ±eeuiucu.  xiioci  — 

nally  in  the  prescribed,  order.  This  is  indicated  on  the  flow  chart  just  prior 
to  START  with  the  box  containing  A  and  B  appended  with  a  broken  line .  Note 

qIc/^  +V<q+  T.r^  T.TT  T1  •r^v*'{v>+  Qvn/^  +  n  "P-^r  +V*£i  r«/~\ln+'{/^'>n  +T.T"?  cs  '  l^k»  nr*  ■?  r*  +  /-\ 

o-oow  01X0  0  J.11  o  01X0  J.OG IX  oj. X.  one::  ow-i-l4.ox.wii  ow-lwc*  x.ixx.o  X.O  owxxe  ow 

emphasize  that  if  one  specifies  a  place  for  the  solution  (as  we  do  in  our 
prograim)  then  the  solution  res\ilts  are  available  in  two  places.  The  two 
prints  correspond  to  the  solution  references  in  the  two  distinct  places.  The 
program  and  output  for  this  example  are  listed  in  FIGURE  35- 
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23  JULY, 63  PAGE  1 

PROEi  C906  ROMANELLI  C.L.  45107  EXAMPLE  17  1 

BLOC  (A1,1-A4,5)X1-X4  2 

START  ENTER(S.N.E.)Al,l)4)Xi)DA  3 

PRINT-F0RMATIK)-<X1  =  >Xl<  X2  =  >X2<  X3  =  >X3  4 

CONK  X4  =  >X4<  DET  A  =  >DA*  ENTERCPRINT  B)  5 

PRINT-F0RMAT(K)-<X1  =  >A1,5<  X2  =  >A2,5<  X3  =  >A3,5  6 

CONT<  X4  =  >A4,5<  DET  A  =  >DA*  GOTO(N.PROB)  7 

K  FORME  12-2-4) 1-4) 12-3-6)2  8 

Al,l  DEC  (1)1)1)1)10)1)2)1)1)12)1)1)3)1)16)1)1)1)4)22  9 

END  GOTO  (START)  JLO 

MAY. 23,63  BRLESC  FORAST  F62 

PROB  C906  H.J. ROMANELLI  45107  EXAMPLE  17  * 

XI  --=  1.  X2  =  2.  X3  =  3.  X4  =  4.  DET  A  =  6.0 

XI  ==  1.  X2  =  2.  X3  =  3.  X4  =  4.  DET  A  =  6.0 


FIGURE  55. 


EXAMPLE  18. 


In  this  example  ve  vlll  illustrate  another  means  of  obtaining  the  solution 
to  the  problem  of  EXAMPLE  IJ.  Since  the  A  matrix  (the  coefficient  matrix)  is 
symmetric,  less  storage  space  is  required  and  computation  time  is  proportionally 
diminished.  For  lov  order  systems,  the  space  and  time  savings  may  not  be 
appreciable;  however,  for  high  order  systems,  both  space  and  time  savings  m.ay 


be  appreciable. 

Similar  to  the  general  case,  the  symmetric  case  also  requires  that  the 
given  elements  be  recorded  internally  in  consecutive  spaces  by  rows.  In  parti¬ 
cular  the  required  order  is  as  follows: 


a, ,  ,a 
XX' 


^  ^  I  •  • 


^  in’  1'  22’  23 


b  -  .a  .b 
n-X'  nn'  n 


The  symmetry  condition  is  specified  in  both  the  BLOC  statement  and  the  ENTER 
statement  as  follows: 

BL0C(A1,1  -  A4,5/SY.)X1  -  X4 
ENTER (SY . SNE)A1 , 1 )4 )X1 )d 

Again,  as  in  the  general  case,  in  obtaining  the  solution,  the  original  A  matrix 
is  replaced  by  the  inverse  of  A;  however,  the  solution  is  recorded  only  in  the 
specified  place  XI,  and  is  not  recorded  in  the  space  previously  occupied  by  the 
augmented  vector  as  in  the  general  S.N.E.  The  flow  chart  for  this  symmetric 
case  is  similar  to  that  of  the  general  case. 
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Start _ 

ENTER (SY.SNE) 


PRINT: 

XI  X2  =/^X5  =-^XA  w 


FLOW  CHART  FOR  EXAI-IPLE  l8. 

FIGURE  36. 

The  program  and  results  for  this  example  are  listed  in  FIGURE  37- 
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23  JULY, 63  PAGE  I 

PROS  C906  M. J.ROMANELLI  A5107  EXAMPLE  18  1 

BLOC  ( A1 , 1-A4,5/SY. ) XI-X4  2 

START  ENT£R(SY.SNE)A1,1)4)XI)D  3 

PRINT-FORMAT ( Kl-<Xl  =  >Xl<  X2  =  >X2<  X3  =  >X3  4 

CONT<  X4  =  >X4<:  DET  A  =  >0%  GO  TO  { N.  PROB 1  5 

K  FORM! 12-2-41 1-4) 12-3-6)2  6 

Al,l  DEC  ( 1) 1) 1) 1) 10)2) 1) 1) 12 )3)1 ) 16)4)22)  7 

END  GOTO (START)  8 

MAY. 23, 63  BRLESC  FORAST  F62 

PROS  C906  M. J.ROMANELLI  46107  EXAMPLE  18  * 

XI  =  1.  X2  =  2.  X3  =  3.  X4  =  4.  DET  A  =  6.0 


FIGURE  57- 


TTY’AKfPT.Tr  IQ 


This  example  illustrates  the  ENTER  statement  which  applies  Simpson's  rule 
to  approximate  the  definite  integral 


I  f(x)dx  . 
k 


The  statement  has  the  general  form 


ENTER ( S . INTE )f)X )FX )l ) A ) B )E 


■vdiere : 


S.INTE  is  an  ahhreviation  for  Simpson  integration; 

F  is  the  location  name  where  the  integrand  f(x)  is  explicitly 

defined; 

Since  the  process  is  designed  to  provide  the  approximation  for 
arbitrary  f(x)j  it  is  the  responsibility  of  the  programmer  to 
define  the  f(x)  of  interest.  Further,  the  statement  or  state¬ 
ments  which  define  f(x)  must  be  terminated  with  GOTO(S.I.FF) . 
This  enables  the  computer  to  evaluate  and  weight  the  integrand 
as  many  times  as  is  necessary. 


FX  is  the  .name  of  the  integrand,  f(x); 

I  is  the  name  of  the  resulting  approximation; 

A  is  the  name  of  the  lower  limit,  B  the  upper  limit; 

E  is  the  name  of  a  "relative  error  bound; 


The  program  is  designed  to  obtain  the  approximation  with  a  minimum  number 
of  evaluations  of  the  integrand.  Successive  approximations  are  obtained  using 
successively  smaller  AX  until 


or  uii'fcil 


1025. 


When  the  latter  occiirs,  a  "run  error"  print  includes  the  numerical  value  of 


-r  (  AV\ 


i(Ax; 


(For  approximations  which  tend  to  zero,  defining  the  original  integral  as  a 
s\im  of  integrals  may  suffice  to  produce  a  satisfactory  approximation) , 


in(2)  = 


J 


dx 


0 


1  +  X 


3T  =  4  V 

1  r\ 

y  u 


dx 


1  +  X 


REQUIRED: 

Use  Simpson  Method  of  Numerical  Integration  to  obtain  approximations 
for: 


LN(2) 

with 

an  associated 

E  =  .0001  ; 

If 

It  tf 

E  =  .0001  ; 

tf 

It  It 

E  =  .000001  : 

rnii^x 

A'R'rr* 

nnu 

T^Tnr«TTTmr« 

X  X  xr  1  nro  u  xix  o . 

FLOW  CHART  FOR  EXAlvIPLE  19 . 

FIGLIRE  58. 

Note  the  "broken  lines"  to  the  evaluation  of  the  integrands  and  the  "solid"  line 
from  the  evaluation  back  to  the  Simpson  integration.  We  used  the  "broken  lines" 
because  the  progrcunmer  does  not  write  an  explicit  statement  which  says  GOTO(EVINGD) . 
Recall  that  we  specify  this  connection  in  the  ENTER  statement. 
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We  used  the  solid  line  returning  from  the  evaluation  since  the  programmer 
does  write  GC)TO(S.I.FF)  . 

Note  the  solid  lines  from  the  "integration  hoxes"  to  the  "print  boxes" 
which  are  used  to  indicate  that  the  computers  go  to  the  statement  immediately 
following  1116  h:[\i'Ph!W  s'ba'beniGn'fc  stflGi*  tliG  s'ta't-GniGn'b  is  complGlGly  6X6011161 

Hence,  evaluation  of  an  integrand  associated  with  an  ENTER  statement  should 
not  immediately  follow  the  ENTER  statement. 

The  program  and  output  for  this  example  are  listed  in  FIGURE  39* 
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PROS  C906  M. J.ROMANELLI  45107  EXAMPLE  19 
START  A-0%  B=l%  EPS=.0001 

ENTEK(S.INTE)EVH4GO)X)FOFX)LN2)  A)B)EPS 


23  JULY, 63  PAGE 


EVINGD 

SECOND 


THIRD 


PRINr<LN(2)  =  >LN2 


EPS  =  >EPS«  GOTO(SECOND) 


FOFX=l/( 1+X) X  GOTOIS.I.FF) 
ENTER! S. INTE)C)X)F)PI )A)B)EPS 


PI=4*PI-g  PRINT< 


PI  =  >PI< 


EPS  =  >EPS%  GOTO(THIRD) 


F=l/( 1+X*X) %  GOTO(S.I.FF) 


Z=. 000001  *  ENTER! S. I NTE ) D ) Y ) H ) K ) A ) B ) Z 


K=4*K*  PRINT< 


PI  =  >K< 


EPS  =  >1%  GOTO(N.PROB) 


H=l/ !  1  +  Y« Y) i  GOTOIS.I.FF) 
END  GOTO! START) 


MAY. 23, 63  BRLESC  FORAST  F62 

PROB  C906  M. J.ROMANELLI  45107  EXAMPLE  19 


LN!2)  =  69314765 

PI  =  31415925  1 


EPS  =  10000000-03 

EPS  =  10000000-03 


0000001 

0000002 


PI  =  31415927  1 


EPS  =  10000000-05 


0000003 


FIGURE  59- 


EXAMPLE  20. 


This  example  illustrates  how  the  ENTER  statement  is  used  to  obtain  a 
numerical  solution  of  a  system  of  first-order,  ordinary  differential  ecjuations, 


1 

dt 


y[  =  fi(t,yj^(t),  y2(t),. .  .,y^(t)) 


y.  (t  )  =  y. 


1  o 


lO 


i  =  0,1,2, . . U  . 


The  Runge— Kutta— Gill  method  is  applied  to  obtain  an  approximate  solution.  The 
statement  has  the  general  form: 


where : 


ENTER (R.K.G)  At  )n)EVDS)Y)Y‘ )Q 


R.K.G.  is  an  abbreviation  for  Runge-Kutta  Gill; 

At  is  the  name  of  the  incremental  value  of  the  independent  variable; 

The  numerical  solution  is  determined  at  discrete  values  of  the 
independent  variable,  t.  i.e.,  when  initial  conditions  are 
specified  corresponding  to  t  =  t^,  the  solution,  y^(t),  is 
determined  at  t  =  t^  +  At.  The  computer  replaces  the  given 
values,  with  the  new  values,  y^l'tQ  +  At).  These  new 

values  then  serve  as  the  given  values  for  determining  the 
solution  at  the  "next"  discrete  t,  namely  at  t^  +  2At.  Each 
subsequent  entry  produces  the  solution  at  the  next  t; 


is  an  integer  (or  .name  of  an  integer)  which  specifies  the  number 
of  equations  in  the  system,  i.e.,  n=N+l; 


The  method  requires  a  zeroth  equation 


..  =  ^  =  1 


dt 

which  is  treated  in  the  same  manner  as  the  other  equations . 
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This  provides  for  adjusting  the  value  of  the  independent 
variable  as  required  by  the  method. 


EVDS  is  the  naine  of  the  location  where  the  derivatives  are  explicitly 
defined; 

Similar  to  the  Simpson  integration  of  the  previous  example,  this 
statement  is  designed  for  arbitrary  f^,  hence,  again  it  is  the 
responsibility  of  the  programmer  to  write  the  statements  which 
evaluate  the  f^  of  interest.  These  statements  must  be  designed 
to  record  the  evaluated  derivatives  in  consecutive  places  called 
Y* .  Further,  this  sequence  of  statements  must  be  terminated  vrith 
G0T0(R .K. GD) .  Specifying  the  name  EVDS  in  the  ENTER  statement  and 
terminating  this  evaluation  of  the  derivatives  with  GOTO(R.K.GD) 
enables  the  computer  to  evaluate  and  weight  the  derivatives  as  many 
times  as  required  by  the  method. 


Y  is  the  name  of  the  first  of  the  N+1  functional  values  (or  equiva¬ 
lently,  the  name  of  the  "zeroth"  functional  value); 

To  avoid  confusion  in  (n)  and  (N+1)  and  the  inclusion  of  the  in¬ 
dependent  variable  as  one  of  the  (n+I)  functional  values  and  deri¬ 
vatives  it  is  suggested  that  the  system  be  considered  and  identi¬ 
fied  as  follows: 


dt 


y’ 


f 


1 


dt 


dt 


y' 

“'n 


-r 

X 

n 


l44 


To  identify  and  reserve  space  for  the  recjuired  (luantities  it  is  suggested 
that  the  programmer  use  the  following  or  equivalent  BLOC  statement 

BLOC(Y  -  Y  n  )y'  -  Y'  n)Q  -  Q  n) 

As  stated  above,  Y  is  the  name  of  the  first  of  the  (N+l)  consecutive  functional 
values.  These  correspond  to  the  names  of  the  given  conditions  and  correspond¬ 
ingly  to  the  names  where  the  resulting  values  will  be  recorded. 

Y'  is  the  name  of  the  first  of  the  (N+l)  consecutive  derivative  values; 

It  is  in  these  places  that  the  EVDS  statements  instruct  the  computer 
to  record  the  derivative  values.  As  identified  above,  y'  =  1 
dt 

corresponds  to  ^  =  1. 

Q  is  the  name  of  the  first  of  (N+l)  quantities  that  are  required  for 
intermediate  values  in  the  computations. 

For  subsequent  steps,  i.e.,  for  solutions  at  t^  +  kAt  where  k  =  1,2,5,..., 
etc,,  it  is  not  necessary  to  re-specify  the  parameters  as  in  the  original  EKTEH 
statement.  Indeed,  for  subsequent  steps  one  need  only  write  GOTO(R.K.Gl) . 

GIVEN: 


dt 

=  y’  -  1 

where  y  =  t,  and  at  t  =  t  =  0 

^Yl 

dt 

=  y{  =  (a/h)y2(t) 

y(t.)  =  0 

dy2 

dt 

=  yA  =  (-■b/a)y^(t) 

o 

II 

o 

■p 

rH 

dv_ 

"5 

dt 

=  y^  =  (d)y^(t) 

II 

o 

p 

y3(to)  =  c 

where  a,  b,  c,  and  d  are  constants 


i45 


REQUIRED :  For  a  =  2,  b  =  3»  c  =  d  =  -1  , 

Determine,  print  (and  identify) 

yj(t)  for  t  =  0,  .1,  .2,  ..  ,  1.0  . 

Use  a  computation  step-size  At  =  .01 


The  flow  chart  for  this  example  is  given  in  FIGURE  kO  below;  the  program  and  output  in 
FIGURE  4. 


FLOW  CHART  FOR  EXAMPLE  20. 


FIGURE  i+O. 


EXAMPLE  21. 


This  example  is  given  to  illustrate  addition^  subtraction,  multiplication 
and  inversion  of  matrices.  We  use  ENTER  statements  for  muJ.tiplication  and 
inversion  whereas  we  use  general  formxilas  and  COuifT  statements  for  addition 
and  subtraction.  The  illustrations  given  in  this  example  do  not  include 
permissible  options  such  as:  references  to  symmetric  matrices;  references  to 
the  transpose  of  a  matrix;  references  to  internally  stored  matrices  w’hose 
elements  are  uniformly  spaced  different  from  unity;  or  "accummulative  multi¬ 
plication",  i.e.,  adding  the  matrix  C  to  the  product  matrix  AB  and  recording 
the  result  at  C.  For  ease  in  i-eferenee,  it  is  convenient  to  identify  matrices 
as  two  dimensional  arrays.  This  is  readily  achieved  through  the  use  of  BLOC 
statements . 


GIVEN: 


REQUIRED: 


A,  B  and  C,  where: 

A  is  a  (2x2)  matrix; 

B  is  a  (2  X  2)  matrix; 

C  is  a  (2x3)  matrix. 

Print  and  identify  the  following  matrices: 


1. 

A 

6, 

D  =  B 

X  C 

2. 

B 

T- 

E  =  D 

+  C 

3. 

/-f 

8  • 

H  =  A 

f 

- 

4. 

A-1 

Q. 

✓ 

M 

II 

> 

■1a 

5. 


The  flow  chart  for  this  example  is  given  in  FIGURE  k2. 
The  program,  input  and  output  are  listed  in  FIGURE  43. 
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OETA  =  -30000000  1 

OETA*  =  - 

33333333 

OOOOOOl 

A  MATRIX 

B  MATRIX 

C  MATRIX 

0000002 

1«0000  2.0000 

3.0000 

4.0000 

7.0000 

8.0000 

9.0000 

2.0000  1.0000 

5.0000 

6.0000 

3.0000 

2.0000 

l.OOOO 

A  INVERSE 

A  INV. 

INV. 

D  =  B*(: 

0000003 

-  .3333  .6667 

1.0000 

2.0000 

33.0000 

32.0000 

31.0000 

.6667  -  .3333 

2.0000 

l.OOOO 

53.0000 

52.0000 

51.0000 

lOENTI lY 

NULL 

E  =  D+C 

0000004 

1..0000  -  .0000 

.0000  - 

.0000 

40.0000 

40.0000 

40.0000 

.0000  l.OOOO 

-  .0000  - 

.0000 

56.0000 

54.0000 

52.0000 

FIGURE  k3. 


To  form  a  sxim  of  matrices,  (or  differenc-e),  we  write  a  general  equation 
using  an  index,  then  by  means  of  an  initial  setting  of  the  index  and  an 
appropriate  COUIW  statement,  the  sum  or  difference  is  obtained.  As  illustrated 
on  the  flow  chart  and  the  program,  for  the  sum,  E  =  D  +  C,  we  initially  set 
the  index  J  =  0.  The  general  equation  is  then  expressed  as: 


SUMMAT 


E1,1,J  =  D1,1,J  +  C1,1,J  <jo 


So  that  J  takes  on  the  required  integer  values,  we  write 

COUNT ( 6 )IN(J )G0T0( SUMMAT ) 

Note  that  J  was  set  equal  to  zero  on  line  ^  and  hence  the  first  evaluation  of 
the  general  definition  corresponds  to  simming  the  first  elements  of  the  matrices. 
Note  also  that  we  set  M  =  2  and  N  =  5  on  line  4.  This  was  done  only  to 
illustrate  later  references  to  dimensions  by  name  rather  than  by  explicit  in¬ 
tegers. 

Next,  to  obtain  the  inverse  we  could  have  written  the  appropriate  ENTER 
statement  immediately  following  line  5.  However,  as  in  the  case  of  solving 
linear  equations,  the  given  matrix  is  replaced  with  the  resulting  inverse. 

Hence,  since  we  want  the  original  A  matrix  for  future  operations,  we  move  A  to 
F  and  retain  the  original  at  A.  We  will  instruct  the  computer  to  invert  F 
and  hence  after  the  inversion  is  completed,  A  ^  will  be  recorded  at  F.  To 
obtain  the  inverse  of  F,  'we  write 

ENTER ( MAT . INV )F1 , 1 ) 2 ) DETA 

where : 


MAT. INV  is  an  abbreviation  for  matrix  inversion; 

F1,1  is  the  name  of  the  first  element  of  the  matrix  to  be 

inverted,  ( successive  element  in  consecutive  spaces  by  ro'ws); 


2  indicates  the  dimension,  i.e.,Fisa2x2  matrix; 

DETA  is  optional  and  when  specified  as  above  is  the  name  of  the 
value  of  the  determinant  of  the  matrix  to  be  inverted. 
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22. 


In  the  next  examples,  several  methods  for  determining  approximations  to 
real  roots  of  continuous  functions  are  illustrated.  This  exam.ple  illustrates 
a  "bisection"  method. 

GIVEN: 


f(X)  =  X-  -  X  -  1  ; 

X  and  X,  such  that  f (X  )f(Xj 
o  1  o  1 

e  >0 


0  ; 


REQUIRED: 


Determine  X_,X  ,X,  , .  .,X.,  such  that  |f(X.)|<e  ; 
^  ^  ^  j  J 

Print  and  identify,  and  f(X^)  for  i  =  0,1,2, 


In  the  bisection  method,  the  "next"  approximation  for  the  root  is  defined 


as 


X  = 


X  +  X 
_S _ E 


where  f (X  )  <  0  and  f(X  )  >  0.  The  function  is  evaluated  for  this  mean  X 
'  n  p 

and  a  test  is  applied  to  determine  if  the  magnitude  of  the  function  is  less 
than  the  given  £.  If  the  magnitude  is  less  th.an  e,  X  is  the  desired  root.  If 
the  magnitude  is  not  less  than  e,  then  the  interval  in  which  the  root  lies  is 
diminished  by  replacing  X^  or  X^  by  X.  We  replace  X^  by  X  if  f(X)  <  0,  or 
X  by  X  if  f(X)  >  0.  The  next  mean  X  is  determined  and  the  process  continued. 
The  flow  chart  for  this  example  is  given  in  FIGURE  44, 
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Since  the  process  will  rec^uire  an  indefinite  nuinber  of  evaluations  of 
the  function,  we  wrote  the  general  definition  for  arbitrary  X  and  terminated 
the  general  definition  and  printing  of  X  and  f(X)  with  a  "variable”  exit,  E. 
Prior  to  each  entrance  to  the  evaluation,  we  establish  values  for  the  arbi¬ 
trary  X  and  the  variable  exit,  E. 

The  function  is  first  evaluated  for  X  =  with  exit  E  =  1.  At  1.  we 

record  f  =  f(X  )  and  establish  either  X  or  X  depending  on  whether  f(X  ) 
o  o  p  n  o 

was  positive  or  negative . 

The  function  is  next  evaluated  for  X  =  X^  with  exit  E  =  4.  At  4.  we 

record  f^  =  f (X^  )  and  establish  either  X  or  X  ,  again  depending  on  whether 
11  p  n ' 

f(X^)  was  positive  or  negative. 

At  6.  the  exit  E  is  set  equal  to  J.  for  all  future  exits  from  the  general 
function  evaluation.  Next  we  apply  a  test  to  determine  if  the  initial  con¬ 
dition  is  satisfied;  i.e.,  if  f(X  )f(X  )  <  0.  If  the  condition  is  not  satis- 
fied,  this  indication  is  printed  and  the  computer  is  directed  to  the  next  pro¬ 
blem.  If  the  initial  condition  is  satisfied,  the  initial  X^  and  X^  required 
f'or  the  general  process  are  established  and  recorded.  Hence,  we  direct  the 
computer  to  the  definition  of  the  mean  X  and  initiate  the  general  process. 

The  program  and  results  for  this  example  are  listed  in  FIGURE  45.  Note 
that  this  program  may  be  modified  for  other  f(X)  and  corresponding  initial 
conditions  by  replacing  cards  5  and  15  accordingly. 
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X 

F(  X) 

0000001 

10000000 

I-IOOOOOOO  1 

0000002 

20000000 

1  50000000  1 

0000003 

15000000 

1  87500000 

0000004 

12500000 

1-29687500 

0000005 

13750000 

1  22460938 

0000006 

13125000 

1-51513672-01 

0000007 

13437500 

1  82611084-01 

0000008 

13281250 

1  14575958-01 

0000009 

13203125 

1-18710613-01 

0000010 

13242187 

1-21279454-02 

0000011 

13261719 

1  62088296-02 

0000012 

13251953 

1  20366507-02 

0000013 

13247070 

1-46594883-04 

0000014 

13249512 

1  99479097-03 

0000015 

13248291 

1  47403682-03 

0000016 

13247681 

1  21370716-03 

0000017 

13247375 

1  83552438-04 

0000018 

13247223 

1  18477852-04 

0000019 

13247147 

1-14058747-0.4 

0000020 

13247185 

1  22094948-05 

0000021 

FIGURE  45. 


The  flow  chart  for  this  example  is  given  in  j^'IGukE  46. 


FLOW  CHART  FOR  EXAMPLE  85. 

FIGURE  46. 

The  program  and  results  are  listed  in  FIGURE  47.  Note  that  this  program 
may  be  modified  for  other  f(X)  by  replacing  cards  3  and  5  accordingly.  Note 
alsO;  this  process  may  not  converge,  indeed,  it  may  diverge  if  f'(x)  tends  to 
zero. 
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EXAMPLE  2k. 


This  example  illustrates  the  "Regula -Falsi"  method 
approximation  for  a  real  root  of  a  function,  f(X)  =  0. 
next  approximation  for  the  root  is  defined  by: 


_  \ 
cl  .  J 

X  ,  such  that 

<  0  ; 

b.) 

X  ,  such  that 

P 

f(X  ) 
P 

• 

0 

A 

c.  ) 

f(X  ) 
n 

d. ) 

\ 

H.A  ; 

P 

for  obtaining  an 

n  X-iixto  liic  oiiuu.^  oiit; 


In  general,  the  new  approximation. 


X  = 


f(XJ 


X  f(X  ) 
p  n 

'Wi 


n' 


T 


GIVEN: 

X^  and  X^  such  that  f(X^)  f{X^)  <  0,  f(X)  =  X^  -  X  -  1,  and  e  >  0 

REQUIRED: 

Determine  X  such  that  f(X)  <  e.  Print  and  identify  each  X  and 
corresponding  f(X) 

The  flow  chart  for  this  example  is  given  in  FIGURE  48.  Note  the  similarity 
to  the  bisection  method  flow-chart.  The  program  and  results  are  listed  in 
FIGURE  49. 
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» 


X 

F(  X) 

0000001 

lOOOOOOO 

1-10000000  1 

0000002 

20000000 

1  50000000  1 

0000003 

11666667 

1-57870370 

0000004 

12531120 

1-28536303 

0000005 

12934374 

1-12954209 

0000006 

13112810 

1-56588487-01 

0000007 

13189885 

1-24303747-01 

0000008 

13222827 

1-10361850-01 

0000009 

13236843 

1-44039499-02 

0000010 

13242795 

1-18692584-02 

0000011 

13245320 

1-79295919-03 

0000012 

13246391 

1-33630103-03 

0000013 

13246645 

1-14261375-03 

0000014 

13247038 

1-60474995-04 

0000015 

13247119 

1-25643  798-04 

0000016 

13247154 

1-10873904-04 

0000017 

13247169 

1-46109160-05 

0000018 

FIGURE  h9. 


EXAMPLE  25. 


This  example  illustrates  a  "constant  secant"  method  for  approximating  a 
real  root  of  a  function,  f(X)  =  0.  Again  we  assume  that  and  X^^  are  given 
such  that 

f(X^)  f(X^)  <  0. 

The  slope 

f(X^)  -  f(X^) 

“  "  ~T,  He 

1  o 

is  computed  and  the  "next"  approximation  to  the  root  is  defined  similar  to  the 
Newton -Raphson  method,  i.e., 

f (X  ) 

X.  =  X. - - 

1+1  1  m 

GIVEN: 

f(X)  =  X^  -  X  -  1 

X  and  X  such  that  f(X  )  f(X_,  )  <  0  ;  e  >  0. 

0  1  '  o  '  1'^  ^ 

REQUIRED: 

Determine  X  ,X  ,X,  , . . . ,X . ,  such  that  |f(X.)|  <  e  ; 

Print  and  identify  X^  and  f(X^),  i  <  0,1,2, ... ,j 

Tli0  x*low= ciis.jr't  foi*  "tills  iii0'tltO  d  is  shown  in  FIGURE  50.  The  program  and  output 
obtained  are  listed  in  FIGUEIE  51' 


172 


FLOW  CHART  FOR  EXAMPLE  25- 


FIGURE  50. 


23  JULY, 63  PAGE  1 

PROB  C906  M. J.ROMANELL I  45107  EXAMPLE  25  1 

START  X0=1  %  Xl=2  %  EPS=.OOOOl  2 

PRINT<  X  F{X)>  %  ENTER(PR1.\T  B)%  X  =  XO«  SET(E=1.)  3 

eve  F=X**3~X~1  t  PRINT(X)F  %  GOTO(,E)  4 

1.  F0  =  F%X=XJL?SEr(E  =  2.  )GOTO(  EVE)  5 

2.  IF(F*FO<0)GOTO(3.)  XPRINT<CONDiTIONS  NOT  S  AT  I SF  I  ED>?GOTO  (  N  .  PROB  )  56  6 

3.  M=:(  F-FO)  /  (  X-XO )  %  SFT(E=4.)  7 

4.  IF-ABSlF<EPS)GOTO(N.PROa)?  X=X-F/M*  GOTO(EVF)  8 

END  GOTO(STAKT)  9 


MAY. 23, 63  BRLESC 
PKOB  C906  M.J 


FORAST  F62 

ROMAIMELLI  45107  EXAMPLE  25 


» 


X 


F(X) 


0000001 


10000000 
20000000 
11666667 
12631173 
13044269 
13185779 
13229171 
13241950 
1 3245665 
13246741 
13247053 
13247143 
13247169 


1-10000000  1 

I  5U00GU00  1 

1-57870370 

1-24785752 

1-84906120-01 

1-26035534-01 

1-76669149-02 

1-22292803-02 

1-64576564-03 

1-18685691-03 

1-54051212-04 

1-  15633  700-04 

1-45217501-05 


0000002 

0000003 

0000004 

0000005 

0000006 

0000007 

0000008 

0000009 

0000010 

0000011 

0000012 

0000013 

0000014 


FIGURE  51. 


HiAATOt'litlj  dD. 


This  example  illustrates  another  iterative  method  for 
approximation  to  a  real  root  of  a  function  f(X)  =  0.  Here 
f(X)  =  0  can  be  expressed  in  an  equivalent  form 


determining  an 

4  4-  4/*.  4-1^^4- 

1.0  l.i3  cai3i3LUU<rU.  OilClU 


X  =  F(X) 

The  general  iteration  is  then  expressed  in  the  form 

X.^l  =  F(X^) 


where  we  assume  an  initial  X  =  X^  is  given  to  start  the  iteration. 


GIVEN: 

f(X)  =  X^  -  X  -  1  j 

^o  ^ 

e  ■>  0 


REQUIRED: 

1/5 

Use  the  iteration,  X.  ^  =  (X.  +  l)  to  determine  X.,  such  that 

|f(Xj)  I  <  e. 

Print  and  identify  X^  and  f(X^),  1  =  0,1,2, 

The  flow-chart  for  this  example  is  given  in  FIGURE  52. 
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STAET  EVT’ 


FLOW  CHAkt  i<'Ok  kXAMPLE  26 
FIGURE  52. 


The  program  and  output  obtained  are  listed  in  FIGURE  55* 
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1 

PKOB  C906  M. J.KOMANELLI  45107  EXAMPLE  26  1 


START 

PRIMTC  X  F( X) >^ENTER{PRINT  B)1  X=U  C=l/3* 

EPS=. 00001  2 

lEVF 

F  =  X**3-X-nPRlNT(  X)F*1F-ABS<F<EPS  )GOTO(N.PROB) 

3 

X=|X+1)**C  €  GOTO(EVF) 

4 

END 

GOTO( START) 

5 
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» 

X 

F(  X) 

0000001 

10000000 

1-10000000  1 

0000002 

12599210 

1-25992105 

0000003 

13122938 

1-52372787-01 

0000004 

13223538 

1-10059982-01 

0000005 

13242687 

1-19149254-02 

0000006 

13246326 

1-36388070-03 

0000007 

13247017 

1-69123259-04 

0000008 

13247149 

1-13129931-04 

0000009 

13247174 

1-24939947-05 

OOOOOlO 

FIGURE  53. 


EXAJ.IPLE  27. 


This  example  illustrates  a  "least-squares"  method  for  "fitting"  a 
polynoniial  to  given  tabular  data .  Although  there  are  packaged  programs 
available  for  this  purpose,  this  and  the  next  two  examples  are  given  to 
illustrate  how  the  ENTER  statement  is  used  to  generate  the  normal  equations 
required  for  least-squares  solutions.  To  illustrate  we  consider  the  follow 
ing  problem: 


GIVEN; 

Xf,  y(X,)  for  i  =  1,2,3,  . 

Assume  this  data  on  cards  in  standard  form,  one  pair  per  card. 

Since  an  indefinite  number  of  pairs  are  given,  we  will  terminate 
this  data  with  two  blank  cards. 


REQUIRED; 

Determine,  print  and  identify,  C^,  and  such  that 

S  =  ^  (y(X^)  -  y(X^))^  is  a  minimum, 

i 


where 


y(X.)  =  C,  +  C_X.  +  C,X.^  . 
1  1  2  1  -3  1 


The  necessary  conditions  for  a  minimum  are: 


Ss 


-  n 


= 21  *  21  *  21 


As 


=  0 


'  "  21  *  21 
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yZ 

“1-5 


_z: 


Y~^r(Y  '> 


To  obtain  the  solution,  we  must  generate  and  solve  the  above  system  of 
linear  equations  in  the  three  unknowns,  C^,  C^,  and  C^.  Note  that  the  matrix 
of  coefficients  of  the  unknowns  is  symmetricj  i.e.. 


Note  also  that  the  system  to  be  generated  and  solved  is  of  the  form  illustrated 
in  EXAMPLE  l8,  i.e.,  of  the  form 


AC  =  b 


where: 

A  Is  symmetric 

C  is  the  imknown  vector,  ; 
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Hence,  the  following  elements  must  be  generated  and  recorded  internally  in 
consecutive  spaces  by  rows. 


To  generate  the  above,  the  FORAST  language  includes  a  convenient  EIHER  state¬ 
ment  of  the  form 

ENTER(F.N.E.)Al,l)n)Pl)W  ^ 
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F.N.E.  is  an  abbreviation  for,  "Form  Normal  Equations"  ; 

Al,l  is  the  name  of  the  first  element  in  the  first  row  of  the 

symmetric  system  to  be  generated; 

(space  should  be  reserved  for  this  system  with  a  BLOC 
statement) 

n  is  an  integer  (or  name  of  an  integer)  'which  specifies  the 

number  of  equations  in  the  system  to  be  generated: 

(for  our  example,  n  =  3) 

is  the  name  of  the  first  of  (n+l)  consecutive  elements,  in 
particular. 


General  form. 


dC^ 


(x,,y,) 


P 

n 


n 


P  .,  =  V. 
n+±  ”  1 


For  our  example. 


^1  =  1 


P2  =  Xi 


p, 

5  1. 


^4  = 


i8l 


d.BiiO'tcs  cvslu.S'fci on  of  "bliG  p3-i''txsT  dGi'ivs.'tivG 


at  (X^^y^)  ;  (again^  space  should  he  reserved  for  these  (n+l)  quantities  with 
a  BLOC  statement). 


W  is  optional  and  if  specified  is  the  name  of  a  wei^t  value  to 
be  applied  in  the  sums,  if  unspecified,  the  weight  value  is 

ocixv'^ti.  ciCD  u-LiJ.  oy  • 


Corresponding  to  a  given  point,  (X^,y(X^),  we  instruct  the  computer  to 
compute  and  record  the  (n+l)  quantities  at  P^,  P^,  ...  ,  Then,  the 

above  ENTER  statement  tells  the  computer  to  make  the  necessary  contributions 
to  the  sums  in  the  normal  equations.  When  all  points  have  been  processed  in 
this  manner,  we  instruct  the  computer  to  solve  the  symmetric  system  generated 
through  the  use  of  the  ENTER  statement 


ENTER(SY.SNE)A1,1)5)C1  io 
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rnV\  ^  ^/-vnn+*?/-vv^  •i  ^  +V^/-vv^  n  ^  T  r«  T-\  T  /-k  T  ^  r«  >-n  r  <^1 1  Vv  n<^rt 

Xiicr  av-»_LU.u -LVJii  j.o  c>iic:ix  ctvctj.j.cxuj_c  x  v^x  Jjxxxiuxii^  *-'x  x  wx  cxiij  o  u.uc>cvivxcxi  u  vxo^ 


desired. 


The  flow-chart  for  this  example  is  given  in  FIGURE  5^* 


Clear  the  aiigmented 
symmetric  A  matrix  area 

,•  „  _ U.,TT  _  r\ 

j. .  e  .  t>e  o  aj-x  >>■  ^  =  <j 


Read: 


hr^ 


X*  y 

j  ^  \  one  iiuiiiuei' 

of  quantities 
read^ 


-  Yes 

.1  =  0  ?|  ^  i 


FLOW  CHART  FOR  EXAMPLE  27. 


FIGURE  54. 


The  program,  input  data  and  results  obtained  for  this  example  are  listed  in 


FIGURE  55. 


H 

(XI 
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PROB  C906  M. J.ROMAMELLI  4bl07  EXAMPLE  27  1 

BLOC ( A1 , l-A3,4/SY. ) Pl~P4)Cl-C3)  2 

START  CLEAK(9)iN0S.AT(A1,  1  3 

REAUPT  REAU(X)Y:6  I  F- I  NT  (  09  =  0  )  GO  TO  (  SOL  •  N  )  4 

Pl=l^  P2  =  X?  P3=X*X5:  P4  =  Y^  5 

ENTER(F.ivl.E.  )A1,13  3)Pli  GOTO(KEAOPT)  6 

SOL'N  ENTERCSY.SfMOAl,!)  3)C1«  7 

PRlNr<Cl  =  >C).<  C2  =  >C2<  C3  =  >C3«  GOTO  ( (M .  PROC )  8 

END  GOTO (START)  9 


(j.O  3.0000 

.1  2.9850 

.2  2.9402 

.3  2.8660 

.4  2.7632 

.5  2.6327 

.6  2.4760 

.7  2.2945 

.8  2.0901 

,9  1.8643 

l.O  1.6209 
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Cl  =  30077021  1  C2  =  -10693520  C3  =  -12891375  I  0000001 


FIGUEE  55. 


EXAMPLE  28. 


In  the  previous  example  we  obtained  a  "polynomial"  fit  to  some  tabular 
data.  In  this  example  we  will  include  a  computation  and  printing  of  the 
"residuals”;  i.e.,  the  difference  between  the  given  functional  value  and  the 
approximating  functional  value  at  each  point.  We  will  also  compute  and  print 
the  "root  mean-square"  error. 

First,  since  we  need  the  given  functional  values  to  compute  residuals, 
we  will  have  to  make  these  functional  values  available  after  the  solution  for 
the  approximating  function  is  obtained.  (Recall  that  we  discarded  the  data 
points  as  soon  as  their  contribution  to  the  normal  equations  was  completed). 
Hence,  to  make  the  given  data  available  for  residual  computations,  we  will 
first  read  and  record  all  of  the  given  data  and  retain  it  internally  for  as 
many  future  references  as  desired. 

To  provide  for  a  maximvim  of  say  500  points,  we  will  reserve  space  for 
1000  values  of  given  data,  hence  we  will  use  a  BLOC  statement 

BL0C(D1  -  DIOOO). 

Further,  since  only  two  data  values  are  recorded  per  card,  we  will  specify  a 
format  which  departs  from  standard  and  indicates  the  desired  departure. 

The  READ  statement  to  accomplish  the  reading  and  recording  of  a  maximum 
of  1000  values,  two  per  card  has  the  form 

READ-F0RMAT(f/2)-(1000)K0S .AT(D1 ) 
where  the  format  F  which  specifies  the  particular  form  desired  is: 

F  FORM  (10-12)10-12)  2  ^ 


185 


To  identify  and  refer  to 
of  the  form 


oh-B  glVGi 


's  and 


wc  wx  . 


Q  CVT\T 
d  0x11 


D  OCX  OCUiCii  O 


rrtnvr/ vi  t\o\  H 

=  l)1  j  ll  =  UC.)  ^ 

Hence^  to  refer  to  the  given  data  we  may  refer  to  the  X's,  Y's  or  D's  bearing 
in  mind  that  the  X's  are  two  units  apart  in  internal  stcarage,  and  likewise, 
the  Y's  are  two  units  apart  in  internal  storage. 

The  flow  chart  for  this  example  is  given  in  FIGURE  56. 
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FLOW  CHAET  FOR  EXAMF’LE  28. 
FIGURE  56. 
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PROB  C906  ri.  JI.ROMANELLI  45107  EXAMPLE  28  1 

BLOC ( A1 , 1-A3,4/SY. ) P 1-P4 ) C 1-C 3 ) D l-D 1000  2 

SYN  ( Xl  =  J1 ) Yl  =  U2 ) 2J  =  09  )  3 

START  CLEAH(9)U0S.AT(A1,1)^  RE AO-FORMAT ( F/2 ) - ( 1000 ) NOS . AT ( D1 ) 4  SET(I=0)  4 

NEXTPT  Pl=r^  P2  =  X1.,1[:«  P3-=P2*P2%  P4  =  Yi,l  5 

ENTEK(F.N.E.  )A1,1)3)P1)?:  COUNT  (  2  J/2  )  I N  (  I  )  GOTO  (  NE  X  TP  T  )  6 

ENTEi<(SY.SNr)Al,l)3)Ci)i;  7 

PRINT-FORMAr(FO)-<Cl  =  >C1<  C2  =  >C2<  C3  =  >C3«  ENTER(PRINT  B)  8 

PRINT<  X  Y  YBAR  R>7  ENTER(PRI.NT  B)  9 

SET (1=0)^  S=0  10 

RESID  YBAR-Cl  +  Xl ,  I  (C2+C3*X1 ,  I  ):«  R=YEAR-Y1,U  S  =  S  +  R*R*  11 

PRINT-FORMAT  (  FU  )  -  (  X  L  ,  1)  Y 1  ,  I  )  YBAR)  RSf  COUNT  (  2J/2  )  IN(  I  )  GOTO  (  RES  ID  )  12 

ENTERICVI TOF )09)N?  RMS=S CRT ( 2 *S/ N ) X  ENTERIPRINT  B)  13 

PkINIT-FOKMAir(FO)-<RKS  =  >RMSX  GOTOtN.PROB)  14 

F  FURM( 10-12) 10-12 )2  15 

IFO  F0RM(  12-4-10)3-2)  1-4)  2  16 

END  GOTO(START)  17 

O.O  3.0000 

.1  2.9850 

.2  2.9402 

.3  2.8660 

.4  2.7632 

.5  2.6327 

.6  2.4760 

.7  2.2945 

.8  2.0901 

.9  1.8648 

1.0  1,620:> 


i 
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3.0077 

C2  = 

.1069 

C3  = 

1.2891 

X 

Y 

Y8AR 

R 

.0000 

3.0000 

3.0077 

.0077 

.1000 

2.9850 

2.9841 

- 

.0009 

.2000 

2.9402 

2.9347 

- 

.0055 

.3000 

2.8660 

2.8596 

- 

.0064 

.4000 

2.7632 

2.7587 

- 

.0045 

.5000 

2.6327 

2.6320 

- 

.0007 

.6000 

2.4760 

2.4795 

.0035 

.7000 

2.2945 

2.3012 

.0067 

.8000 

2.0901 

2.0971 

.0070 

.9000 

1.8648 

1.8673 

.0025 

1.0000 

1.6209 

1.6116 

— 

.0093 

RMS  ==  .0056 


OOOOOOl 


FIGURE  57. 


EXAMPLE  29. 


In  the  previous  example  of  "least-squares  curve -fitting",  the  approxi¬ 
mating  function,  (the  polynomial),  was  linear  in  the  unknown  coefficients. 
This  example  illustrates  "least -squares  curve -fitting"  where  the  approxi¬ 
mating  function  is  non-linear  in  the  unknowns.  The  method  illustrated  is 
often  referred  to  as  the  method  of  "differential  corrections".  To  illustrate 
we  consider  the  following  problem: 

GIVM; 

X,,  y(xj,  (i  =  1,2,5,...)  ; 

X  X 

(Here  we  will  assume  the  given  data  is  recorded  on  cards  in 

standard  form  with  six  values  per  card),  ,  X^,  y^,  etc. 

A  and  B  ,  initial  estimates  of  the  unknowns  A  and  B  in  the 
o  o 

approximating  function 

REQUIRED: 

Determine,  print  and  identify,  A  and  B  such  that 

S  =  \  (y(X^)  -  y(X^))^  is  a  minimum. 


where 

y(X)  =  A(B“)  . 

Compute,  print  and  identify  residuals  and  RIvlS  as.  in  the  previous 
example . 

The  necessary  conditions  required  to  minimize  S  with  respect  to  A  and  B 
would  lead  to  a  system  of  equations  which  is  non-linear  in  the  iinknowns  A  and 
B.  Hence,  we  first  expand  y(X,.A,B)  in  a  Taylor  series  about  A^.B^.  We  obtain 
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then 


y(X,A^  +  M,  +  AB);5J  y(X,A^,B^)  +  ^  ^ 

o 


AB  =  y(X,AA,AB) 


an  approximating  function  which  is  linear  in  the  unknowns,  M  and  AB.  (The 
notation  I  denotes  evaluation  of  the  partial  derivative  at  A  ,B  ).  Note 

|o  o"  O'' 

that  we  have  substituted  A^  +  AA  for  A,  and  B^  +  AB  for  B.  The  necessary 
conditions  for  minimizing 


v(X.)  -  y(X.) 
"  '  1'  1' 


with  respect  to  M  and  /SB  are: 


as 


=  0  =  2 


Z[  ]i-l. 


_  n  -  p 


Z[ 


or  equivalently. 


Svl  \  ,  'n]  ,by 

5b  J  Z_  '5b, 


)  AB  = 


U  'U 


z. 


“ /,ir  A  \  \ ^y 

y^A.,a^,..o^^5B 
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o  j  o  ucm  xii  iiict  ux  liuuduxuii  ux  ox^t;  xui  m 


CD  =  E 


where: 


C  = 


2 

'v  ,?h7 

r\\r 

) 

> 

) 

o 

N.  “ 

i 

2  V  -N  - 

o 

2 

o 

% 

)  >  (% 

) 

UI5 

o 

y  u£> 

o  — 

■> 

o 

1 

; 


the  symmetric  matrix  of  coefficients; 


D  = 


M 

AB 


the  unknown  vector; 


E  = 


(y(x,)  -  y(x.,A^,B^))f 


>  (y(Xj  -  y(X,,A^,Bj)^| 

. 


KJ  yy 


.  .  .  ,  „  ^  _ _ _  _ _  ^  ^  -1  ^  h  ^  ^  ^  ^  X  X-»  U  'IV  U  I  ^  U  *1^  {  |j '  TVT  Ij '  ^  X  X'  'r^  X  lo  j-N 

cib  ill  oiitf  jjitfvxuut)  ^  wtr  ud.ii  ut)t:  uuc  iDi'ixiAi\\^x'  •I’l.iii.  /  uu  x  ux  iii  oiiu 


x^u  ^ 


normal  equations  by  defining: 


X. 

=  B  ^ 

n 


IQP 

•*-.✓*- 


O 


o 


X. 


=  A  X.B 
o  1  o 


=  A  X.P,/b 
o  1  1'  o 


P_  =  vfX.  )  -  vfX.  .A  .B  ■)  =  y(X.  )  -  A  P. 
X  */  \  V  '  o'  i'  o  1 


After  forming  the  normal  equations  we  can  use  the  EIWER ( SY . SNE )  to  obtain  the 
corrections ;  and  We  then  form 


A=  A  +  M 
1  o 

B,  =  B  +  AB 
1  o 

which  serve  as  "new"  estimates  for  A  and  B  for  repeating  the  entire  process. 

To  terminate  this  correction  process  we  will  test  two  conditions: 

(l.)  if  both  (aA.  j  <  e  and  |AB|<e  we  will  consider  the  problem  solved 
and  direct  the  computer  to  the  computation  and  printing  of  the 
residuals; 

(2.)  if  after  10  corrections,  (l.)  above  is  not  satisfied,  we  will  assume 
that  the  process  is  diverging  and  no  residuals  will  be  computed. 

In  either  case,  (l.)  or  (2.),  we  will  repeat  the  process  for  different  initial 

estimates,  A  and  B  . 

o  o 

The  flow  chart  for  this  example  is  given  in  FIGURE  58* 

The  program  and  results  obtained  are  listed  in  FIGURE  59* 
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FLOW  CHART  FOR  EXAJ4PLE  29- 
FIGURE  53. 


VO 

av 
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PROS  C906  M. J.ROMAMELLI  45107  EXAMPLE  29 
BLOC  I A1 , 1-A2, 3/SY. ) P1-P3)C1-C2) Dl-DlOOO) 

SYN  (X1=D1)Y1=D2) 

START  RtAOC  iOOO)NOS.AT(Dl)SE  2J=09 

NEXTC  READ(AO)BO)«  PR  I NT-FORMAT  (  F  ) -<  AO  =  >A0<  BO  =  >BOSi  ENTER(PRINT  B) 

A  =  A04;  d  =  BG4:  SET(K  =  0) 

NEXTIT  SET(I=0)i  C LE AR ( 5 ) NOS . AT ( A 1 ,  1  ) 

EPOS  Pi=B**Xl,I?  P2=A*XI,I*P1/B*  P3=Y1,I-A*P1 

ENTER  (F.N.E.  )  Al,  l)2)Pn  COUNT(2J/2)  IiM(  I  )  GOTO  (EPOS) 
ENTER(SY.SNE)A1, 1)2)C1*  A=A+C1«  B=B4C2 

PRINr-FORMAT( F)-<A  =  >A<  B  =  >B 

IF-ABS(C1<.001  )  ANL)-A0S(C2<.001  )GOTO(  RCAP  ) 

COUNT (10) IN(K) GOTO (NEXTIT) 

PBAC  ENTER(PRINT  B)*  GOTO(NEXTC) 

RCAP  SET(I=0)«  S=Q%  eNTER(PRlNT  B)* 

PRlNr<  X  Y  YBAR  R>%  E:NTER(PRINT  B) 

RESIO  Y6AR=A*B**X1 , I*  R=YRAR-Y1,I*  S=S+R*R 

PRINT-FORMAT ( F)-( XI , DYl , I ) YBAR)R 
C0UNT(2J/2) IN( I )GOTO(RESID)*  ENTER ( C V  I TOF ) 2 J ) N 


F 

RMS=SCRT  (  2*S/N)  PRINT 

FORM! 12-4-10) 3-2)  1-4)2 

END  GOTO (START) 

-FORMAT( F )-<RMS 

=  >  RMS« 

GOTO! P0AC ) 

0.0 

3.005  0.5 

4.223 

1.0 

5.98  3 

1.5 

8.5  2.0 

11.99 

2.5 

17. 

-1.0 

• 

! 

C 

• 

2.180 

1 

• 

o 

.747 

1 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 


4.0 

l.U 


1.5 

1.0 
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AO 

=  4.0000 

BO  = 

1.5000 

A 

=  3.0115 

B  = 

1.9496 

A 

2.9963 

R  = 

2.0025 

A 

2.9972 

8  = 

2.0016 

X 

Y 

YBAR 

R 

.0000 

3.0050 

2.9972  - 

.0078 

.5000 

4.2230 

4.2404 

.0174 

l.OOOO 

5.98  30 

5.9993 

.0163 

1.5000 

8.5000 

8.4878  - 

.0122 

2.0000 

11.9900 

12.0085 

.0185 

2.5000 

17.0000 

16.9895  - 

.0105 

— 

l.OOOO 

1.4700 

1.4974 

.0274 

— 

.5000 

2.  1800 

2.1185  - 

.0615 

— 

2.0000 

.  7470 

.  7481 

.0011 

RMS 

=  .0253 

AO 

=  1.0000 

BO  = 

1.0000 

A 

=  4.6050 

8  = 

4.4132 

A 

=  1.7647 

B  = 

3.8917 

A 

=  2.0481 

B  ••= 

2.5747 

A 

*  2.8476 

H  = 

1.9399 

A 

2.9959 

0  = 

2.0066 

A 

=  2.9972 

B  = 

2.0017 

A 

=  2.9972 

B  = 

2.0016 

X 

Y 

YBAR 

R 

.0000 

3.0050 

2.9972  - 

.0078 

.  5000 

4.22  30 

4.2404 

.0174 

l.OOOO 

5.9830 

5.9993 

.0163 

1.5000 

8.5000 

8.4878  - 

.0122 

2.0000 

11.9900 

12.0085 

.0185 

2.5000 

17.0000 

16.9895  - 

.0105 

— 

l.OOOO 

1.4700 

1.4974 

.0274 

— 

.5000 

2.1800 

2.1185  - 

.0615 

- 

2.0000 

.7470 

.7481 

.0011 

RMS  =  .0253 


0000004 


0000005 


FIGURE  59 
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Some  problems  require  random  numbers  to  simulate  errors  or  other  random 


nru-^r^  ^  ^  tTlU  ^ 

j.iixi3  c^a.iii_^xc:  xxx u.^  ux  <a  oc: c>  cl  pc>cuu.o -x  cLtiu.um  iiuiiiucx  ^ciicxcloljx  •  xiic 


numbers  are  called  "pseudo-random"  since  their  genesis  is  known  and  any  set 
produced  may  easily  be  reproduced.  Further,  for  large  enough  samples,  they 


■'rMT^rorrh  "i  r^nci  1  t^v'/^'Kq  cj>^H 

V-.X1.  X>>.4>  _W  i.  \JJ  <_>.  X  X\-*. 
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use  in  various  problems.  In  this  example,  we  illustrate  a  generator  of 
normally  distributed  pseudo-random  numbers.  Two  ENTER  statements  are  used  for 
ceneratinc  these  numbers;  thew  have  the  <7eneral  form 


ENTER ( NRNOSl ) A1 ) n ) B1 
ENTER ( NRN0S2 )A1 ) n )B1 


NRNOS  is  an  abbreviation  for  Normal  Random  Numbers; 

The  two  distinct  entrances,  NRNOSl  and  i©N0S2,  are 
essential  for  reproducibility.  The  initial  entry, 

NRNOSl,  always  yields  the  same  set  of  numbers.  Sub¬ 
sequent  entries  at  NRN0S2  yield  sets  which  differ  from 
preceding  sets. 

Al  is  the  name  of  the  first  of  n  standard  deviation  values; 

It  is  assumed  that  they  are  recorded  in  consecutive 
locations  in  internal  storage .  For  standard  deviations 
CT^,  i  =  1,2, ...,n,  the  pseudo-random  numbers  produced,  X^, 
lie  in  the  interval 


-1+CT.  <  X.  <  4a. 

11  1 

n  is  an  Integer  (or  name  of  an  integer)  which  defines  the 

nunibcr’  or  psoudo-r&ndoDi  numtiors  dosirsdj 

is  the  name  of  the  first  of  the  n  pseudo -random  numbers 

cTCtr\  o  T'Q  "h 


1  nft 

j-yyj 


B1 


In  general^  "the  programmei’  should  establish  the  values  in  prior  to  the 
entry  to  the  generator  and  should  reserve  space  for  these  and  the  resulting 
values  which  will  be  recorded  in  B. ,  i  =  l,2,...,n. 

In  the  example  which  follows,  we  will  generate  one  pseudo-random  number 
with  each  entry  to  NRWOS.  Further,  we  will  specify  a  standard  deviation  a  =  1. 

To  illustrate,  we  will  generate  sets  of  pseudo-random  numbers  with  set  (sample) 
size  equal  to  N,  where  N  =  200,  ^4-00,  600,  800,  1000  ,  2000,  ...  ,  5000* 
Corresponding  to  each  set  N,  we  will  determine,  print  and  identify  the  following: 

MN  the  number  of  negative  numbers  in  the  setj 

WP  the  number  of  positive  numbers  in  the  set; 

LX  the  largest  number  in  the  set; 

SX  the  smallest  nimiber  in  the  set; 

MO  the  number  of  X^  in  the  range,  -k  <  £  -2.9999  j 

Ml  the  number  of  X^  in  the  range,  -2.9999  <  <  -1.9999  ; 

M2  the  number  of  X^  in  the  range,  -1.9999  <  “  *9999  ; 

M5  the  number  of  X^  in  the  range,  -  .9999  <  X^  <  0  ; 

PO  the  number  of  X^  in  the  range,  0  <  X^  <  .9999  > 

PI  the  number  of  X^.  in  the  range,  .9999  £  1.9999  ; 

P2  the  number  of  X^  in  the  range,  1.9999  £  2.9999  t 

P5  the  n’omber  of  in  the  range,  2.9999  <  <  2.9999  ; 

51  the  percentage  of  N  which  lie  within  1  standard  deviation  of  the  mean; 

52  the  percentage  of  N  which  lie  between  1  and  2  standard  deviation  of  the  mean; 

53  the  percentage  of  N  which  lie  between  2  and  3  standard  deviation  of  the  mean. 

The  flow  chart  for  this  exanqile  is  given  in  FIGURE  60. 

The  program,  input  and  output  are  listed  in  FIGURE  6l. 
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PKOB  C906  K. J.ROMANELLI  45107  EXAMPLE  30  1 

BLOC  (M-MiO)  *  2 

SYN  {P=M4)Sl=M8)  %  3 

START  PRINir<  N  NN  NP  LX  SX  M3  M2  Mi  M0>  4 

CONT  <  PO  Pi  P2  P3  Si  S2  S3>*  EImTER(PRINT  B)«  5 

SET(iN  =  200)Ni\l  =  0)NP=0)  J  =  Oi  SX=LX=0:«  D=1  6 

LNTEk  (iMKNOSl)  0)  1)X  %  CLEAR(83  NOS.  AT  (M)  %  7 

B1  IF  (X<0)  GO  TO  (B2)  INT  (NP=NP+1)  t  IF  (X>LX)  GO  TO  t  B4 )  8 

B6  ENTER  (CVFTOI)X)  I  P,I=P,1  +  L  %  GO  TO  (AATJ)  9 

82  InT  (NN=inN+1)  %  IF  (X<SX)  GO  TO  (B3)  10 

B5  X  =  -X^  ENTER  (CVFTODX)!)  *  I  NT  {  K=  3- 1  )  ^  M ,  K=M,  K  +  l*  GO  TO  (AATJ)  11 

B3  SX=X  %  GO  TO  (85)  12 

84  LX=X  X  GO  TO  (B6)  13 

AATJ  COUNT(N)  IN(J)GOTO(MOKNOS)«  ENrEK(CVITOF)N)N'*  SET ( I =0 ) GOTO ( T A8 )  ife  14 

MORNOS  ENTER  ( l^l<^NO  S2 )  0 )  1  )  X  X  GO  TO  (Bl)  15 

TAB  JimT(K  =  3-I)«  S  1  ,  I  =  I  00  (  M  ,  K +  P  ,  I  )  /  N  '  «  COUNT  (  3  )  I  N  1  I  )  GOTO  1  T  AEU 'S  16 

PRLmT-FOxMAT  (F)-  (N)NN)NP)LX)SX)  (  IDNOS.ATIM)!?  17 

COUNT  (1200/200)  IN  (N)  GO  TO  (MORNOS)  18 

COUNT  (5001/300)  IN  (N)  GO  TO  (MORNOS)?  OOTO  (  Nl .  PROP  )  ?  19 

F  FORM (4-6) 1-3)3-1)12-1-5)1-2)3-1)1 1-3-4 )  1-3 ) 3-1 ) 1 1-4-5 ) 1-2  )  3- 1 ) 11-3-4 1  20 

CONTl-2) 3-1) 1-1) 11-2-3) 1-4)?  21 

END  GOTO (START)  22 
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CONCLUSIONS 


Specific  examples  have  been  employed  as  a  vehicle  to  enable  the  novice  to 
program  for  the  ORDVAC  and  BRLESC  Scientific  Computers  and  to  obtain  appropri¬ 
ate  solutions.  While  the  full  generality  of  the  FORAST  language  has  not  been 
presented,  a  sufficient  varity  of  approaches  are  made  available  to  indicate 


the  flexibility  and  application  of  this  method.  As  sophistication  is 
developed,  both  mathematically  and  in  programming  ability,  the  student  should 


X  tx  CX  »-/<-»  -I-  • 

■»  m 


e  complete  discussion  of  FORAST.  Additions  to  the 
language  subsequent  to  the  publishing  of  reference  j^l^  are  found  in  Appendix 
A  oi  this  report.  Detailed  information  on  available  subroutines  for  plotting 
results  can  be  found  in  reference  [a.]  . 
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APPENDIX  A. 


SOME  ADDITIONS  TO  THE  FORAST  LANGUAGE 


The  following  additions  to  the  FORAST  language  are  now  available . 


1.)  GOTO, I  (A)b)c) . i 


where  contents  of  I  must  be  an  integer  1,2,3)  . i- 


The  transfer  of  control  is  directed  to 

It  It  It  It  It  tt  It 

It  tl  II  tt  tt  tt  It 

A,B,C,  etc.,  may  themselves  be  indexable  names. 
2.)  GOTO, (I  +  inc)(A)B)c)  .  ^ 


A  if  I  =  1 
B  if  I  =  2 

V 

C  if  I  =  3)  etc. 


Same  as  above  with  provision  for  incrementing  the  existing  value  of  I  by 
a  positive  or  negative  increment. 


NOTE;  These  additions  reduce  the  size  of  the  ORDVAC  SYN  Table  from  64  to  55* 
Compiler  still  checks  for  synonym  full  at  64. 


3.)  ENTER  (PLOT)r)X^)XM)Y^)YM)h)ix)iy) 
where 

X  =  Address  of  1st  X 
o 

XM  =  Address  of  last  X 

Y  =  Address  of  1st  Y 
o 

YM  =  Address  of  last  Y 


h  =  Handler  number 


ix 

iy 


distance  between  X  entries) 
distance  between  Y  entries) 


Taken  as  1  if  not  entered 


If  r  >  0,  tape  is  rewound  before  return  to  program  is  made. 
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This  subroutine  is  designed  to  be  used  when  a  "quick  look"  at  some  data 
is  desired.  It  is  not  necessary  to  determine  beforehand  the  scales^  maximum 
and  minimum  values^  etc.,  since  the  subroutine  scans  the  data  and  computes  the 
necessary  values. 


The  subroutine  produces  (on  tape)  the  information  for  a  26  by  26  inch  * 
plot  of  the  data  and  a  plot  of  the  quantities  XMAX,  XMEN,  XSCALE,  YMAX,  YMIN, 


ySCALE  as  ^  digit  (+  or  -) 


integers  wj.th  2  digit  (+  or  -)  10' s  exponents. 


WARNING;  Processing  data  for  plotting  uses  a  lot  of  tape.  If  there  are 
many  large  variations  in  Y  throughout  the  plot,  one  entry  to  this  subroutine 
may  easily  use  one -half  of  a  reel  of  tape. 

4.)  ENTER ( SET. TI)u)E.T  )b)BMAX  ^ 


optional  after  first  entrance  for  this  tape  unit. 

u  is  tape  unit  integer  (eff.  address  itself  is  used.) 

l<u<5  or  9<'r<l4 

E.T.  is  optional;  if  it  is  zero  (or  blank)  then  the  routine  goes 
to  N.PROB  when  the  END  TAPE  sentinel  is  read. 

If  specified  (not  zero),  then  the  routine  jumps  to  that 
address  when  the  END  TAPE  sentinel  is  read. 

B  is  the  initial  address  of  a  block  of  core  storage  that  is 

large  enough  to  hold  the  largest  block  on  the  tape  being  read. 

BI''IAX  is  the  last  address  in  the  storage  block  for  "tS-pG  • 


SKT.TI  allows  a  program  to  read  data  on  magnetic  tape.  It  sets  the  computer 
so  that  subsequent  READ  statements  (or  A. READ  or  READBL  subroutines)  will  cause 


data 


the  tape  unit  specified.  Each  80  characters  on  tape  is 


considered  to  be  a  "card"  by  this  routine.  (Easier  use  of  "formated"  hi-speed 
printer  tape's  might  be  allowed  in  the  future.)  As  many  as  six  input  tapes  may 
be  used  in  one  program  by  entering  this  subroutine  at  different  times  in  the 
program,  the  data  will  continue  with  the  "card"  that  follows  the  last  "card" 


*  This  may  be  changed  to  26  x  13  by  means  of  a  manual  switch  at  the  plotter  board. 
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that  was  read  from  that  tape.  Each  unit  should  have  its  own  block  if  the 
program  ever  re-uses  that  unit  because  a  part  of  a  block  may  need  to  stay 
there  while  another  unit  is  being  used.  When  entering  with  a  unit  that  was 
prevlousl;/  used,  it  is  not  necessary  to  specify  the  storage  block  addresses; 
if  specified,  they  will  be  ignored.  It  is  not  possible  to  change  the  storage 
block  once  it  has  been  already  assigned.  The  storage  block  may  be  longer  than 


any  block  on  the  tape  but  must  not  be  shorter  than  the  longest  block  that  is 
read  from  the  tape.  (A  '*card'*  requires  8  words  of  memory  storage.  The 
storage  addresses  may  be  larger  than  OkOOO  in  the  large  memory.)  The  tape 
block  length  can  be  variable  and  if  any  block  is  longer  than  the  storage 
allocated,  the  rest  of  the  block  will  be  ignored.  All  tape  reading  is  parity 
checked  and  re-read  five  times  before  causing  the  erroneous  "card"  to  be 
punched  and  a  RUN  ERRCR  card  saying  "PAR.ERRORu" . 


The  "E.T."  (end  tape)  address  should  be  zero  unless  it  is  actually  needed. 
If  it  is  zero,  this  tape  unit  is  rewound,  the  computer  is  set  to  read  cards  and 
control  goes  to  N.PROB.  If  an  address  is  specified,  then  these  things  should 
be  done  by  the  program  before  going  to  N.PROB. 


It  is  desirable  that  a  standard  end  of  tape  sentinel  be  used  by  everyone. 

It  is  also  nice  to  have  a  standard  end  of  reel  sentinel.  This  routine  uses 
"ENDbTAPEbb"  (b  is  blank)  as  the  end  of  tape  sentinel  when  it  appears  as  the 
first  ten  characters  at  the  beginning  of  a  block  and  the  next  ten  characters  do 
not  say  "ENDbREELbb" .  When  the  next  ten  characters  do  say  "ENDbREELbb " ,  then 
it  assiomes  that  there  is  another  reel  to  be  read  on  this  same  unit,  so  it  re¬ 
winds  the  tape  and  halts  at  08l  so  that  the  operator  can  mount  the  new  reel. 

(The  unit  no.  is  in  the  B  address  of  the  halt  Oidei'.  )  Standard  BRLESC  output 
will  have  the  END  TAPE  sentinel  if  the  "rewind  tape  8"  switch  was  properly  used. 
When  making  tapes  off-line,  an  extra  block  of  one  card  with  this  sentinel  should 
be  added  at  the  end  of  all  the  data . 
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wiicii  i  ca.u.±ii^  ucipe  ^  one  oiii  cc  iieciu.ex'  u^x'u.t)  oiici  o  wei  e  pi'uuueeu.  ±ii  ix'uiio 

of  previous  FORAST  or  FORTRAN  output  are  automatically  skipped.  (It  checks 
for  "bbBRLESCbb"  in  characters  11-20.)  A  dictionary  will  not  be  automatically 


oVt  "^11+  ■?  + 

0.11.  V-4.  ^  ^ 


c  Q  1  i.TCnrc?  "Pr^l  ~\  r\i.7a/^  "Kir  ■P/^n>»  "KT  Q  y^lr  r»c3  o  Viov^no  +V>i=i  A  T?TT*  A 

uO  <.U._I.intO>jrO  .L'o.^.i.'ortov-*.  J.  \./«.4.X  »  •  ■«..■■  r  %  txo  U4.11tX  llN.^1101^  1^110  >  ■  ■'-«  ■  • 


routine  may  be  used  and  your  program  can  check  for  the  blank  cards  (you  need 
only  check  the  first  word)  before  starting  to  read  actual  data. 

A  word  within  this  routine  is  named  SKP.TL  and  it  may  be  used  to  "skip 
tape  lines".  If  it  is  set  to  an  integer  (not  fl.pt.),  then  the  next  tape 
read  will  skip  that  many  "cards".  (if  the  skip  Includes  the  "header  cards", 
then  they  must  be  included  in  the  integer  that  is  put  into  SKP.TL). 


Th,G  inpu.'t  ds-tis,  insiy  s.l’tGrn.s.'tG  bG't.WGGn  t.Q,p0  3,nd  cs-pcis  s.t  will .  Th.G  usg 
of  ENTER ( SET . Cl )^  will  set  the  computer  for  reading  cards.  (SET. Cl  is  a  small 
subroutine  within  the  SET.TI  subroutine.  If  tape  No.  6  is  being  used  instead 
of  card  Input,  then  it  sets  for  tape  6  input.  If  cards  or  tape  6  is  being 
used  at  the  time  SET. Cl  is  entered,  it  does  nothing.) 

There  are  three  possible  error  prints  in  this  subroutine.  They  are: 

SET.TI  6  Tried  to  use  more  than  six  tape  units. 

SET.TI-BLTj  Have  negative  buffer  length.  (B  >  BMAX)(or  length  >16,585) 

PAR.ERRORu  Parity  error  on  tape  unit  u.  Is  preceded  by  the  "card" 

that  contains  the  error. 


It  is  permissible  to  set  for  tape  unit  u  when  the  same  unit  u  is  being 
used  at  the  time  SET.TI  is  entered. 

5.)  goto(c.prob)^ 


n  T>Ti  /^"n  j 


.  J.T_  - 


o.rxwjiD  ±b  biiu  nariiti  ui  a  buui'uuoxiit?  uiiau  a±xuwfc>  uiit;  e:uiiip±xa uxuii  aiiu. 

running  of  several  programs  without  stopping  between  them.  It  is  similar 
to  N.PROB  and  does  everything  that  N.PROB  normally  does  except  for  putting 


t~\  4  T  /^i  1  4*  1  4*  4*  ./-s  n  ^  v»  j-A  4  4*  ^ 

ct  X  XJLC  iiictx  JL  uii  oa,_^c;  v-*  cxiivo.  V.-JVXli^  XWX  XCWX1XU.X11^  '-J  • 


f 

\  X-  '-'X 


purposes  of  tape  6  input,  tape  8  output  and  for  operator  control  of  the 
computer,  several  programs  combined  by  using  C.PROB  will  still  be  con¬ 


sidered  as  one  problem. ) 


210 


Either  ETflER  or  GOTO  (or  any  jimp  order)  may  be  used  to  enter  C.PROB. 

Any  new  program  that  is  compiled  after  going  to  C.PROB  must  have 
a  PROS  card  as  its  first  card.  It  is  permissible  to  allow  a  READ  state¬ 
ment  to  read  this  PROB  card  as  the  sentinel  to  indicate  that  the 
previous  problem  is  done  r’onning  if  the  name  C.PROB  was  used  somewhere 
in  the  program.  (if  C.PROB  is  in  the  dictionary^  then  reading  a  PROB 
card  causes  control  to  go  to  C.PROB  subroutine  and  the  PROB  card  is 
also  used  as  the  first  card  of  the  next  program  to  be  compiled.) 


In  any  set  of  programs,  the  last  one  should  not  use  C.PROB.  (if 
this  is  done  and  tape  6  input  is  used,  then  the  next  problem  will  be  done 
as  a  continuation  of  your  problem. ) 


The  "date"  will  be  propagated  through  all  the  programs  if  it  precedes 
the  first  one. 


When  using  C.PROB,  the  "card  counter"  (067)  and  the  input-output 
options  will  remain  set  to  what  they  were  at  the  completion  of  the 
previous  program.  (The  permanent  constant  block  (p^O-OTL)  is  not  re¬ 
read  before  the  next  program  is  compiled.  However  ERROR  (066)  and  M.DUMP 
(058)  are  reset  to  N.PRQB  ) 

After  a  RUN  ERROR  print,  control  goes  to  N.PROB  (not  C.PROB)  unless 
the  program  has  used  ERROR  as  a  location. 

If  tape  input  has  been  used  (SET.TI  subroutine),  C.PROB  will  set 
for  "card  input"  before  compiling  the  next  program. 

6.)  ENTER(MAX.  )A)b)c) . io 


MAX  is  a  subroutine  that  finds  the  largest  floating  point  number  and 
stores  it  in  the  last  address  specified  in  the  ENTER  statement. 
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MIN.  Is  a  subroutine  that  finds  the  smallest  floating  point  number 
and  stores  it  in  the  last  address  specified. 

8. )  ENTEE(MAX.I)i)j)k) . ii 

MAX. I  is  a  subroutine  that  finds  the  largest  integer  and  stores  it 
in  the  last  address  specified. 

9. )  enter(min.i)i)j)k) . ^ 


MIN. I  is  a  subroutine  that  finds  the  smallest  integer  and  stores  it 
in  the  last  address  specified. 

In  each  of  the  above  four  subroutines  (actually  four  entrances  to  one 
subroutine),  the  number  of  arguments  is  variable  and  is  determined  by  the 
number  of  addresses  written  in  the  ENTER  statement.  (There  must  be  at 
least  two  arguments.)  The  last  address  is  always  the  store  address  and 
is  not  used  as  an  argument.  Each  argument  is  a  single  number,  not  a 
block  of  numbers . 


Care  must  be  exercised  in  using  negative  integers  as  arguments  for 
MAX. I  or  MIN. I  because  a  "negative  integer"  in  an  index  register  will 
appear  as  a  large  positive  integer.  (The  comparison  is  made  on  a  full 
word  basis . ) 

The  next  order  following  one  of  the  subroutine  entrances  must  not  be 
a  NOP.  (it  won't  be  unless  you  write  one  as  an  assembly  order.) 

10.)  biwer(mod. )a)b)c  io 

This  subroutine  will  compute  C  =  A(mod  B)  where  A,B,  and  C  are  the 
respective  addresses  of  floating  point  numbers  as  indicated  in  the  ENTER 
statement  above.  It  is  defined  as  C  =  A  -  WH0LE(A/b)*B  in  terms  of  a 
FORAST  formula. 
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the  elements  of  the  rows  and  colrunns  of  the  matrices  to  be  stores  in 
equally  spaced  memory  positions. 


ENTER(MAT.MP)A^^)B^^)C^^)i)j)k)z)r^)c^)r^)c^)r^)c^)% 


where 


^11’®11^^11 


are  the  addresses  of  the  first  elements 
of  the  respective  matrices. 


Scxmc  cLc>  cLuvjvc 


i,  j,  k  are  the  dimensions  of  the  matrices 

(j  is  always  the  commop  dimension. ) 

z  is  a  3  digit  sexadecimal  number  (Od^d^d^) 

where  indicates  the  options  applying 

to  matrix  d^  to  matrix  B,  and  d^  to 

matrix  C. 


d.  =  0  means  to  use  the  matrix  as  stored, 
1. 

not  augmented 

d^  =  1  means  to  use  the  transpose  of  the 
matrix  stored. 


d^  =  2  means  the  matrix  is  augmented 


d^  =  3  means  both  1  and  2  apply 

if  d^  =  4,  5z  or  7  it  means  to  accumulate 
in  C  as  well  as  above  options . 


r  is  the  distance  between  the  first  elements  of  each 
a 

row  of  matrix  A. 


(ie)  address  Ap^ 


11 
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c  is  the  dista.nce  between  the  first  elementb  uf 
a 

each  colimn  of  matrix  A. 


(ie)  address  -  address  A^^ 


^b'^b  ^  same  meaning  for  matrices  B  and  C. 

Pairs  of  arguments  may  be  omitted  from  the  right  ^ 

Zero  arguments  for  r  and  c  are  not  valid,  vmen 
r  and  c  is  specified,  the  augment  portion  of  z 
for  that  matrix  is  ignored. 


EXAMPLE: 


Given: 

^11 

\l 

^11 

^12 

\2 

c,  „ 
±'^ 

a,  , 
i;) 

b,  , 

-l;) 

"i3 

^21 

V 

=21 

^22 

^22 

=22 

“23 

“23 

‘^23 

^31 

^31 

""31 

^ZO 

^32 

""32 

^53 

"33 

EM'ER(MAT.MP)a^^)b^^)c^j^)5)3)5)010l+)l5)5)l5)5)l5)3)5^ 


X  X  X  X 

X  X  X  X 


X  X  X  X 


would  produce 


A  B  +  C  in  C 


i.e.,  c^^  -  c^^  +  +  agl  ^21  ^31^31  ' 
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